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DYADIC PROGRAMS AND SUBDUAL METHODS*t 


A. Charnes 


Northwestern Technological Institute 
and Purdue University 


W. W. Cooper and Merton Miller 


Carnegie Institute of Technology 


INTRODUCTION 

The past few years has witnessed an increasing amount of research and more efficient 
methods for solving important types of linear programming problems. Guided in part by 
reports and experiences in practical applications, some of this work has been directed toward 
evolving suitable transformations and interpretations so that problems from one area might 
thereby be made amenable to efficient methods of attack which were originally designed for 
problems in other areas. E. H. Bowman [4], G. T. Bishop [2], A. S. Manne [28] and others! 
have, for example, shown how features of "transportation methods" might be adapted to prob- 
lems of production scheduling. Other efforts —e.g., W. Prager [31], L. R. Ford, Jr., and D.R. 
Fulkerson [22}° — have been devoted to improving and extending these methods, while simul- 
taneously showing how additional problem areas may then be brought within the ambit of these 
improved versions of standard methods. Even problems which appeared to be nonlinear in 
nature — e.g., because of probabilistic elements *— have thus been found amenable to this mode 
of attack. Finally, efforts have been made to utilize some of these features in a more general 
framework and some attention has been devoted to research in the joint use of dual and direct 
methods of attack which might lead to a "pincers" approach* for improved general computa- 
tions, or at least as a means for ascertaining bounds to the solutions in early stages of a large 
calculation. The "'primal-dual algorithm" of G. B. Dantzig, L. R. Ford, Jr., and D. R. 
Fulkerson [26a] and E. M. L. Beale's "method of leading variables" [1]° both lie athwart this 
line of attack. 





*Manuscript received 14 May 1958. 

Released as ONR Research Memorandum No. 21, December 1957. Part of the research 
underlying this paper was undertaken for the project, Methodological Aspects of Management 
Research, at Purdue University; and part for the project, Planning and Control of Industrial 
Operations, at Carnegie Institute of Technology. Both projects are under contract with the 
Office of Naval Research. Reproduction of this paper in whole or in part is permitted for any 
purpose of the United States Government—Contract N-onr-1100-(05), Project NR 047-016, and 
Contract N-orn-760-(01), Project NR 047011. 

lSee [9]. 

2The "flooding techniques" of A. W. Boldyreff [3] should also be mentioned because of their 
close resemblance to the Ford-Fulkerson approach. Also, while less rigorous and complete, 
Boldyreff's methods show an origin in management-type problems and thus lie more clearly in 
the kind of evolution which is here being described. 

3Vide, G. B. Dantzig and A. R. Ferguson [21]. 

a ve discussion of this possibility in A. Charnes and C. E. Lemke [6] and [27]. 
Cf [34]. 
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Each of these methods has added further power to potential programming applications 
by either (a) providing ways to identify and relate the structural properties in different prob- 
lems so that they can thereby be made amenable to attack by one or more of the available 
algorithms or (b) adding new methods of solution which are suited to problems with special 
structural features. Considered collectively, this work suggests the possibility of one or 
more general approaches which may be considered, in some respects, an alternative to the 
use of general purpose methods (such as the simplex and dual methods). 

For example, one such alternative might evolve into an "algorithm for generating 
algorithms" over certain classes of problems. Such an algorithm would take the form of a 
set of rules or guiding principles for (a) identifying and transforming problems (or parts 
thereof) into one or more standard types of models and (b) synthesizing solution methods 
which take advantage of the special structural properties which are thus revealed. 

The dyadic transformations and the subdual methods presented in this paper are 
intended as a start in one possible approach to the search for such algorithms. It builds on 
past experience in that each of the methods for solving linear programming problems (which 
were cited in the opening paragraph) can be generated by utilizing the principles which will be 
outlined in this paper. As will be shown by an example, it is also possible to utilize this 
approach to effect further extensions of some of these algorithms to situations which were not 
comprehended in the original formulation. 

On the other hand, the objective ''an algorithm for generating algorithms" is not 
achieved. It will become evident during the discussion that the methods for effecting trans- 
formations and generating specialized algorithms suited (more or less) to particular problems 
is less automatic than might be desired.© To compensate for this deficiency, observations 
will be :nade on the kinds of situations where the indicated approach might usefully be con- 
sidered. But these observations are not as sharply pointed as should ultimately become pos- 
sible (it is hoped) when more experience is gained in identifying significant model types.’ 

The topics to be treated in this paper can be summarized as follows: First, a means 
for transforming every linear programming problem into a dyadic equivalent is supplied. 
Second, a general characterization of subdual methods is provided which are, essentially, 
elaborations of the dual method of C. E. Lemke in that the direct problem constraints are suc- 
cessively tightened at each succeeding stage of the calculations. Third, the usual requirements 
on the rowrank — or basis requirements — associated with the dual method is relaxed in order 
to ensure that the possible range of solution methods is not encumbered from this quarter. 
This is followed by a discussion of "directed subdual methods" (a special instance of the more 
general subdual methods) along with a general means for handling both upper and lower bounds 
on any or all variables. Finally, the preceding theory is illustrated (as much as a single exam- 
ple will allow) for the reasons that have already been indicated. 


DYADIC FORMULATION OF LINEAR PROGRAMMING PROBLEMS 
Figure 1 shows schematically the kind of structure typically encountered in a "blending 
model."' The nutmix example [12], the models for blending aviation gasolines [13], and still 





6The question of a "best" or "most efficient'' model and algorithm for any particular prob- 
lem is not treated in this paper. 
TIncluding approximations to such types. See [9] for further discussion. 
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The following expressions may be regarded 
as a detail drawn from one of the rectangles in 
such a model. Together with its associated ele- 


ments from the functional, this detail may be writ- Figure 1 - Schematic for the structure 
of a simple blending model 


oe 
~ 


ten in the form of a standard linear programming 
problem as: 


max Cy Aq + CoAg + CQ rg 
subject to 


M1444 + AyQAQ + 24g Ag = ay 
(1) 
By inserting new (double subscripts) variables joined by the condition AI = roi , for each 


rj = 1, 2, 3, the model (1) can be written: 


max Cy4q1 + Corgg + Cgrqg + CyAQy + CgAQe + €3%o3 
subject to 
814441 + 22% 72 + 243443 = a 
491491 * 29QAQ9 + 293493 = Ag 
(2) Ay -Agy *5 
Ay2 “Ago sil 
A143 -Aog = 0 
all A.. =0 


This is a dyadic equivalent of (1). Since the constants, c. and aii, in these two problems are 
the same, it is clear that solutions to one are also solutions to the other with only a scale dif- 
ference in the value for the functional. 

In this case a reduction to the echelon-diagnosed arrangement which is characteristic 
of a transportation model? is secured. But this cannot always be expected, as is apparent 





8As explained in [9], the concept of model types is not confined to exact replication of the entire 

model, It allows for (a) model approximations and (b)the synthesis of different models from com- 

binations of standard types. See [7] for a case in which blending patterns are embedded in a 

larger model for simulating (via an optimizing principle) the flowof traffic over a street network. 
See [9] for further discussion. 
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when an attempt is made to apply this same kind of transformation to an extended arrangement 
for a model of the kind shown in Figure 1. 


It is useful, now, to show the general dyadic form which can be written as follows:!° 


m non 
max )) ) “ij ij 
i=1 j=l 
subject to 
(3) (i) + me air = 4 
j=1 
m 
(ii) ); bi Xj = b; 
i=1 


as iL. sk. = &. 
(iv) (a) Some at with the same i are equal. 
(b) Some ij with the same j are equal. 
Here, without loss of generality, the problem is written in a form which provides access to the 
bounded variables techniques discussed in [16]. That is, the Li and U;; are lower and upper 


bounds for each variable, with the values Li, = 0, Ui = © representing the special case 
A., 20. The standard statement of the general linear programming problem, 


1) 
n 
max | oe 
j=l 
subject to 


(4) 


is comprehended in the following specialization of the general dyadic form (3): 





10phis formulation provides access to economies associated with the "bounded variables" 
techniques [16]. 
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max 


subject to 


6) E ai jij = Fi 


(As in (1) and (2) the only alteration is in the value, m » Cj rij - 
j 

is evidently also a specialization of (4), it follows that the two are equivalent in the sense that 

they yield the same solutions. 


for the functional.) Since (3) 


The kind of arrangement that might be expediently em- 
ployed in calculations for a dyadic formulation is presented in Cj | bij Uij 
Figure 2. This arrangement consists of a typical cell which is 
assigned to each rij: The row coefficients, aij? and column 
coefficients, bis, for each such variable are drawn from the Qij 
constraint pairs (ii) and (iii) of (3) and entered in the cell so that 
the corresponding row and column sums may be calculated and 
summed ina convenient manner from some Suitable tabular ar- 
rangements of such cells. The lower and upper bounds, Lip Ui» 
are also entered in the cell for ready reference and so is the 



































Figure 2 - Celldetail in 
a dyadic arrangement 
criterion element, Ci) which appears out of the way, in the upper 


left-hand corner. 

Nonlinear problems!! can also be handled in this kind of arrangement by devices such 
as the introduction of new cells and special methods of record keeping. Piecewise linear func- 
tionals which can be expressed in bounded variables form can be dealt with in a more com- 
pressed manner. Subdivision of the Cij . Lj , and U;, may be employed within the cells of 
Figure 2 and values for the variables recorded therein may be distinguished according to their 
role, relative to an active basis, by suitable notational devices. !¢ 


STRATEGIC CONSIDERATIONS IN DYADIC FORMULATIONS 

The following comments may be helpful in identifying situations where a dyadic formu- 
lation may prove helpful. In general, advantages may be expected to accrue when the dyadic or 
semidyadic!* counterparts afford special insights into the model. This will usually occur only 
when the original problem can be transformed into one with a few constraints applicable to each 
variable. Two general constraints on each variable plus upper and lower bounds on each is an 
ideal case. 





1lsee [9], [10], and [11]. 

l2This kind of compression is employed, for example, by G. B. Dantzig and A. Ferguson 
[21] in their method for routing aircraft under conditions of uncertainty. 

Cases in which transformations to dyadic form are not carried through to completion as 

when, for example, known properties ofcertain constraints makes it inadvisable to alter them. 
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Figure 2 shows that an advantage in producing specialized algorithms can be secured 
from the fact that the "virtual" dual variables!+ "w"— where w' = ¢' B -, B for basis, is 
composed of the collection of vectors for the direct problem { P|iel} which are actively in the 
current basis, and ¢ is the mx1 vector of corresponding functional components — are split 
into two groups: a row group and a column group. In many cases of interest these two groups 
of variables may interact in a simple way to forma Zi = z=. Cij xij 
uate the cells for each Aj e-8-, in a manner akin to the one used for generalized distribution 
(or transportation) models as in [9]. 

In "management-size problems,'' substantial advantages may sometimes be secured 
from the "folding over" that occurs when, for example, access to a transportation type equiv- 
alent can be secured. The compression resulting from this folding over is convenient both for 
the study of structure and the execution of calculations. The dyadic form "essentially" uses 
this same kind of folding over. Note, for example, that the proof offered for the equivalence 
between (3) and (4) made no essential use of the column constraints (3.ii). 

The number of variables encompassed in (3) is roughly mn for m+n general con- 
straints plus 2mn individual bounds. Strategic considerations thus turn on whether an incentive 
is apparent for introducing additional variables which individually participate in fewer con- 
straints than the (possibly fewer) variables in the original formulation. Some of the examples 
cited in the introductory section of this paper are instructive instances. All network analyses 
which are reducible to representation by "incidence matrices"!° are cases in point. This 
includes not only the standard versions (e.g., electrical networks) but also the generalized 
incidence matrices required for accounting networks 16 and spread-sheet analysis as well.!? 
Finally, multicopy networks [7] as well as their dynamic and game equivalents are also pointed 
in this direction. 


which can be used to eval- 


Because the dyadic constraints are equations, their corresponding dual variables are 
(a priori) unrestricted in sign. Where most of the c4;'8 are zero (or can, equivalently, be 
reduced to zero), it may be advantageous to replace the usual "row-column number calcula- 
tion," 18 Z;; = aij Ri + bij K. by its difference counterpart, Zi = aij Ri - bij K, . This device, 
employed but without the ''w"' recognition by Prager [30], for example, makes it possible to 
obtain a ready comparison between the aij Ri and the bij K. when evaluating program possi- 
bilities. Of course, the difference method is less advantageous when most ci; 0. 

In preparation for the next section, it may finally be observed that the theory under- 
lying C. E. Lemke's dual method, [17] and [27] may often provide a useful approach for gen- 
erating algorithms in such situations because (a) the "rectangular" array of cells formed from 
Figure 2 may afford a natural directional ordering of the constraints and (b) the relatively sim- 
ple interaction between "row numbers" and "column numbers" makes it attractive to look 
toward these dual variables as an avenue of attack. 


SUBDUAL METHODS 


It may be recalled that Lemke's general dual method [27] proceeds, iteratively, as 
follows. At each stage a solution to the dual problem is secured along with a corresponding 





14The virtual variables become the actual dual variables at an optimum. 
15 
See [10]. 
See [14 
17See 


25], first edition, pp. 389 and 162. 


18See [9]. 
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direct program. When the latter does not satisfy the initially stated constraints (including non- 
negativity), new dual solutions and corresponding direct programs are arranged until, even- 
tually, both sets of constraints are satisfied. When this occurs, optimum solutions to both 
problems are at hand.!9 

As Lemke has shown, this method may be regarded as operating on extreme points of 
the dual problem.“ Consequently, it may also be viewed (because the dual constraints are 
satisfied) as also operating at each stage with w vectors which reflect optimality for less 
restricted direct problems. 

The relation between these w values to direct problem extreme points and other 
"critical" points has been elaborated by Charnes and Lemke.*! At each succeeding stage the 
w values are optimal for a direct problem which is more restricted than the proceeding one. 
A particular case of this "'tightening"' occurs when each succeeding w implies that the values 
dj = 0 secured at the preceding stage remain nonnegative in the next solution which is 
arranged. 

All such special instances of procedures — which make essential use of the dual method 
and which also arrange to satisfy previously satisfied direct constraints — will be called ''sub- 
dual methods."' The following sections are devoted to showing in more detail where efficient 
use of subduals might be devised, with special reference to problems in dyadic form. 


CONSTRUCTION OF SUBDUAL ALGORITHMS 

To secure increased generality, it is desirable to relax the usual basis requirements 
when the dual (and simplex) method is employed.“* For this purpose, consider the standard 
version Eq. (5) of the linear programming problem, this time in matrix form: 


(6) max c'A with PA = y, and A = 0, 


with P an mxn matrix and ), c each nx 1 vectors which are primed when presented in 
transposed form.*> Its dual is 


(7) min +, with p'P=c' 
and yp otherwise unrestricted. Suppose that a basic dual solution (B, w') is available. That 
is, the solution consists of a selection B of column vectors in P which form a basis for the 


collection of vectors (Pj ‘ -) together with an m <1 vector of values w' satisfying 


(8) w'P=c'*4 





19 See Figure 2 and the related discussion in [27]. 

The dual method therefore offers considerable flexibility for exploring alterations in 
original stipulations, addition or deletion of constraints, etc., matters which are often of con- 
siderable interest in applied work. See [14] and [15] for further discussion. 

elsce [17] and [27]. 

22In particular, it is desired to eliminate the need for a basis in a space with dimensions 
equal to the number of rows in the original matrix. 

23To simplify the presentation, it is assumed that nondegeneracy is present. Little is 
thereby lost since available perturbation devices may be employed (a) to reduce degenerate 
cases to this one and (b) to resolve computational difficulties. See [5], [17], and [27]. 

24Note that there may be an infinite number of choices for a, e.g., the situation where P 
is an incidence matrix. See [10]. 

























A. CHARNES, W. W. COOPER, AND M. MILLER 


and determined from the relation w'B = ¢', where ¢' consists of c;'s associated with 
B. 


Since B is a basis for the vectors (P; ‘ P,); there is a unique solution to the direct 
problem which may be written 


(9) Bn= P 


where (possibly), however, the requirement 7 = 0 may not be satisfied. 
It will now be shown that 


(10) w P; =%, 
is a solution to a dual problem which corresponds to an optimal solution of a (possibly) less 
restricted direct problem. This may be done, moreover, without necessarily assuming that 
the rows of P are linearly independent, or that B is square. 

By definition, 


(11.1) zj = bs CX 





iclI 

when 

(11.2) P, = Zz. Px,;- 
ic! 


It is assumed that B's columns are composed of vectors {P; |icI} where i ranges over one of 
the collections from the set I in P which can be used for a basis. Now let 


(12.1) P. =)° ea 


where e . is the unit column vector with nonnull entry in row s. Then, also, 





m m 
(12.2) P; = >»; ); €545i Xj = 7. e, 7. 455% 
icI \s=1 s=1 icl 
Since the vectors es, $= 1, 2, ..., m, are linearly independent, the coefficients asi « 
in Eq. (12.1) are unique so that 


(13) »: 855 %j = agi 
icl 
By definition, w is selected to satisfy 


(14) Yr weag; = cj» iel. 
s 
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(z ssi) Xij 


z We (2 ‘a =>. Ws j- 


icl - 


Since, by Eq. (10), z; = w'P; the dual conditions w'P =c' really require z, - cj =0 and the 
condition for optimality maintains in a direct tableau. 

It follows that with each such basic dual solution, w, there is associated uniquely an 
optimal basic solution, 7, to a (possibly) less restricted direct problem. When, also, 7 =0 
all of the direct constraints are satisfied and an optimum solution to the entire direct problem 
is secured since the maximum for a less restricted problem cannot be lower than a more 
restricted one (in the same variables). The assertion made in the penultimate paragraph of 
the preceding section is thus verified. *° 

The characteristics of P which make for convenience in recourse to subdual methods 
may now be considered with reference to the special instance of a "directed" case. A directed 
subdual method commences with a basic dual solution (p(t), w)) with an associated value 
nf for the direct problem. If n) is not optimal for the original problem — i.e., some 
nw < 0—a new program (B 2 ,Ww 2 ) is eonnene with an associated 7 2 having nonnega- 
tive entries in the same components as 7 ) plus a nonnegative entry in at least one additional 
component. 

Clearly, any procedure which preserves nonnegativity in this fashion must yield optimal 
direct and dual solutions in at most m iterations if the rank of P is m. The efficiency of this 
procedure is apparent. It is also of interest because most of the methods previously referred 
to make implicit or explicit use of one or more phases of these directed subdual procedures. 

Synthesis of such algorithms will be easiest for matrices, P, with structures which 
can be arranged in echelons of submatrices that (1) encompass the full matrix and (2) increase 
in size to yield a roughly block triangular array in the nonzero entries of P. Dyadic formula- 
tions suggest possibilities of this character which are not evident from the standard form. 


BOUNDED VARIABLES 

Simplification is possible by removing constraints involving upper and lower bounds in 
Eq. (3) from consideration in an "active" basis. The following way of extending the usual 
"Bs - c;" criterion for this purpose is intended to accomplish this in a manner which (a) encom- 
passes both upper and lower bounds on the variables and (b) provides a route for ready incor- 
poration in such special algorithms as may be devised. 

Consider once more the direct problem, 


(16) max c'd with PA = P,, A =U, -=-L, 





*5For a discussion of the use of such subdual methods in establishing bounds (from above 
and below) in pincer-like fashion, see [9]. 
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where I is the nxn identity matrix and U and L are, respectively, vectors which consist of 
the upper and lower bounds in (3.iii). No restriction is imposed on the vector A other than 
those shown in (16). 

The dual to (16) is 


min p'P, + T'U-p'L 





(17) subject to 
p' Pp + 7° wai B ' - c' 
Ty B =0 
and p otherwise unrestricted. 
Given any basic solution, 
(18) w B=c', 


it is evident from (17) that values 7 and 8 may be selected so that a basic solution to 
w'P+7'-£'=c' is readily arranged. Their column vectors being the negatives of each other 


(and hence linearly dependent), the condition 7 B, = 0 must clearly be fulfilled in any basic 
solution. Moreover, 


7; 0< w'Pi<c¢ and tye O<a> fe 24, ) 
(19) 

B; > 0 <=> P; Cj and Bj =0 << w' P; =f, 
(where ''<=>"' means "if and only if''), as is evident from w' P +7' -8'=c' 


The theorem of the alternative in linear programming asserts@® that strict inequality 
cannot, at an optimum, simultaneously obtain for (i) a dual variable and the corresponding 
direct constraint or (ii) a direct variable and its associated dual constraint. Therefore, at an 
optimum the result 7; >0O implies dj = U; and By > 0 implies AY = L,. Because each con- 
straint associated with a lower or upper bound involves only a single variable, its value is 


uniquely determined (at an optimum stage) when the corresponding dual variable is positive. 
Let 


J ={j | Tj m 0} 
(20) 


K = {k|p, > 0}. 


Then 7 (for the direct problem) may again be uniquely associated with w and B via 


(21) Bn = P, - >. PU, - 2: P, Ly. , 


jeJ keK 





26vide, [33]. 
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Now to show how the promised extension of the ''z. - c."’ criterion may be used to 
eliminate the bounds from an active basis (and thereby compress the tableau), define 


(22) Zz, = ». c,X,, Where P, = iB P.x,,- 


icl icl 
Then, again, 
(23) Zp = w' Py 
and the criterion for optimality becomes 


Zz = C, for L, < A, <U, 


(23) Zz <¢, for rj = U; 


Z. > Cy for Ay * Ly). 


Hence, by noting only these conditions, supplemented possibly by simple memorandum 
entries,*’ it is possible to ascertain the status of the variables in the extra dimensions 
associated with the bounds. If N= A; satisfies Eq. (23) and L, =A, = U; for all ieI then 
an optimum has been obtained. If, on the other hand, some dj violates these bounds, a new 
value of w is secured in a manner which preserves such conditions as are already satisfied 
and adds at least one more to the collection. 


AN ILLUSTRATION 

An example will now be given which shows how these procedures may be employed ina 
specific case. A directed subdual method, originally devised by W. Prager [31] in connection 
with warehousing problems, will be extended by means of the preceding theory. In particular, 
this theory will be used to generate a directed subdual algorithm to handle financial (liquidity) 
constraints which were not considered by Prager.*® Interpretative aspects of this problem 
(especially the dual variables) have been dealt with elsewhere [14]. Only computational 
aspects will be examined herein. 

Table 1 contains the data and Figure 3 contains the model for this illustration. Table 2 
shows the optimum direct variable solution (reached in nearly a dozen simplex iterations) and 
Table 3 contains the values for the dual variables which were, of course, available as a by- 
product of these (simplex) calculations. 





27See [9]. 
28See [31]. 
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TABLE 1 
Data for Warehousing Model With 
Financial Constraints 











Period Costs, $c; Prices, %; 
1 25 20 
2 25 35 
3 25 30 
4 35 25 
5 45 50 

















B = 200 tons, warehouse limit 
A = 100 tons, initial inventory 
M, = $1500, initial cash position 
M = $1000, minimum cash position 
t = 1 period lag in cash receipts 























Direct 
Variables Buying Selling 
Stipulations 
oom Xy Xp %3 %q%5) ¥y Yo ¥g Yq V5 
Variables 
Buying u, 1 -1 B - A= 100 
Constraints Uy 11 4 4 B - A= 100 
Us » £ -1 -1 -1 B - A= 100 
Uy : % & 4 -1 -1 -1 -1 B - A= 100 
Us 2 2° ss -1 -1 -1 -1 =-1 B - A= 100 
Selling Wy 1 A= 100 
Constraints Wo -1 1 1 A= 100 
W3 -1 -1 1 1 1 A= 100 
Wy -1 -1 -1 1 1 1 1 A= 100 
Ws -1 -1 -1 -1 1 1 1 1 1 A= 100 
Financial vy 25 M, - M= 500 
Constraints Vo 25 25 -20 M, - M = 500 
V3 25 25 25 -20 -35 M,, - M= 500 
V4 25 25 25 35 -20 -35 -30 M, - M= 500 
V5 25 25 25 35 45 | -20 -35 -30 -25 M, - M=500 
Market data 25 25 25 35 45 20 35 30 25 50 














Figure 3 - 5-period warehouse model with financial constraints 
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Direct Program Solutions: Warehouse Model 
With Financial Constraints—Direct Program 















































































Buy Sell 
Cash 
Period Receipts on 
Quantit Cc: Outla Quanti : 
. J , Y Pj Accrue Com | = 
0 100 -- -- 1500 
1 20 25 500 20 1000 
2 25 120 35 4,200 1000 
3 168 25 4200 30 4,200 1000 
4 35 25 1000 
5 45 168 50 8,400 1000 
6 -- -- -- -- -- 8,400 9400 
Total 288 4700 288 12,600 12,600 
TABLE 3 
Dual Values: Warehouse Model With Financial Constraints 
Buying Selling Financial 
Period Constraints Constraints Constraints 
(u) (w) (v) 
1 u, = Wi= v= 4/5 
2 Uy = Wo = 20 Vo = 
3 Ug = Ws = V3 = 4/7 
4 _— on 
5 
Total 0 70 9/5 




















Easy access to a subdual algorithm is secured through the following transformation. 
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so that g. represents goods remaining at the end of period j and h. the amount available for 
sale in period j+l. Hence Y, +8, = ho = A, the initial inventory and every 0 <h, =< B, so that 
the goods available cannot exceed the warehouse limit. The latter constraints are thus placed 
in bounded variable form. 

Similar interpretations apply to R. and 35 , the funds component, viz., Ry = M and R. 
and S541 are, respectively, the funds remaining at the end of period j and available for j+l. 
Finally, the objective in the transformed problem becomes max a3 - M where, in this case, 
Sn+l ” 56 : 

The cells of Table 4 are now arranged for computation in accordance with the detail 
suggested by Figure 2.*? The variables defined by Eq. (24) are inserted in the boxes (as shown) 
to show how the direct constraints are accommodated. Row 1 is to be interpreted as 
¥, +8,= A= 100, row 2 as C,X, + Ry = 25 X, + Ry = 500, etc. Column 1 conforms to the sum 
PY, + Ry - 35 = 20 y, + R - 35 = 0, Column 2 to 8, +X, - hy = 0, etc. The numerical mag- 
nitudes inserted in the squares of each cell are to be rogarded as multiplying the variables in 
each cell only for the direction (row or column) in which the square is oriented — this is to be 


interpreted as the same direction in which the corresponding equation is to be. When no square 


appears in the direction (row or column) being examined, the coefficient is always unity. Thus, 
Row 3 forms the equation Yo + 89 - hy = 0 while Column 3 is 35 Yo + Ro - 3g = 0. Each row 
and column sums to zero (except the first two rows)°2 so that rim entries are not required. 

All variables are nonnegative so that their negatives are nonpositive. These conditions 
may be derived by substituting the relations Eq. (24) for those shown in Figure 3 (including 
nonnegativity) and effecting the indicated cancellations and eliminations. It will then be found, 
also, that h; = B= 200 for every i= 1, 2,...,5 so that -h; = - B, as shown in Table 4. As 
explained in the section entitled ''Bounded Variables,'"' these variables can then be handled by 
the bounded variables techniques by means of the ''extended Zi; - Cj criterion" in a way which 
will become apparent when the dual variables are discussed. 

Before turning to the dual, it might be noted that no cells appear for 85 Or Xz. Witha 
5-period horizon, it is known that the values for these variables must be zero at an optimum 
(at least so long as prices are positive) so that they might as well be eliminated from consid- 
eration at the outset. Finally, the coefficient, -1, in the cell for -Fe is placed in the corner. 
The objective is to maximize Fs (in the transformed problem).?! Hence, this coefficient is 
drawn from the functional and accorded a different status than the other numerical magnitudes 
(which are drawn from the constraints).°¢ 

The dual to this problem is shown in the lower left corner of Table 4. The variables 
are entered (with altered sign where convenient) along the rims of Table 4 for ease in solving 
these dual constraints. Thus, o7 2 Py¥y> Le., o,- 20 v4 = 0, occurs where these two vari- 
ables intersect at a cell with the value P= 20 oriented in the column direction directly below 
the variable, Yy> which it multiplies. Recalling that a unit coefficient is assumed for direc- 


tions (row or column) in which no rectangle enclosed amount appears, it is seen that o, - 6, 20 


1 1 





<I vide [14] from which this table was drawn, for a corresponding network analogue. 
Note. These may also be interpreted as yy, + 8) - A=0 and 25 x) + Ry - 500 = 0, 
respectively. 


More exactly the objective is to maximize Fe - M but the constant may be eliminated 
because it does not affect the solution. 
32Cf£ (24). 
33Cf remarks on use of "difference counterparts" for "row" and "column numbers" in the 
section entitled ''Subdual Methods." 
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TABLE 4 


Arrangement for a Subdual Method” 
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subject to 
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33 4 


Note: 


Ao, +M6, + B(u, + Hy + Ug + U4) 


re 
_ 
IV 


1 
_ 
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*Table 1 is the source of data for this table. 
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states the next dual constraint. This is followed by 5) -%y 2 0 and C15, - >, = 0, viz., 

25 5y - 2 0, in that order. Notice also that at the "So cell,"’ the dual constraint ¥17- 592 0 
is covered by the expedient of associating the minus sign for the direct variable in that cell 
with each of its dual counterparts. 

A natural directional ordering ** is available for solving the dual as will now be shown. 
Commence with the tail at the lower right and end with the head at the upper left. Notice that 
56 = 0 by the conditions of optimality so that Y5 2 1. Since minimization is the objective set, 
“¥,* -1, the amount in the rectangle for cell 3g which not only reduces this variable to its 
lowest permissible level but also makes it possible to lower 5 and 55 to levels O5 = 50 Y= 50, 
55 "%* 1. From the latter expression it is found that y 4 can be lowered to unity also. That 
is, the lowering of these variables makes it possible to lower their predecessors (working 
toward the upper left). Thus, moving backwards both a natural beginning and order of solution 
is arranged, since setting each dual variable at its lowest level does not increase its prede- 
cessors in the constraint chain shown at the bottom of Table 4. 

The value - %4= -1 makes it possible to set 54 = 1 and O4= 25 which then satisfies 
the dual constraints 64 2 y 4? % = 4 %4= 25 y 4° Notice now that the highest value assignable 
to 4 is then the lower of P4 54 = 35 and O4= 25. Setting -O4= -25 it is now found that 
Hg = % - co) 4" 50 - 25 = 25. By reasoning as in the concluding portions of section 6 the value 
of h 4 is therefore at its bound —i.e., the value -h 4= -200 is implied for this direct constraint 
because the corresponding dual variable, LU 4? iS positive. 

At each place where equality holds in a dual constraint, the vector corresponding to the 
direct variable is entered in the basis. Each such case is indicated by the mark " Y shown 
in the cell which corresponds to each such vector. Where the direct variable is not bounded, 
e€.g-, X4, and its dual constraint holds as a strict inequality this (basis) mark is absent. (The 
direct variable must be zero by the theorem of the alternative.) It is not necessary to exhibit 
the values for Hy explicitly in the bounded variables case.°> When the dual constraints appli- 
cable to such a cell are satisfied as equations, this is treated just as any other variable and 
entered in the basis. (Cf the cell for hg-) When strict inequality holds this variable is set at 
its upper bound. Its vector is not counted as part of an active basis and the symbol b, is 
entered in its cell. (Cf the cell for h 4 -) 

Continuing in this manner all of the dual constraints are satisfied with the values shown 
in Table 4. The value for the functional corresponding to this dual solution is 


(25) 20 x 100 + 1 x 500 + 200 (15 +5 + 0+ 25) = $11,500. 


The implications of these w components for the direct problem are indicated by the " Y and 
% notations in the cells, as already explained. 

The "Stage 0" results being available, the next step is to see whether the direct con- 
straints can be satisfied under the conditions implied by this dual solution. If satisfaction is 
available, an optimum is at hand. If the direct constraints are not all satisfied, a directed 
subdual procedure may be employed. A new dual solution is then arranged which leaves previ- 
ously satisfied direct constraints unmolested while making it possible to satisfy at least one more. 





34Cf discussion in the section entitled "Strategic Considerations in Dyadic Formulations." 
The equation My =%%441] - $4 derives from the constraint My + O54 - 7542 2 0 plus the natu- 
ral directional ordering and the minimizing objective. 
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To start the Stage 1 calculations, consider the dual implications for that portion of 

Table 4 which is reproduced in Table 5. Notice that when -h, = -200, as implied by the Stage 0 

dual relation, it is impossible to obtain a solution which satisfies the four direct constraints >© 
100 = yy, + 8 

(26) 200 = 81 

y,,8,-9 


TABLE 5 
Detail for Stage 1 (from Table 4) 
































-35 
ft ache “a 
ey 
A= 100 100 0, = 26 35 
y yy 8 
M = | 
M = 500 25] 20 6,=1¢ 
1 . ’ 
& Rly X,| -Yg=-1 Oo = -25 
5 
1 
Y -120 Oy = 35 
by 
-20 “hy /2 Yo 82 
0 25] by =1 
Y “39 Y Ro} Y Xo 
03 = 30 
by 
-20 “ho 





















yY -33 








36Notice that a solution could be secured by eliminating the nonnegativity requirement on Yj: 
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It is necessary, therefore, to eliminate the "bo implication" for hy . This can be done by 
altering -4 = -20 to -$, = -35. The first value for this variable is therefore struck out in 
favor of the second. But the remaining dual constraints then require, minimally, that 01 = 35 
and 54 = 7/5. No other dual variables need adjustment to satisfy the dual constraints, so that 
the remaining Stage 0 dual values are left undisturbed because of (a) the natural directional 
ordering and (b) the minimizing objective. 

The new dual solution now requires an alteration in the basis for the direct problem. 
The "Stage 0" check marks are eliminated from the cells for yy and Ry and new checks 
labeled " Y/" are entered in the cells for Xy and hy - A solution to this portion of the direct 
problem may now be arranged with 8, = 100, X= 20, -hy, = -120 and 39 = 0. 

The direct problem breakthrough at Stage 1 should, of course, be driven as far as pos- 
sible. Thus Yo= 120, Ro = 0, 33 = -4200 are values which can be arranged but trouble arises 
in the relation between Xo and “ho. The former nonnegative variable being pinned to a value 
of zero by 39 and nonnegativity on Ro cannot be used to accommodate the relation 


(27) -ho + Xp + Bo = 0, 





since -hy = -200 (as implied by the current dual solution) and 8 is not in the basis (and hence 
S> = 0).37 


Before adjusting to the next dual solution it is interesting to calculate the functional 
value, 


(28) 35 x 100 + 7/5 x 500 + 200 (0 + 5 + 0 + 25) = $10,200 





for the new solution to the dual which has just been arranged. This is lower by $1300 than the 
value in Eq. (25). As was shown in preceding sections, each of these dual solutions corre- 
sponds to the maximum for a direct problem which has fewer constraints than the one origi- 
nally stated. The alterations effected in the dual when passing from Stage 0 to Stage 1 tight- 
ened the corresponding zero stage direct problem to a Stage 1 set of conditions and, hence, 
lowered the optimum value for the functional. 

A Stage 2-to-1 solution can be arranged in the same way as Stage 1-to-0 —i.e., without 
disturbing previously satisfied direct program constraints. Since Og = 30, a new value -$9 = -30 
is selected to release ho from its bound. The only other dual variables which then require alter- 
ation are 55 and Y 1 which are each raised to 6/5. This removes the cell Ro from the basis 
and enters a " 2/" mark in the cell for ho. Adjoined to the previously secured direct program 
values (left undisturbed) are then the new values Yo= 120, Xo = 0, -hy = 0 and 3g = -4200. 

The functional for this new dual solution is lower by $1000 than its predecessor, since the 
alterations replace Uy = 5 by My = 0 without affecting any other variables in the functional. 

Table 6 is prepared for the Stage 3 calculations with the Stage 2 alterations in place. 

That is, the old values for the dual variables are replaced by their new ones, and the Stage 2 
(and, hence, Stages 1 and 0) basis entries and direct solution values are entered as indicated 
in Table 5 (and accompanying discussion). 


z 
E 
& 

F 
E 
og 
gs 
PS 








37Separate tables and accompanying notation herein are designed for exposition rather 
than efficiency in calculations. 
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TABLE 6 
Detail for Stage 3 



































70 
A= 100 100 |o, = 35 
& yVyY & 14/5 
64 = TAK 
M = 500 25] 20 -10/5 ~~ 
SB ly, %y\-%_=-¥  -oQ= -30 
Y,,~ 129 [35] 
120 
200. HY ~ Ye 
0 
Y -3, & %, 
- 4200 
Y -933 





















































An ambiguity now appears in that many values might be used to satisfy Y3; Ra; Xq; ~hg 
and -3 4° This ambiguity may be eliminated by noting that all operations are on basis sets so 
that unique solutions are required for each direct problem.?°® Jumping over to the b, noted in 
the cell for hy and tracing the implications of -hy = -200 back to the cell for Xg; it is seen 
that the bounding condition on hy must be eliminated. Even at Rg = 0, the price Cg = 25 makes 
it possible to purchase at most Xg = 168 units. 





38see (17) and (18) and related discussion. 
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To start the adjustment process, the Stage 0 value for %4 is replaced by - 4 = -50. 
This alteration makes it necessary to adjust as follows, O4= 50 and 5 4= 10/7, in order to 
secure satisfaction of the corresponding dual constraints. The ¥4 and Ry cells are, therefore, 
eliminated from the basis and the x 4 and h 4 cells entered in their place. The dual condition 
7,= 6; +] ROW requires an upward adjustment in ¥3- The old value -¥3= -1 is replaced by 
-Y3 = -10/7 and the cell for 34 thereby held in the basis. This new value for V3, then, re- 
quires upward revision in 93 and 5g. Before making these adjustments, however, it should be 
noted that the pair -3 = -25, O4= 50 creates a Ug = 25 and brings -hg to its lower bound of 
-200 with X3 (as before), being unable to absorb this amount because of the 3g limit of $4200. 
Hence, $3 = -25 is altered to - bg = -50. Necessarily, then, to satisfy the dual constraints 
03 = 50, 5g = 10/5 are required. Cells ¥3 and Rg are thereby removed from the basis and 
cell 83 enters in their stead.2? New values must then be secured as follows: ¥_* -10/5, 
-%o = -50, 49 Jy = 70 and 55 = 10/5. This, in turn, requires the new values -y , = -10/5, 

-$, = -70, o, = 70 and 6, = 14/5.") 

The solution to the entire direct problem is now readily arranged. The values 8, = 100, 
X, = 20, -hy = -120, Y= 120, 33 = -4200, Xo = 0, -hy = 0 are, of course, left undisturbed in 
this directed dual procedure. 3 - ho = 0 yields 83 = 0. 25 Xg - 3g = 0 gives Xq = 168 so that 
-hg = -168 and - Y =0. Therefore, &4= 168, X4= 0 and -hy = -168. Since, also, 35 = 0 it 
follows that Y5= 168 and Re = 0 so that 36 = -8400. 

The solutions in x and y clearly check with those shown in Table 2, and 01 = 70, 

6, = 14/544 clearly correspond to the cumulative values for the dual variables shown in the 
bottom row of Table 3. Noting that all w= 0 since all Fi41 = % the evaluation (cumulative) 
for the warehouse capacity is the same as that shown in the u-column of Table 3. 

The value 3g = $8400 may now be compared with the cumulative difference between 
cash receipts and disbursements shown in the bottom row of Table 2 (i.e., $12,600 - $4700 = 
$7900) and the final cash position of $9400. However, the discrepancy is easily reconciled. 
The value at the final sink must be adjusted by M = $500** to yield 3g - M= $7900 as the 
correct dual optimum. Alternatively the value for cash on hand at the conclusion of business 
must also include the amount of $1000 which was impounded in reserved cash to obtain 
$8400 + $1000 = $9400. Finally, the values shown for H;, 0,» ¢; and y, should be compared 
with those given in Table 3. It will then be seen that the special algorithm generated by the 
principles given in this paper effects its saving over the simplex method without losing this 
byproduct information which is often important in applied work. 





39 This creates a degeneracy which is here ignored. 

The change to -¢) = -50 is necessary to hold cell h2 in the basis. 

This gives the marginal compounded discount factor for internal yields discussed in [14]. 
It shows that by the end of Period 5 the amount to which an additional dollar invested will 
accumulate is $2.80. When the original dollar is "paid back," the net marginal incremental 
yield is 9/5 (as shown in Table 3). 

See footnote 31. The preceding discussion of dual and relaxed direct optima might also 
be altered by this amount, if all niceties are to be observed. 
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A large-scale experiment was designed to obtain the intuitive 
evaluations of intelligence reports from a large group of experienced 
Army officers. An interval scale of quality of the intelligence reports 
was constructed with a view to providing some mathematical models 
that describe in a simple fashion the way in which the components of an 
intelligence report combine to determine the value of the report, The 
results of the study indicate (a) that officers judge such intelligence 
reports consistently, (b) that they use criteria in judging intelligence 
estimates different from those they use in judging courses of action, 
and (c) that the judgments depend to some extent on the combat experi- 
ence of the officer, The different techniques of obtaining rating scales 
were found to produce essentially the same results. 











INTRODUCTION AND SUMMARY 

In order to design a good battle-area surveillance system, one would like to know the 
relative contributions of various kinds of intelligence obtained by the system to the over-all 
system performance. 

This study is designed to determine, by constructing an example, the possibility and 
feasibility of building a linear model (and a fortiori, to provide a procedure for doing this) to 
describe the relative importance of various kinds of information and mis-information that can 
be produced by an intelligence system. 

A large-scale experiment was designed to obtain the intuitive evz!uations of intelligence 
reports from a large group of experienced Army officers. Two such sets of quality judgments 
were obtained by different techniques and were then compared. These intelligence reports 
were also evaluated directly and compared with the values obtained when the stimuli judged 
were courses of action based on the intelligence reports. On the basis of these judgments, an 
interval scale of quality of the intelligence reports was designed with a view to constructing 
some mathematical models which describe simply the way in which the components of an intel- 
ligence report combine to determine the value of the report. 

The results of the study indicate that officers judge such intelligence reports consist- 
ently, that they use criteria in judging intelligence estimates different from those they use in 
judging courses of action, and that the judgments depend to some extent on the combat experi- 
ence of the officer. The different techniques of obtaining rating scales were found to produce 
essentially the same results. Using the intelligence estimate data taken from combat officers 
only, a linear model was constructed to evaluate intelligence estimates in terms of the 





“Manuscript received 22 May 1959. 
TThis experiment was conducted under support from Project MICHIGAN for the U.S. Army 
Signal Corps under Department of the Army Prime Contract Number DA-36-039-SC-52654, 
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contributions of its components. This model was found to be a good predictor when validated 
against additional data. This result indicates that the linear model technique may be useful in 
the general problem of evaluating intelligence systems. 

The particular examples of this study are basically concerned with a World War II type 
(preatomic) divisional front. This study is intended to be methodological in nature; it is rea- 
soned that the techniques are most appropriately developed in a situation which has been 
encountered and dealt with in the past. This is particularly pertinent since the basic data used 
in this study are judgments of Army officers who have had extensive experience in World War 
II and Korea, i.e., nonatomic warfare. It is assumed, however, that the methods developed 
here will be extendable to atomic ground warfare and also to target selection problems in an 
atomic ground-to-ground or air-to-ground bombardment. The procedures by which such 
extensions can be made will be considered in the last part of this article. 


BACKGROUND AND HISTORY 

In designing an intelligence system, the designer must make some decisions regarding 
the relative priorities of the different kinds of information his system is going to produce. For 
example, the information may concern infantry, artillery, or armor, etc., and may describe 
numbers of units, location of units, identification, etc. This priority scale must be conditioned 
by the degree to which each of these kinds of information determines the adequacy of the sys- 
tem for a particular purpose. A sensing device of such a system, which is intended to respond 
to the presence of tanks, could conceivably be constructed to be sensitive enough to register in 
the presences of a large-sized stone. A device of such sensitivity would seldom fail to register 
if an enemy tank appeared in its range. However, it would register many times when there was 
no tank present. This kind of erroneous report is called a "false alarm.'' On the other hand, 
if the sensitivity of the instrument is too low, the device might never report a tank unless a 
tank is actually present, but many times when a tank is present it will fail to report it. This 
kind of error will be called a ''miss."' Any enemy unit thatis reported but mislocated willbe said 
to be "'displaced,"' and a correct identification and location will be called a "hit.'' In order to 
construct the intelligence-gathering system in a proper fashion, it is necessary to be able to 
evaluate the output of ..1e system, and, knowing the value of the output of the system, to be able 
to evaluate the components in terms of their contribution to the value of the system. The 
rationale for such a view has been described in a previous report [1]. 

In brief, it is believed that the output of a system can be described as being composed 
of various kinds of information about the enemy (in the case of the present system, these kinds 
of information are hits, false alarms, displacements, and misses, for front line and reserve re 
infantry and armor) and furthermore, that the contributions of each of these kinds of informa- ; 
tion to the quality of the whole system can be added together to describe the quality of the sys- 
tem. That is, a linear combination of the components is expected to be a good approximation 
to the (unknown) way in which the components actually combine. In the above mentioned report 
[1] it was shown that such a linear model was feasible in an extremely simplified situation, 
where only a region of regimental size was considered and only infantry and artillery units 
were displayed. 

It was decided to attempt to construct a linear model in a situation which was suffi- 
ciently realistic to be considered akin to an actual battle situation by experienced Army offi- 
cers. There was some question as to the kind of criterion that was used by the officers in 
their role as judges in the previous experiment. They were instructed to judge an intelligence 
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estimate with regard to its quality "considering the good of the Army as a whole." It is well 
known that different kinds of intelligence have different relative importance at the various 
command echelons. For example, it is very important (to pick an extreme case) for an infantry 
soldier to know the location of a particular aggressor soldier so that he can fire at him and 
avoid that man's fire; the location of that soldier, however, makes very little difference to a 
division commander except insofar as he knows that his troops have that kind of information 
and that therefore they become a more efficient fighting force. In order to overcome the dif- 
ference in importance of intelligence to the different echelons, the officers were asked to eval- 
uate the reports from the point of view of "'the good of the Army as a whole." It was hoped that 
a high-ranking Army officer, having previously operated at lower echelons, would be able to 
judge the value of an intelligence estimate to all echelons, and that this would then free the 
quality judgment from dependence on the echelon context in which it was obtained. Since the 
majority of the officers in the previous study had had little or no previous command experience 
at a high echelon,! it was doubtful that they were able to assume this grand, over-all view of 
"the good of the Army as a whole."' Furthermore, there was some doubt as to whether officers 
could make use of their military experience in judging situations which were as over-simplified 
as the ones in that experiment. Consequently, it was decided to conduct this experiment using 
officers with sufficient combat command experience at a sufficiently high rank to increase the 
probability that they would be able to make judgments from the point of view of the good of the 
whole Army in the context of the situations that were used. 

Probably the most difficult problem which had to be considered prior to conducting the 
present experiment was the question of obtaining an appropriate criterion for measuring the 
value of an intelligence estimate. It was concluded that ultimately it would be necessary to 
depend on the opinions of experienced Army officers on the assumption that they were experts. 
This was taken to be the most reasonable solution to a problem which did not seem to have any 
totally adequate solutions. Presumably, the quality of an intelligence estimate depends on the 
extent to which such an estimate will increase the effectiveness of an Army unit. It was argued 
that one way of measuring the effectiveness of the estimate in helping a commander win a battle 
is to evaluate the quality of a course of action which is drawn up on the basis of such an intelli- 
gence estimate. For this purpose, officers might be shown the true agg: essor disposition 
together with the course of action drawn up on the basis of various intelligence estimates of 
the aggressor disposition. Their evaluations of each of these courses of action would then also 
serve as a measure of the quality of the corresponding intelligence estimate. In response to 
this argument, part of the effort of this experiment was devoted to constructing an optimal 
course of action for a particular intelligence estimate and then getting officers to evaluate the 
courses of action. This would constitute another measure of the quality of the intelligence 
estimate. The officers could also be asked to assign a value to the intelligence estimate 
directly. This provides for each intelligence estimate a judged value of its quality and a judged 
value of the course of action designed on the basis of that estimate. These two measures could 
then be compared and thereby provide an estimate of their similarity. The values assigned to 
the intelligence estimates are taken to be the primary measure of the quality of the system that 
produced the intelligence. This orientation results from the assumption that an intelligence 
estimate ordinarily contains information which, although not used in constructing a course of 





1 These officers were for the most part students at the Infantry School at Fort Benning or at the 
Artillery School at Fort Bliss. 
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action, is nevertheless of great value, particularly to lower echelons, since it increases their 
efficiency as fighting units. 

When values are assigned to courses of action, the interpretation of the authors is that 
only a relatively small amount of the available intelligence is being judged so that intelligence 
estimates of different qualities (as defined above) may give rise to identical courses of action. 
Consequently, the primary measure of intelligence in this experiment is the judgment of the 
intelligence estimates directly, and the focus of the analysis will be on this measure. 

It is to be noted that a more realistic measure of quality of intelligence would involve 
the adequacy of intelligence estimates in actual use [3,5]. Since this avenue of experimentation 
did not appear a satisfactory one, we chose as a second best the judgments of Army officers 
regarding their expectations of how good such intelligence would be in actual battle. 

Since the data indicates good agreement among the judges, we can derive some inferen- 
tial comfort in the assumption of a similarity of the measure we used as compared to a meas- 
ure derived from a real war. 


OUTLINE OF EXPERIMENT 
The description of the experiment is divided into two major categories: (a) the con- 


struction of the experimental materials and (b) the experimental procedure for the actual col- 
lection of data. 


Preparation of Test Materials 

Enemy Situations - A small group of officers on the faculty of the Command and General 
Staff College (C & GSC) at Fort Leavenworth was assigned by the college to aid project person- 
nel in the construction of test materials. These officers were asked to select three different 
terrain instances which were, however, of the same general West European type.“ The three 
different terrain instances were requested so that the results of the experiment would have 
generality at least with respect to a particular type of terrain. On the basis of a series of 
conferences with project personnel, the C & GSC officers constructed seven different enemy 
situations on each terrain which satisfied the following conditions: 

1. No two enemy situations were to give rise to the same optimal course of action on 
the part of the friendly forces. This was required so as to permit a comparison of the intelli- 
gence estimates by the method of comparing the courses of action that were derived from them. 

2. The situations were to be militarily reasonable. 

3. Friendly forces in approximately division strength are in each case attacking to gain 
some fixed objective. 

4. For the enemy situation, aggressor order of battle was to be used. 

5. The only enemy units to be displayed were infantry and armored units. Artillery 
was not to be displayed but was to be assumed present in strength consonant with the strength 
of the armor and infantry units displayed. (It was the consensus among these officers that 
enemy artillery had little or no effect on the choice of the optimal course of action.) 


6. Fixed obstacles (mine fields, entrenchment, etc.) were to be kept constant on a given 
terrain instance. 





2Map numbers of the three terraininstances are: CGSC 50-119 Hof Lobenstein; CGSC 50-124 
Fulda Hershfield; CGSC 50-134 Pont-a-Mousson-Nancy. 
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7. On each terrain instance, one additional intelligence estimate of an enemy situation 
was taken to be just the terrain map with all the fixed enemy obstacles. This will hereafter be 
referred to as the "blank terrain." 

An example of one of the seven enemy situations used on Terrain Instance I is shown in 
Figure 1. 


Figure 1 - One of the seven enemy situations used on Terrain Instance I 
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Courses of Action - Corresponding to each enemy situation (including the "blank ter- 
rain") constructed as above, the officers then constructed an optimal course of action for the 
friendly forces. To assure that these courses of action were optimal, they were then reviewed 
by a panel of officers from the faculty of C & GSC (the "Murder Board"). A requirement on 
the courses of action was that for a particular terrain instance, the length of the (friendly) 
division front was to remain constant. (Normally, the size of a sector assigned to a division 
commander varies with the difficulty of his task, i.e., the terrain and the strength of opposi- 
tion.) The friendly forces would be reinforced or have units detached in accordance with the 
principle of economy of force. The courses of action stimuli were presented as avenues of 
attack and disposition of friendly forces overprinted on a military map, with an extract from 
Corps operation order (including organization for combat) attached. This extract served to 
describe the essential details of the planned operation other than those that could be indicated 
by sketches on the maps. Each of the seven courses of action which had been constructed as 
being optimal against a corresponding aggressor disposition was presented as a possible 
course of action against each aggressor disposition. 

The course of action designed for use against the enemy situation, appearing as an 
intelligence estimate in Figure 2, is given in Figure 3 as a possibility for use against the 
aggressor disposition of Figure 1. . 


Se oe 





—— 


Figure 2 - An intelligence estimate of the situation of Figure 1, which 
was provided to the officers in the form of a transparent overlay 
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Figure 3 - The course of actiondesigned for use against the enemy 
situation of Figure 1, which was constructed on the basis of the 
intelligence information about that situation indicated in Figure 2 
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Intelligence Estimates - The intelligence estimate stimuli were constructed by printing 
one aggressor disposition on a terrain map, printing another on plastic, and overlaying the 
plastic on the map. This procedure resulted in 42 stimuli corresponding to each terrain 
instance. Each of the seven aggressor dispositions was paired with each of the remaining 
ones, except that no such disposition was ever paired with itself. In addition, for each disposi- 
tion, a plastic was prepared showing only the permanent features such as barbed wire, mine 
fields, etc. This was also used as an intelligence estimate and will be referred to as the "blank 
terrain." 


Figure 2 shows an overlay which was used as an intelligence estimate when overlayed 
on the enemy situation illustrated in Figure 1. The intelligence estimate is that situation 
against which the course of action indicated in Figure 3 was designed as optimal. 

A questionnaire was also constructed which presented verbally the outputs of intelli- 
gence systems in much the same sense as the map materials presented the outputs visually. 
This was done in the hope that it might—if the results agree with those of the map exercise— 
serve as a simplified procedure for evaluating intelligence or surveillance systems. 


Experimental Subjects 

The officers used as subjects in this experiment were colonels and lieutenant colonels, 
students at the Army War College. All had previously attended the Command and General Staff 
College. At the time of the experiment 181 people were available, of whom 139 were from the 
combat arms (armor, artillery, infantry, engineers, and marines), and 42 were from other 
sources, including civilians, Air Force, and Navy. 


Procedure 


The Army War College made its student body available for the present experiment for 
a period of 2 days. Approximately 1 week preceding the experiment proper, each officer 
received by mail each of the three terrain maps, together with scenarios, and introduction to 
the experiment, and instructions to study the material. The morning of the first day included 
an initial orientation to the procedure for the remainder of the experiment. Each officer was 
given a schedule describing his tasks for the remainder of the experiment and stating where 
and when he was to perform them. The experiment proper followed immediately after that. 

The morning of the first day, half of the officers were asked to evaluate three sets of 
intelligence estimates (one set per terrain), and the other half evaluated courses of action. 
During the afternoon those officers who had evaluated courses of action in the morning evalu- 
ated intelligence estimates, and similarly those who had evaluated intelligence estimates in the 
morning evaluated courses of action in the afternoon. 

The task of evaluating a set of intelligence estimates involved considering a set of 
maps, each with the same aggressor disposition overprinted on it, but each having a plastic 
displaying a different one of the remaining six dispositions on that terrain overlaid on it. In 
addition an overlay showing the "blank terrain" was also included. The instructions to the 
officer were to assign to each of the experimental stimuli a value rating between 0 and 100 
which would describe the quality of the intelligence estimaie that the overlay represented, con- 
sidering the situation printed on the map as "'true."' When the officer had finished rating these, 
he proceeded to another place where he was given a similar task for a different terrain 
instance. When the officer completed the second set he then proceeded to the third set. When 








EVALUATING INTELLIGENCE SYSTEMS OUTPUT 33 


all three sets were completed, there was an hour lunch period, and the officers returned after 
lunch to finish the day's experiment. 

Evaluating the courses of action involved considering a set of maps, each with the same 
aggressor disposition overprinted on it, but each having a different course of action overprinted. 
These courses of action were the ones that were designed as optimal for the various aggressor 
dispositions. The instructions to the officers were comparable to those they were given for 
judging intelligence estimates. 

For each terrain each officer received one true situation, and each had to evaluate 
seven intelligence estimates of that situation and also seven courses of action fought against 
the same situation. 

The morning of the second day was devoted to comparing pairs of intelligence estimates 
and judging in each case which of the pair was the better intelligence estimate of its true situ- 
ation. All of these intelligence estimate stimuli (not including the ones having a blank terrain 
plastic) had been roughly rank ordered according to quality of the intelligence by a group of 
officers from Fort Wayne, as well as by some officers stationed locally. From these, 60 stim- 
uli were selected in such a way as to cover approximately the complete range of quality as 
determined by the preliminary rank order and to be fairly evenly spaced in order of quality. 
The stimuli in this ordered set of 60 were then paired so that each stimulus appeared in a pair 
with each other stimulus which was eight or fewer steps removed from it in the order of the 60 
stimuli. There were 424 stimulus pairs constructed in this manner. In addition to this set of 
stimuli, some of the intelligence estimates at the extreme ends of the rank order were used in 
the same pairs more than once, and one pair of stimuli was used a great number of times to 
evaluate sequence and position effect of placement of the stimuli. In all, 500 pairs of stimuli 
were prepared. The stimuli were placed on tables in sets, with 10 pairs per set. Each officer 
was then assigned to judge the pairs belonging to 10 of these sets, that is, he judged 100 pairs 
in all. Since 500 stimulus pairs had been constructed, there were 50 of these sets altogether, 
but the scheduling was so arranged that, with very few exceptions, no two officers judged more 
than one set in common. 

In all cases the assignment of officers to tasks was random within the group of combat 
officers and also random within the group that was considered to be noncombat officers. Fur- 
thermore, in the assigning to the officers of the true situations for them to consider in the task 
of Day 1, the assignments were made so as to ensure as far as possible that different officers 
had not more than one true situation in common, and that the sequence of tasks was different 
for different officers. 

The afternoon of the second day was devoted to a debriefing session and to the conduct 
of a supplementary experiment which involved obtaining the judgments of these officers as to 
the quality of intelligence estimates that were only verbally described. This involved present- 
ing the officers with questionnaires that contained statements as to the number of units in each 
of several categories that were hit, missed, overestimated, or displaced by a particular intel- 
ligence estimate. The results of this questionnaire study are not included in the present paper. 

After the morning session of the first day and again after the afternoon session of the 
first day, the officers were requested to fill out questionnaires which were intended to provide 
some insights into the considerations that were involved in the making of such judgments by 
Army officers. In general, these questionnaires confirmed the preconceptions of the experi- 
menters but provided few additional insights, and no analytic use has been made of them. 








34 CALDWELL, COOPER, SCHOEFFLER, AND THRALL 


RESULTS AND DISCUSSION 
Throughout this experiment, the records of the officers from the combat arms and those / 

from the noncombat arms were kept separately. It was hoped that it would prove possible to 

combine the data from these two groups and thereby make use of a larger sample. However, 

various statistics implied that the judgments made by the two groups of officers were not from 

the same populations. Consequently, the results of the experiment, as reported here, deal pri- 

marily with the judgments of the combat officers only. . 
In judging the course of action stimuli, each officer was presented with a set of seven 

different friendly courses of action, each fought against the same aggressor disposition. A 

course of action consisted of avenues of attack overprinted on the map together with an abstract 

of the Corps order which was stapled to the map (Figure 3). First, the officer rank ordered 

these in terms of their respective quality, and then he assigned a value between 0 and 100 to 
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each. Subsequent analysis° indicated that: (a) when ratings of all the courses of action fought 4 
against a particular aggressor disposition were averaged, the value obtained was approximately 
the same for all of the aggressor dispositions, and (b) some courses of action tended to be E 


rated as good regardless of the disposition of the forces.against which they were fought, and 
others tended to be rated as poor. In general, the courses of action that were rated as good 
were the ones that included greater friendly strength. That is, a particular course of action 
was judged to be good or poor seemingly independently of the aggressor disposition against 
which it was to be used. 

These results imply that the officers tended to rate the courses of action in a set in 
comparison to one another, rather than rating them on some sort of absolute scale, and that 
the numerical values assigned to the courses of action were therefore not to be interpreted as 
absolute measures of their qualities. The results also imply that in this case the officers were 
not giving high priority to economy of force as they had been instructed to do in the introduc- 
tory briefing. An alternative explanation is that perhaps the officers interpreted ''economy of 
force" to mean adequate development of troops which were assigned by Corps. In this case, 
the implication is that the officers assumed that their concern was only with their own division 
front, and that shortages of troops at other points were beyond their cognizance, and a further 
implication is that the officers considered the additional forces more important to the success 
of the mission than other features of the course of action. bs 

The interpretations of the results regarding the evaluations of courses of action must, 
however, be considered in the light of one major procedural defect. For all of these judgments, 
the officers were forced to rank order and evaluate seven courses of action in the very short 
period of approximately 1 hour. This was reported by some of the officers to be an impossibly 
short time for adequate performance. 

Essentially the same effects were observed when the officers were asked to rate the 
respective values of seven different intelligence estimates of the same true aggressor disposi- 
tion. The interpretation of Effect 1 is the same as above. The reason some intelligence reports 
received a high average rating and others a low average rating when they were used as esti- 2 
mates of each of the aggressor dispositions is not obvious, in view of the fact that each of the 
aggressor dispositions was used as an intelligence estimate of each of the other dispositions. 








3The statistical basis for these conclusions are presented in Report No, 2144-272-T, The Uni- 
versity of Michigan, Willow Run Laboratories, April 1959. The reader is also referred to 
this report for additional detail regarding other aspects of this experiment. 
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As has been described above, each intelligence estimate also served as a "true" situa- 
tion, and for each of these, the Command and General Staff College constructed an optimal 
course of action. Part of the experiment was designed to permit a comparison of how well the 
officers' rating of an intelligence estimate agreed with the officers' rating of the quality of the 
course of action that was constructed on the basis of this intelligence estimate. The correla- 
tion was found to be essentially 0 (r = -0.094). Since it had already been discovered that a 
consideration of the amount of friendly strength available was of overwhelming importance in 
evaluating a course of action, another comparison was made which only considered courses of 
action that had the same strength. This new measure again showed no important relationship 
between the rating given to a course of action and the rating assigned to the intelligence esti- 
mate on which it was based. 

Apparently, the fact that an officer judges an intelligence estimate to be a good one does 
not imply that he will be able to construct a good course of action on the basis of that estimate. 
Furthermore, it appears that, in general, the quality of the course of action is determined pri- 
marily by such factors as the amount of strength available and the nature of the terrain, 
although the particular situation does determine the quality to some extent. A major value of 
an intelligence estimate must be something other than the basis of a course of action. One 
possibility is that when an officer says that an intelligence report is good, he is referring to 
the extent to which the report will help his subordinate units fight more effectively rather than 
to the extent to which it will assist him in deciding on a particular course of action. 

Each officer was also presented with 100 sets of two stimuli, each consisting of an 
aggressor disposition and an overlaid intelligence estimate of it. He was asked to choose the 
better of the two intelligence estimates in each case. These choices will be referred to as the 
pair comparison judgments. On the basis of preliminary testing, a tentative rank order was 
made of the stimuli in terms of their quality as intelligence estimates. Of these, 60 were picked 
to represent the whole range of quality, and these were each paired with each of the 16 adjacent 
stimuli in the tentative rank order scale for pair comparisons. The proportion of times a 
stimulus was preferred to one it was compared with, was taken as the measure of the differ- 
ence in quality of the two intelligence estimates. These proportions were transformed to nor- 
mal deviates, and the resulting matrix was solved to yield a scale of quality by the procedure 
outlined by Gulliksen [4]. Evidence for the adequacy of this scaling procedure was available 
from a high correlation (r = 0.717) that was obtained when these scale values were compared 
to the values obtained from the direct rating procedure described above.* 

The scale values arrived at in this way were tentatively assigned a 0 point and a unit of 
measurement by fixing the values of two stimuli of the set on the basis of the value ratings they 
had received in the part of the experiment described above. This also served to place the val- 
ues of the stimuli in the range of 0 to 100. Each score was then divided by 100 to compress the 
range to 0-1, and a model was constructed to describe the values of the components of the 
intelligence estimates according to the philosophy described in an earlier report [1]. 


The Model 
Since we had some a priori knowledge of the nature of the contribution of certain com- 
ponents of the system to the quality of the intelligence report, viz., that correct identifications 





4The reliability coefficient ofanofficer's judgments of such material is estimated to be in the 
neighborhood of 0.8 from data obtained in a previous study [1]. 
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and displacements [1] made a positive contribution to the value of the intelligence estimates 
and that false alarms made negative contributions, we had available a minimum criterion for 
the acceptability of a model, viz., whether it gave rise to the proper algebraic signs on the 
values of the components of the intelligence report. A series of different models of varying 
degrees of complexity was attempted, and the following one was selected on the basis of ade- 
quacy of fit and degree of simplicity. 


p = -0o1i * 81°11 * 22°11 | *3°2i + A421 * 85°21, B6°si * *n%i * *8%si 
y thi i 33 





in which 


the value of the i-th intelligence estimate computed from the model, in distinc- 
tion to b; which is the scale value of the i-th estimate; 


the contribution to the quality of the estimate by the j-th component; 


the number of front line infantry battalions correctly located on the i-th intelli- 
gence estimate; 


the number of front line infantry battalions reported on the i-th intelligence esti- 
mate in excess of the number that actually appeared; 


the number of front line infantry battalions reported but incorrectly located on 
the i-th intelligence estimate; 


total number of front line infantry battalions that appear in the i-th true situation; 


the number of reserve infantry battalions correctly located on the i-th intelli- 
gence estimate; 


the number of reserve infantry battalions reported on the i-th intelligence esti- 
mate in excess of the number that actually appeared; 


the number of reserve infantry battalions reported but incorrectly located on the 
i-th intelligence estimate; 


total number of reserve infantry battalions that appear in the i-th true situation; 
the number of armor battalions correctly located on the i-th intelligence estimate; 


the number of armor battalions reported on the i-th intelligence estimate in 
excess of the number that actually appeared; 


the number of armor battalions reported but incorrectly located on the i-th 
intelligence estimate; and 


. = total number of armor battalions that appear in the i-th true situation. 


(See footnote>.) 





>The values of the t;, Cj,e;, amd s; were determined by a group of officers from C & GSC in 
conference with members of the project staff. They were derived from looking at the ''true"' 
situation and its plastic overlay, which presented the intelligence information in the following 
way. First the total number of battalions of each kind were counted on the ''true'’ map. When 


(Continued) 
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As a result of errors made in the process of printing the maps, of the 60 stimuli which 
were used in the pair comparison matrix, only 52 were used for constructing a scale. This 
scale was constructed according to a method outlined by Gulliksen [4] for solving an iricomplete 
pair comparison matrix. On the basis of ratings assigned to two of these stimuli the scale was 
transformed to a range of 0 to 1 (the stimuli were chosen so that there would be one that 
received a high rating and one that received a low rating and the two were in the same set, i.e., 
were based on the same true situation). On the basis of these transformed scale values (called 
b;, i=1, 2,... 52) a solution was attempted for the a, so as to minimize the quantity 


26 
. + a 
(bo; - bai)", 
ici 


i.e., using the 26 even numbered stimuli from the set of 52. This solution yielded the data, as 
follows. 


ag 0.555 
a, = -0.307 
ay 0.090 
ag = -0.029 
ag = 0.014 


a5 = 
a6 = 
an = 
ag = 


On the basis of this solution, ''predicted" qualities were computed for the odd-numbered 
26 intelligence estimates. A Pearson product-moment correlation of these predictions with 
the observed judged quality of the intelligence estimates yielded r = 0.671. A rank order cor- 
relation yielded p = 0.700. The indication is that the prediction capabilities of this form of the 
equation are relatively good and exceed chance. 

Since the problem of predicting was satisfactorily settled, the values of the a; were 
recomputed using the complete set of 52 stimuli. This is then the best estimate available for 
the set of parameters. The values of the a; and the meaning of the a, are shown in Table 1. 

The signs of the parameters are appropriate for front-line infantry and for armor, i.e., 
hits and displacements have positive value, and false alarms are damaging to the quality of an 
intelligence report. In the case of reserve infantry the signs are reversed. However, for this 
category, all of the parameter values are near 0, and a test of the hypothesis that the set of 
parameters obtained is not significantly different from one which has the same a values, 
except that ag, a4, and a, are equal to 0, yields a nonsignificant value of p.6 Another check on 





(Continued) 


the classification of an infantry unit as ''front line" or ''reserve'' was doubtful, a decision was 
made on the basis of the echelon which was most likely to be in control of the battalion. If it 
seemed that the battalion was attached to a particular regiment it was designated ''front line," 
and if it seemed that the battalion was attached to the division directly, it was designated 
"reserve.'' Further, whenever a unit was mislocated, the officers were asked to judge whether 
such a mislocation still had some positive value, or whether the quality of the intelligence 
report would be enhanced if the unit were not reported at all. If the first alternative was 
judged to be true, the unit was counted as a displacement, If the second alternative seemed 
reasonable the unit was counted as a''miss" in its proper location, and a '''false alarm" in its 
repegted location, 

The hypothe sis that all of the parameters obtained are not significantly different from 0 was 
also tested, and yielded a significant value of F. (F = 6.45 with 9 and 41 degrees of freedom.) 





CALDWELL, COOPER, SCHOEFFLER, AND THRALL 


TABLE 1 
Values and Meaning of the a; 





Amount of Contribution 
per Battalion to 
Quality of the 
Intelligence Estimate 





0.438 Front line infantry hits 
-0.052 Front line infantry false alarms 
0.258 Front line infantry displacement 
-0.038 Reserve infantry hits 


0.047 Reserve infantry false alarms 


-0.119 Reserve infantry displacement 


0.367 Armor hits 


-0.070 Armor false alarms 





0.255 Armor displacement 














the reasonability of the solution is the following consideration. In the case where only hits 
occur and all units are observed, i.e., intelligence estimate is perfect, the computed value of 
the intelligence estimate should equal 1. This requires that the sum of Ag, ag, and ag equal 1. 
Since ag is essentially 0, and ap + ag = 0.805, a reasonably good approximation to the require- 
ments of the equation is obtained. It must be remembered that the original scaling of the val- 
ues into the 0-1 range was somewhat arbitrary and that the solutions are dependent on the 
particular scaling used. A constant multiplication of the scale values produces a multiplication 
of the solution vector by the same constant. If the constant is equal to 1.24, the requirement is 
fulfilled that a perfect intelligence estimate receives a value of 1.’ If we still consider the 
values of ag, and a 4“ and ae equal to 0, we get a new solution vector, as follows. 


ag = 0.543 as 0 
a, = -0.064 
ao = 0.320 


ag = 
an = 
ag = 0 ag = 


a4 0 








?This procedure has the drawback that the new solution is no longer a least squares fit to 
the data, but has the advantage that the actual value of a parameter has greater intuitive 
significance. 
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As long as the scale values are multiplied by the same constant as the solution vector, 
the degree of prediction is not changed, and the correlations between predicted and observed 

values of the stimuli are not changed. 

The actual values of the ay imply the following: 

1. In judging the quality of intelligence systems, officers tend to be unconcerned with 
the information that is obtained regarding the reserve infantry units. 

2. The fact that a system may occasionally overestimate the number of enemy frontline 
infantry battalions and armor battalions does not appreciably degrade the judged quality of the 
system. 

3. It seems most important that the system correctly report frontline infantry units, 
and, almost as important, that the enemy armor be correctly reported. 

4. Information regarding frontline infantry is seriously impaired in judged value by 
having location incorrectly reported, but such information is nevertheless helpful. Information 
regarding armor battalions is also degraded in judged value by having the units mislocated, but 
this effect is not as pronounced. 

For all of these four statements, it must be remembered that the relative measures on 
the infantry and armor are in terms of the particular units that both kinds of forces are meas- 
ured in (in this case battalions), and also that a larger amount of mislocation was permitted 
for armor than was permitted for infantry before the information was labeled a displacement. 
(See footnote ”). 

On the basis of the model just described, it has proved possible to reach some interest- 
ing conclusions regarding the use of intelligence on a particular kind of terrain in an offensive 
battle of preatomic warfare. On the basis of the obtained data, it seems reasonable to assume 
that a similar solution could be obtained for a different kind of terrain, different state of battle, 
and a different size unit. In general, a model of this type should be constructable for any sur- 
veillance system in any kind of warfare so long as there is general agreement as to what con- 
stitutes an expert. In this experiment, the agreement was to consider as an expert any student 
at the Army War College from one of the combat arms. The dependence of the procedure on 
the availability of a corps of experts is a serious drawback when an attempt is made to con- 
struct such a model for a type of warfare where no one feels very confident of his abilities for 
judging the adequacy of an intelligence display. However, it seems that this dilemma is intro- 
duced by the dependence of the estimated quality on an ill-defined objective. Thus, for pre- 
atomic warfare, the intended use of the output of a surveillance system was not stated as a 
1-variable function. The intelligence was intended to direct artillery fire, to avoid meeting 
large concentrations of troops in an advance, to permit preparation against counterattack, to 
avoid roadblocks, etc. To the extent that a single purpose can be described for the output of a 
surveillance system, it becomes possible to develop other techniques for evaluating the quality 

of an output. 

Consider the requirements on a surveillance system whose sole function is to yield a 
plausible basis for weapon assignment in atomic warfare, and that an SOP for such assignment 
is available. Under such conditions a clear criterion of quality is available, viz., the system 
that produces all of the specified input for the SOP is perfect, and the degradation of perform- 
ance is related to the amount of information not provided or provided incorrectly. If the func- 
tion of the system is to yield information which, when used in accordance with the SOP, will 
yield the greatest probability of kill with the assignment predetermined, then the quality is 
measured by the amount of destruction derived from the system. 







CALDWELL, COOPER, SCHOEFFLER, AND THRALL 


Granted that even in these cases problems arise regarding the measuring of such things 
as, for instance, ''amount of destruction," at least the problem of the effect of experience in 
making such judgments is deferred one step. 
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FIXED-COST TRANSPORTATION PROBLEMS* 
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This paper formulates a fixed-cost transportation problemas aninteger 
program, describes some of its special properties, and suggests an 
approximate method of solution. Examples are given to demonstrate the 
approximation technique. 








INTRODUCTION 
A transportation problem can be described analytically as follows: 


M: Minimize >. >». a..x.. 
nimize =; j Aaj *ij 


subject to the conditions 


j= 


for i=1,...,m and j=1,...,n, where it is assumed that a5 e; = i dj; this being a 


necessary and sufficient condition for the existence of a solution. 

In one interpretation of this problem there are m "reporting activities" which are in 
excesses e; of a certain item, and n which have deficiencies d. of that item. The a;; repre- 
sent unit costs of transporting from the jth excess point to the jo deficiency point. The prob- 
lem is, then, one of distributing the items, in order to have no excesses and no deficiencies, at 
a minimum total cost of transportation. 

Many efficient computing methods for solving this special form of linear program are 
available [1, 2,6, 7], and we note that solutions can always be given in Xj j which are nonnegative 
integers [5]. 


DESCRIPTION 

In the linear transportation problem, M, the cost associated with each route (i,j) is 
aij xij , 1.e., proportional to the number of units Xij which are transported along that route. 
Instead, we consider a transportation problem with fixed cost. Namely, the cost of sending no 
units along route (i,j) is zero; but any positive shipment incurs a fixed cost plus a cost propor- 
tional to the number of units transported. Analytically, this can be described in any one of a 
number of ways. Expressed geometrically, we take the curve of Figure 1(b) as the total 


*This work was supported by the Bureau of Supplies and Accounts of the Department of the 
United States Navy, under Contract No. Nonr-2928(00). 
Mathematica, 76 Nassau Street, Princeton, New Jersey, a subsidiary of Market Research 
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transportation cost along route (i,j) (as versus the cost associated with route (i,j) in the 
linear transportation problem in Figure 1(a)). Thus, the cost of sending no units is zero, but 
of sending one or more is bij plus a sum proportional to the number of units sent. 


COST COST 


SLOPE aj, 


SLOPE aj; 














(a) xij (b) xij 


Figure 1 


Our analytic formulation takes the form: Given positive integers e; and dj and con- 
stants aij and bij» find Xij and Vij to 


(1) P: Minimize ) (&,y) = 2i oF a; xj + bi Vij) 
subject to the conditions 

(2) mF Xj = ej, 55 xij = di, 

(3) 0= Vij =i, Vij integer , 

and 

(4) 0= xj = Mj Vij 

for i=1,...,m and j= i, .,n; where 

(5) m;; = min e;; d;). 


It is again assumed that ©. e. = ©. d.. 


Note that conditions (3) imply that the yj are 0 or 1 and, due to (4), if ¥ij= 0 then 
X45 = 0. If ¥ij = 1 condition (4) is innocuous since x;; cannot, in any case, exceed min e;, dj). 
Thus Vij can be thought of as a "'valve'' which controls route (i,j): only when it is positive can 
the route be used. 

Many other analytic formulations of this problem are possible. P seems to be the most 
straightforward; however, the complete statement of condition (3) is not absolutely necessary, 
since for any optimal solution Vij will certainly not exceed 1. We have retained the condition 
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O=y; ij = 1 despite this fact since it adds certain structure to the problem P which will be 
seen to be of interest (see Theorem 4 which follows). 

The program P isa linear program "in integers," for the y,; are constrained to be 
integers. Gomory [4] has devised an algorithm for finding optimal solutions in integers to 
general linear programs. This algorithm proceeds initially to a ‘first rational solution"—a 
solution to the problem in which the integer constraints are ignored. If this solution does not 
have all variables in integers, then new inequality constraints are automatically generated and 
added to the original program. Geometrically, these new constraints "'cut"' into the convex 
body determined by the original program in such a way that no integer points are cut off. The 
procedure is continued until a best integer solution is found. This algorithm, however, is 
based on the use of the simplex method for solving general linear programs and is thus, at 
best, an expensive process when applied to a transportation type problem. Conceivably, 
Theorem 4 may lead to a specialization of the Gomory technique which could be used for find- 
ing an exact solution to P efficiently. However, we shall propose, in what follows, an approxi- 
mation based on the very efficient algorithms available for linear transportation problems, 
which also provides upper and lower bounds on the exact solution. 


THEORY 
THEOREM 1: There exist solutions {x}; 2 vit to P such that the xe, are integers 


Proof: Let (x? , vijh be a solution to P and define the problem (in variables Xj) 
Q: Minimize &, >. c.. x,. 


subject to the conditions 


“ oe es ss 
25 Xj =e, 2 Xj dj, and Xj 0 
where 
= : Om <x 
Cj = aij if Vij = 1 
es : ay: 
=M if Vij = 0 


and M is very large. Then {x; °} must solve Q, and any set {x,)> Yj ih? where x, . solves Q, 
must in turn solve P. But Q is a linear tranepertaiten problem and io integer sotetions [5]. l 
Thus there exists a solution {Xi , vii} to P with Xj integers, as asserted. 

Due to Theorem 1, it is innocuous to add the constraint that the xij are to be integers 
in P: there will always be solutions to P for which the Xij are integers. This fact is noted 
because existing algorithms for "'all- variable" integer programming problems are considerably 
better than those in which only a subset of the variables are constrained to be integers. 


1The matrix of constraint coefficients of a transportation problem Q belongs to a class of 
matrices which have the "unimodular property" and, as such, all vertices of the constraint set 
of Q have integral coordinates. This property also becomes evident inusing the simplex method 
for solving transportation problems in that only addition and subtraction are required in the 
successive pivotal operations. 
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THEOREM 2: Let P be the program P with integer constraints ignored. 1%; ; Vij! 


is a solution to P only if x,. = m; jij 


ij ij" 


Proof: We consider program P. If Vij = 0 then Xij = 0, and thus Xj =m,. y.,;. Other- 
wise, suppose Vij ij > 0 and x, ij <m; ii ij’ Then FF can be decreased without violating the con- 
straints of P but with a decrease - the value of the objective function (1). Thus Xj = mij Vij 
for any solution {X;, ; Vijs to P. 


Theorem 2, we shall see, iis the key to the approximation technique. 
THEOREM 3: Any solution {xj - vijt to P must also solve the program 


R: Minimize AR&) = D i=; C4; Xj 


subject to the conditions 


for Vij 


: fe) 
. + bij for Vij 


We call R a local description of P ina neighborhood of {xij ' vii! 


Proof: We repeat that since R has the same analytic form as M, and is thus a trans- 
portation problem, it has solutions in nonnegative integers [5]. 


Suppose that the above statement is false, i.e., that {xij} is not a solution to Q, but 
that {x4} is an integer solution. We define 


1 if xj; > 0. 


We have, by assumption, Ap = Ap &') - AR &®) <0 and Ap=A &',y')-A «°, y°) =0. 


Now, consider, in — the different contributions made to AR and Ap: 
1. x° =x, Then me 


ij ij’ ij = Vij and no contribution to either difference is made. 


ij ° 
oO 


(a) Xi . Then ai jXij + b.. is contributed to Ap and @ij + bj i) Xj to Ap- 


<x 


1) 


sac ee ; _— 
4j =1- Then a5 0, xij) is contributed to Ap and aj i xij) to Ap: 


Then -Ai jx} > bij is contributed to Ap and -5 5 Xi} to AR: 


Then -a i - x5 ;) is contributed to Ap and “aij Oxi - x5) to 4p: 


Add all contributions to Ap and AR: respectively, and consider the difference 
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- aij Xj , +b; )- >. (a; xj) + by) 


x; ;=0 


o.. 
xjj=1 


(bj; x4) +B) - >. bj, = 0 
Xj =0 


oO 
x4j=1 


But, by assumption, 4p - Ap > 0, and we have reached a contradiction. Thus {x35} is indeed 
a solution to R. 

In intuitive terms this result can be restated: A "'global"’ optimum 1x3) y: i of P is 
always a "local" optimum of P. Also, we point out in passing, Theorem 3 can be used to prove 
Theorem 1 since R, like Q, has solutions in integers. 


THEOREM 4: Any solution {xP vii} to P is a vertex of the polyhedral convex con- 
straint set S of P (the program P with integer constraints ignored). 


Proof: S is defined by 


-1, mij 


for i=1,...,m and j=1, 


A vertex of S is, by definition, the solution to a set of 2mn linearly independent equations 
taken from (6) and (7) as follows: All equations (6) belong to the set; the rest of the set con- 
sists of inequalities of (7) considered as equations and so chosen that the solution satisfies all 
remaining inequalities [8]. It is worthwhile to point out that this definition is equivalent to the 
perhaps more familiar definition of extreme point. 


In an optimal solution 1x4) , vii} to P at most m+n- 1 of the variables Vii are equal to 


1, and the others are equal to zero. For consider program R of Theorem 3, a local descrip- 
tion of P ina neighborhood of {xj , Vii}: There isa aenae” to R, {x iP which has at most 
m+n- 1 positive ce ; and hence at most m +n - 1 of the Vii are outlive in P. 


©} to P we a 


Therefore, “a an optimal solution 1x; ; Vij! 
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m+n _ equations 6) are satisfied; 


mn of the 2mn inequalities Vij = 6, “Vij = -1 are satisfied as equations; 


mn-m-n+l (at least) of the mn inequalities mij ij - x,.2 0 are satisfied as 


ij 
equations; and 


mn-m-n+l (at least) of the mn inequalities xij = 0 are satisfied as equations. 


Thus, at least 3mn- m-n+2 of the relations (6) and (7) are satisfied as equations for an 
optimal solution to P. From these we choose the 2mn equations: 


". oO 7 
~j Xij *.% 
(8a) any m+n-1 equations of (6) 
a a 
ti xij = dj 
Oo _ 7 
Vij = 0 
(8b) all mny.. 
cos J 
Yi = 
(8c) xi; = 0 } all mn-m-n+1l xe, which correspond to (i,j) for which Vii = 0. 


The matrix associated with the x) ; Vii in (8) isa 2mn x 2mn square and non-singular 


matrix. To see this we partition the matrix by permuting rows and columns to get: 


m+n-1 mn-m-n+l m+n-l = mn-m-nil 




















- 
A 0 0 | (ea) m+n-1 

0 I (8b) mn 
0 I 0 (8c) mn-m-n+1 

















where the columns of Cy correspond to the m+n-1 Vii which are positive (if less than 
m+n-1 are positive, then some Vii which are zero are included); the columns of C, corre- 
spond to the other VF; (which are all zero); the columns of C, correspond to those x2. 


1) 
for which Vij are in Cy; and the columns of Cy correspond to the other x3 (which are 


all zero). 


The rows within (8a), (8b), and (8c), respectively, are easily seen to be independent. 
The rows within (8a) and (8b), (8a) and (8c), and (8b) and (8c), respectively, are also easily 
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seen to be independent. Finally, any row of (8a) cannot be written as a linear combination of 
rows in (8b) and (8c) since each row of the square submatrix A contains a nonzero entry. The 
corresponding statements for rows of (8b) and (8c) are trivial. Thus the matrix of (8) is non- 











E singular. This completes the proof. 
4 This theorem may indicate the road to a modification and simplification of Gomory's 
fi technique to finding optimal integer solutions to fixed cost transportation problems, since these 
a solutions must actually be located at vertices of the original constraint set with integer con- 
4 straints ignored. Conceivably, a completely new approach to this very special kind of integer 
linear program may derive. 
APPROXIMATION PROCEDURE 
4 In Theorem 2, we found that for any optimal solution { Xj ' Vi) to P - the program P 
4 with integer constraints ignored - Xij = mij Vij . Thus we can simplify our search for an 
optimal solution to P by imposing the added constraint xij = MG; Vij Substituting Xii/ ij for 
ij in P we get an equivalent program to P. 

e / b 

a si ene 2 5s nee ij 

’ P': Minimize Ap: &) = 5; . + m,, Xj 

a subject to the conditions 
i. oF xij = es 25 Xj = di and xij = 0, 

P Ee , 

4 for i= 1,...,M and j=1,...,8%. A solution (Xj! to P’ automatically results in a solu- 
Ee tion {x.., y:.} to P by taking 

at 1) 1) 

ret 

e i ace Nl ——— 

7 Xj = Xj and Vij = *ij/ ij 

4 and it is easy to verify that the objective functions of P’ and P have the same optimal value, 
‘ i.€., Apr Ge’) = A &,¥). 

g But notice, now, that P' has the same analytic form as M, i.e., P' is a linear trans- 
E portation problem which, as was pointed out above, can be solved by any one of many highly 

: efficient methods of computation. It is interesting to see the relationship between the cost 

\ functions along route (i,j) of programs P and P'. This relationship is given in Figure 2. 

2 Thus the cost C55 &5;) resulting from a shipment of Xij units along route (i,j) in P 

3 

3 C;j&;) = 0 if x,; = 0 

q _ 

7 = aij Xj + bij if Xj >0 

Fi is transformed into the cost 

i 

2 . b.: 
ae P a 1) 
Sy "5 * ae 





when P' is derived from P. 
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, dij 
SLOPE aj, + —4 


'J 


at 








Figure 2 


The optimal solution {xij} to the transportation problem P’ can be taken in integers 
. * * * 
and can, thus, be easily modified into a feasible solution 1X; : Vij) to P (i.e., 1X; , Vij) 
satisfies all constraints of P). Namely, let ; 


. = 0 
™—” 5” 


* ' * . ’ 
xij xij and Vij = 1 if Xij > @. 


Now, suppose that {xii ? Vij! is an optimal solution to P, that { Xj , yijs is optimal for 
P, and {xj} is optimal for P’. Obviously, Apr (x') = \(X,y) = A(&«°, y°), since {xe : vii! isa 
feasible solution to P. Also, \(x°, y°) = A(*,y*) since xi; ; Vij! is a feasible solution for P 


and (xij Vii } is an optimal solution. Taken together these statements give us upper and lower 


bounds on the optimal value for P: 


Ap: &") = A GF) = 2 &, y9) = Aw). 


One might expect in: viii to be a reasonably good estimate and possibly an optimal 
Solution to P. While counter examples show that the latter expectation cannot be rigorously 
justified, the easily computed values A (,y) and A «", y) furnish evidence on the value of the 
approximation. 

To summarize: the approximation procedure is a sequence of three basic steps. 

1. Given P derive P’. 

2. Find an integer solution {x35} to P’ and its value Ap (x') by using some trans- 
portation problem algorithm. 

3. Determine {x;, ; Yijs and compute v«*, y’). 
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AN EXAMPLE 


Consider the fixed cost transportation problem given by the following data: 














20 50 


1. Formulate the derivative transportation problem P’. 





b.. 
ij 


“ii * min , di) 











20 


2. Solve the derivative transportation problem. 





10 0 
20 10 
20 0 

0 20 


' 
Xt 











3. Determine approximate solution to fixed-cost problem. 
The approximate solution to the fixed-cost problem is taken to be the distribution 
schedule {x;,} above. The total transportation cost is 


d(x*, y*) = 360 = \*. 


Finally, it is known that the cost 2 of the true least cost distribution must be between 
the values Ay, A* 


321.7 =A = 360. 


SOME COMPUTATIONAL EXPERIENCE 


Three "real" fixed-cost transportation problems having 8 reporting activities with 
excesses and 12 with deficiencies are solved using the approximation technique described 





50 


earlier in this paper. In each case, the proportional costs a;; are the same, and the quantities 
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of excesses and deficiencies at each of the reporting points are also the same. However, the 


fixed costs bij are increased by ten in each successive problem in order to observe the 


behavior of the approximation method as fixed costs increase. For the sake of comparison, 


the linear transportation problem without fixed costs is also solved. 

















Excesses Deficiencies 
Reporting Points e, Reporting Points d, 
(1) Scotia 15 (1) Boston 20 
(2) Philadelphia 20 (2) New London 15 
(3) Norfolk 45 (3) Newport 20 
(4) Charleston 35 (4) Mechanicsburg 15 
(5) Great Lakes 25 (5) Jacksonville 5 
(6) Long Beach 35 (6) Key West 20 
(7) Oakland 10 (7) New Orleans 30 
(8) Puget Sound 25 (8) San Diego 10 
ye = 210 (9) Port Hueneme 35 
: (10) El Toro 25 
(11) Port Chicago 10 
(12) Clearfield 5 
d. = 210 








mY 


J 








The proportional costs a; 


All Kinds' " [3]. 


J 


. are taken from Table S, 


"Transportation Rates for 'Freight - 





ij’ 





.69 


1.01 


1.05 


1.94 


1.61 


5.29 


5.29 


.64 -71 

75 .88 
1.06 1.08 
1.50 1.56 
1.40 1.61 
5.94 6.08 
5.94 6.08 
6.08 6.08 


.79 


99 


.64 


1.33 


5.29 


5.29 


1.70 


1.50 


1.22 


1.98 


1.68 


5.96 


5.96 


5.96 


2.83 


2.37 


1.98 


2.83 


6.77 


6.45 


2.02 5.64 5.94 


2.26 5.64 5.85 


1.66 5.64 5.91 


1.36 6.99 6.99 


1.54 4.26 4.26 


5.08 0.31 0.21 


5.08 0.55 0.35 


5.08 2.43 2.30 


5.94 


5.62 


5.62 


6.99 


4.26 


0.17 


0.40 


2.33 


5.94 


5.85 


6.99 


4.26 


0.31 


1.81 


7.60 


4.54 


4.54 


3.68 


2.99 


1.53 


1.53 


2.50 
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The fixed costs dij are between 10 and 20 inclusive and are taken from a table of 
random numbers. 





11 16 18 17 10 20 17 13 15 12 14 14 
14 17 17 13 15 13 16 11 20 11 15 10 
12 13 20 17 13 15 16 13 12 13 10 18 
..): 16 19 16 11 15 12 18 12 18 13 13 14 

J 19 18 15 16 12 14 20 19 11 17 16 18 
13 20 320 17 15 12 14 11 12 19 15 16 
11 12 15 10 17 11 11 16 10 18 17 12 
17 10 20 12 17 20 16 15 10 12 16 18 











The solutions to the four redistribution problems follow. 


1. Linear transportation problem with costs per route of aij Xij if Xij is shipped. 





15 15 

15 5 20 

20 15 5 5) 45 
20 15 35 

10 10 5 25 

10 25 35 

10 10 

15 10 25 


4 d.: 20 15 20 15 5 20 30 10 35 25 10 H) 
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True least cost: $266.70. 17 routes used. 


2. Fixed-cost transportation problem with costs per route of 


ie ya cia 
oO 


if a. =< 


aij Xij + bij if Xj >0O is shipped. 


see Na aa 




















d.: 20 15 20 15 5 20 30 10 35 25 10 5 


Total Cost: * = $504.55 
A» = $451.05 < true least cost 
d = $504.55 = A*. 
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If used transportation problem solution (1) as approximation to this fixed-cost 
transportation problem, total cost would be: $530.70. 


16 routes used. 


3. Fixed-cost transportation problem with costs per route of 


0 if Xij = § 


aj 


j*iyj * 5; + 10) if Xj > 0 is shipped. 











20 15 20 15 


Total Cost: »* = $645.65 


A, = $573.05 < true least cost 
d =< $645.65 = a*. 


If used transportation solution (1) as approximation to this fixed-cost transportation 
problem, total cost would be: $700.70. 


15 routes used. 


4. Fixed-cost transportation problem with cost per route of 


0 if x4; = 0 


aij Xj + 5; + 20) if Xj > 0 is shipped. 








20 15 20 15 20 30 


Total Cost: A* = $797.55 


A» = $694.80 < true least cost 
A = $797.55 = * 











FIXED-COST TRANSPORTATION PROBLEMS 


If used transportation problem solution (1) as approximation to this fixed-cost 
transportation problem, total cost would be: $870.70. 


15 routes used. 


This sequence of problems submits the approximation technique to a severe test due to 
the highly degenerate nature of the excesses and deficiencies. Despite this fact, the technique 
seems to behave properly; that is, as fixed costs become higher, fewer routes are used, thus 
avoiding more of the fixed costs. Judging from a small sample of actual redistribution prob- 
lems which are faced by the Navy Supply System, it seems that degeneracy is a phenomenon 
not often encountered unless artificially introduced. On the other hand, the lower bound A, on 
the true least cost of the optimal distribution for the fixed-cost problem, is well below the cost 
of distribution resulting from the (proposed) approximate distribution solution. However, A, 
is also well below the (unknown) true least cost A and is, as such, a rather misleading lower 
bound. Unfortunately, theory in support of this last statement has not been carried through so 
that, at present, no improvement on this lower bound is forthcoming. 


CONCLUSION 

We feel that the approximization technique herein described is very powerful and that 
under fairly general conditions will prove to yield close to optimal solutions to fixed-cost 
transportation problems. We have found, however, a number of techniques designed to improve 
the result obtained through use of the approximation procedure. As a rule, these have led to 
slight improvements; but we have been able to derive no theoretical proofs as to the generality 
of their applicability. 
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A CLASS OF DYNAMIC GAMEs*T? 


William E. Pruitt 


Applied Mathematics and Statistics Laboratories 
Stanford University, Stanford, California 


INTRODUCTION 

There is a large class of business and military situations in which decisions must be 
made on the basis of certain observed random variables. When these situations are described 
as games they often bear a striking resemblance to poker. A fixed-point method for the solu- 
tions of games of this type was proposed by Karlin and Restrepo in [1]. Several illustrations 
of its application to poker models were given in [1] and [2]. One of their examples, known as 
"le Her," indicated the applicability of the method to games in which there is more than one 
relevant random variable observed by each player. It is the purpose of this paper to examine 
some further games in which more than one random variable is observed by each player and 
some of the random variables may be observed by both players. Since the primary interest is 
in the method of solution, poker models have been chosen to illustrate the method because of 
their familiarity. It is hoped that solutions of some tactical problems of interest will be forth- 
coming. 

In the games considered here the strategies for Players I and II will be represented by 
the vectors 


Y (x! y”) 


(oy (x}, y?), iat Pm (Xs y’)») 
and 
(yt, x?) 


C¥ily?, x4). Ha, x’) ; 


1 and y2 are random vectors. x! denotes Player I's hand at the 


time he must make his decision as to his course of action; y- is the part of Player II's hand 
that I is allowed to see at the time I must make this decision; y} is Player II's hand at the 
time he must decide on his course of action and x? the part of I's hand II is allowed to see at 
this time. 

Let x denote Player I's entire hand, and y Player II's entire hand. Then the payoff in 
these games will have the form 


respectively, where x) x2, y 


K(v,y) = JJ P(x, y; o(xt, y?), v(t, xa F(x) dG(y) . 





*Manuscript received 24 January 1960. 
Prepared under Contract Nonr-225(28) (NR-047-019). 
{Written as Technical Report No. 19, December 4, 1959, Applied Mathematics and Statistics 
Laboratories, Stanford University, Stanford, California. 
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It will be possible to represent this in two different forms as follows: 


n 
K(v,¥) = el) +) SS ey (94, 4,07, x) AG (9) dF (x) 
j=l ' 


m 
K(9, 4) = eg) +>) SS cy Wh, x', 97) oj (xt, 9) A Fy (x!) aGa(y’) , 
i=l . 


where Cy is independent of y, the cy. are also independent of Y but are the coefficients of 


the Yi» and analogous remarks apply to Co and the Coe Fy is the marginal distribution of 
i 


x Fo the marginal distribution of x2, etc. In this paper, in the interest of simplicity, the 


cards are represented by values in the unit interval, each such value being considered equally 
likely. Also the various cards are independent and thus Fi; Fo, G,, and Go are uniform dis- 
tributions over unit cubes of the appropriate dimension. 

In order for a pair of strategies 9’, y* to be optimal, it is necessary and sufficient 
that 


7) 


m 
K(o", ¥") = max dep") +) Sf co Ws x1, v9) 9, (x1, y?) Fy (x) dG (y) 
i=l ; 


c n 
K(o",¥) . aaa Diss ey (",y', x*) ¥; (y", x”) dG, (y") dFy(x) 
j=l ; 


IE 
Pie 


ea 8 


Soe 


where the maximum and minimum are taken over all possible strategies for the first and 


second players respectively, which will be all possible ~ and y vectors subject to certain 
constraints. 


y Pelee 
oo) 


The basic technique in the fixed-point method is to use any information available as to 
the form of one of the optimal strategies (which may come from intuitive considerations, 
knowledge of the solution of similar games, examination of the c 0 and c ¥ functions, or 

i j 


possibly from experience in playing the game) to obtain some information about the c Y. or 


Co as the case may be, which will then give some information about the form of the other 
i 


optimal strategy. This process is then iterated until the optimal strategies are obtained. 

Later in this article, three continuous modifications of the two-person version of 
Baccarat are considered. The discrete version has been solved by Kemeny and Snell in [4]. 
Their description of the game is as follows: 

"The game to be considered is a card game with the special feature that 

hands are evaluated modulo 10. A hand will consist of two or three cards. 

Each card from ace through 9 is worth its face value, while each card from 10 

through king is worth 10 points (and hence 0 modulo 10). Thus a hand consist- 

ing of an 8 and a 3 gives a count of 1. A count of 9 is the best possible hand. 
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The banker serves as dealer, dealing his opponent and himself two cards 
each, face down. If either man has a count of 8 or 9, he announces this fact, 
and hands are compared at once. 

If neither player has a count of 8 or 9, then the player has the option of 
taking an additional card. If he elects to do so, this card is dealt him face up. 
The banker may then decide to take an additional card, if he wishes one. The 
hands are then compared. 

When hands are compared, the higher count (modulo 10) wins. If the two 
men have equal counts, the game is declared a draw. 

Since nonplayers may bet on the player's hand, the player's strategy is 
f restricted by the rules of the game. He must draw if he has a count of 4 or 
less, and is not allowed to draw on a count of 6 or more. His only free choice, 
arises whenhe has a countof5. The banker is free to make his own decisions." 
In the versions considered here the restrictions on 8 and 9 are removed, and the player 
ly * is allowed complete freedom in his choice of strategy. Furthermore the first two cards are 
S- ‘ replaced by one since this only changes the distribution involved. In the first model the 
player's additional card (if any) is dealt face down. This is done to give a very simple example 
as an introduction to the method. The second model is the same with this card dealt face up. 
In the third model betting is introduced and again the player's additional card is dealt face 
down for simplicity. In the first two models there is a single critical number for each player, 
which is explicitly obtained in the analysis, such that his optimal strategy is to draw if his card 
is smaller than the critical number and refuse to draw if his card is larger than the critical 
number. Of course, in the case of the banker (Player II) the critical number is a function of 
whether or not the player (Player I) drew and, in the second version, of what card he drew. 


Lo SHEERS Seekers 


4 

The optimal strategies in the third model depend on the ratio of the size of bet to the size of 

Ri ante and bring in the concepts of folding and bluffing in addition to the above. 

Also later in this article, a continuous version of stud poker is analyzed. The descrip- 
By tion is as follows: After an ante each player is dealt two cards, one face up and one face down. 
- Player I then has the option of betting or folding. If he bets, Player II then has the option of 

ie folding or seeing the bet. In the latter case, each player is dealt one add ional card and the 


player with the higher maximum card wins the pot. The form of the optimal strategies depends 
on the ratio of the size of bet to the size of ante. If this is less than 2, for example, there isa 
critical number k, again explicitly obtained, such that: 

If II's "up card" is both larger than k and larger than I's maximum card, then Player I 
should bet with some probability between zero and one which depends on the size of II's "up 
card"; otherwise he should bet with probability one. Player II should fold if his maximum card 
is smaller than I's "up card" and I's "up card" is larger than k; he should bet with some prob- 
ability between zero and one, depending on the size of his "up card," if his "up card" is greater 
than his "down card," greater than k, and greater than I's "up card"; in all other cases he 
should bet with probability one. 

An interesting feature of the optimal strategies of this game in all cases is that I's 


Sema a 


} strategy depends only on the maximum of his two cards, while II's strategy is a function of his 
two cards individually. 
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BACCARAT 

Model I - This game is played as follows: Each player receives a card. After seeing 
his card, Player I has the option of drawing another card. Then Player II, knowing the value 
of his card and whether I has drawn or not, has the option of drawing another card. The player 
who has the greater sum, reduced modulo one, wins one unit from the other. 

The hands are < X4, Xg > for Player I and < ¥4> V9 > for Player II, X14» X9, Vy» Vo 
being independent uniformly distributed random variables. Of course, Xo and Yo will not 
necessarily come into the picture every time the game is played and in any case are unknown 
initially. The strategies are completely described by: 





















~ (x,) = probability I will draw if his card is Xy 


l-@ (x,) probability I will not draw if his card is Xy 
val (y,) = probability I will draw if his card is y, andI drew 
1-4 (y,) 
¥o(¥,) 


1- Yo (y,) = probability II will not draw if his card is yy and I did not draw, 


probability II will not draw if his card is Vy and I drew 


probability II will draw if his card is yy and I did not draw f 


where 0 = y(x,)=1,0= 4, (y,)) =1, and 0=Yoly,)=1. 


Let ) 
1 if r-([r|] > s- [s] 
L(r,s) =< 0 if r-[r|= s-[s| 
|-1 if r-[r|] < s-[s], 





where |r| denotes the greatest integer contained in r. Then the payoff is given by 


; £2 J 3 
K(¢9,¥) = J J J J [oy (x4) ¥, (yy) L(x, +Xo, ¥4+Yo) * Y (x,) {1- ¥1 yp} L(x} +Xp, y,) 
+ {1-9 (x,)} Yo (¥y) L(x,, ¥4+Yo) +{1l-9 (x,)} {1- ¥o (v4) L(x,, yy) | dx, dx» dy, dy 
; 
= Jf leg tt - Hy pi - 29) +11 - @ pt Hp) 2x, - 1) ie 
+ {1 - o(x)} {1 - Vy (vp)} LOX, y,)] dx, dy, x 


i 1 
c, (¢) + J “WY (9, yy) ¥1 (yy) dy; + | “WY (9, y,) ¥ 9 (4) dy, 


1 
= Co (¥) + J col¥, X1) ? (x,) dx, . 
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, order to minimize K(o’, yY), II should choose y 1 (9) to be one or zero according as 
%y, (y" . y,) is negative or positive and similarly for y 9). In the same way I will choose 


vind to be one or zero according as c wv . x4) is positive or negative. 


Yr a Since 
a x 1 x 
a vw » 4) = (2y, - )) t g (x,) dx, , 
a 1 
ba 1 yy < 5 
x 
; y 1 (y,) ™ 
Ay 1 
‘ 0 Vy >=. 
This part of II's strategy is intuitively clear, since if I draws his sum (mod 1) will be uni- 
formly distributed regardless of the first card. 
f Now 
* i * ‘4 * 1 * 
+ ’ y,) = J (2x, = 1) {1 =¢ (x,)} dx, ¥ J, {1 = (x,)} ax, im J {1 -@ (x,)} dx, 
1 
and 
* oa 1 Jj * bs * 1 * 
) Cy (y , x) an + (1 - 2x,) Jo V9 (y,) dy, * ‘ {1 = vo(y,)i dy, + I {1 = vo (¥,)} dy, . 
1 
The first is monotonic nondecreasing in yy and the second monotonic decreasing in Xy- 
Therefore 
X, <a 
* 
YQ (x,) = 
|o xX, > a@ 
where 0 < a < 1. Substituting this in cy (0, y,) yields 
2 
-(1- a)? yy < a 
* 
cy 2? »¥y) = ' 
a -1 - a + 2y Yi > a 


and thus 
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* 
There remains only the determination of @ so that ey ,®) =O0or 


2 2 
1 l+a l+av _ 
-24 1-20) (1427) 41 - La2" 


which is satisfied for @ = 0.567. 

The value of the game is given by K(y’, ¥*) = - 0.0112. 

In summary, the optimal strategies are: I should draw whenever his card is less than 
0.567. Then II should draw whenever his card is less than 0.5 if I drew, and whenever his card 
is less than 0.661 if I did not draw. 

Model II - This game is the same as the first one with the single modification that if I 
draws an additional card, this card is dealt face up so that II knows the size of it before decid- 
ing whether or not to draw. 

The structure of the hands will be the same as above and the only difference in strate- 
gies is that now y, (y,) is replaced by 


vy (y,, Xo) = probability II will draw if his card is y, andI drew Xp. 


The payoff is given by 


aS 2 2 
K(9,¥) = J ‘4 J [9 (x4) ¥ 4 (Vy .Xq) L(x, + XQ, 94 +Vq) + P(X)LL- Hy (Vy, Xq)F L(x] + Xp, Yy) 


Oo .O 
* {1 -@ (x,)} %o(y4) L(x,, yy +Yo) +{1 -@Q (x,)} {1 -%o(y,)} L(x,, y,)) ax, AX, dy, dyo 


11} 
J J, ui Y (x;) ¥ 1 (¥y, X9) (2x, + 2Xp - 1) dx, aX dy; 


- & 
J a Y (x,) vy (y;, Xo) (2x, + 2X5 - 3) dx, dX» dy, 
2 


- 42 3 
J J J [¢ (x,) {1-¥ 1 (yy, Xo) L(x, +Xp, 4) +{1- 0(x,)} ¥ 9 (y,) (2x, - 1) 


{1 - 9(x,)} {1 -y 9 (¥4)} L(x,, y,) | dx, Ax, dy, . 
* 
It will now be assumed that ” (x4) has the same form as in the first model, i.e., 


l xX, <a 


ok 
Y (x,) = 
0 xX, > a 
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* * 
where 0< a <1. After obtaining y , it will be clear that cy »X4) is monotonic decreasing 
* 
and this will justify the assumption about ¢ 
Now 


2 
[a + 20x, -44+4+2+42y) - 2x, Yy~Xg<a-1 


| 2 
a + 20 X5 - 2a a-1<y,- X%)<0 


(y" ) 
c, W Fark *< 
v1 as 2 


| o 


+ 20X9 - 2a + 2y, - 2X5 O<yy-X%, <a 
| 2 
| @ + 20 Xo a < yy - Xp 


and thus vi, Xo) has the form 








KX 


2 @ 
ae 
2 


ae- 


* 
where the shaded area is the region where y 1 (y;, Xo) = 1, the nonshaded area the region 


i ie | 

* : 

where ¥ 1 (y,; Xo) = 0. It is interesting to note here that J J Vily, Xo) dy, dX» = 1/2 which 
oO 


indicates that the probability of II drawing when I has drawn is 1/2 as it was in the original 
model. 


+ 
Since cy (y , y,) is the same as in the previous case 
2 


F 


¥o(¥,) = 





W. E, PRUITT 


Substituting y* into K(9,¥) yields 


a 2 











x 2 
- xf +) - 2x, + ax, - a?x, - 4 S48 x <@ 
(y* )= 2 a a2 _a? 202 a<x else? 
2 2 
joa ~ Oey +m > Oxy + - #,% -t < Xy- 


* 
This is monotone decreasing and thus g~ does have the assumed form. Nowa must be 
determined so that c,(¥ ,®) = 0, or 


2 
«wv « £2 Rs Bad 
4 3 3 


which is satisfied for @ = 0.580. In this case the critical number for J, i.e., (1 + a”)/2, is 
0.668. The value of the game is - 0.0200. 


Model III - This is the first game with the addition of betting. Both players initially 
ante b units. After seeing his first card, I may either fold (forfeiting his ante to Il), or bet 
a- b units. If he bets, he may draw an additional card (face down) without further payment if 
he so desires, Then II has the option of folding or seeing the bet, and of drawing an additional 
card if he sees the bet. In the event that II sees the bet, the hands are compared as in the first 
game, the player with the higher hand winning the pot. 

The hands have the same structure as before, but now the strategies have additional 
components. The strategies in this case are described by: 


aT (x4) = probability I will fold if his card is x, 
Pq (x) = probability I will bet and draw if his card is x, 
Pa (x4) = probability I will bet and not draw if his card is Xy 
y 11 Vy) = probability II will fold if his card is yy and I drew 
¥o4 (y,) = probability II will see and draw if his card is yy and I drew 


Ya (y,) = probability II will see and not draw if his card is y, and I drew 


y 12 (y,) probability I will fold if his card is Vy and I did not draw 
y 22 (y,) = probability II will see and draw if his card is Vy and I did not draw 


yY 32 (y,) = probability II will see and not draw if his card is Vy and I did not draw 


where these are subject to the restrictions ¢;(x,) = 0, Vij (y,) =0, 
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. Q; (x4) 1, : Yi; (Vy) 21, j=1,2. 


The payoff is given by 


K(9,¥) = f f f f [- by (xy) + by (xy) Hy (¥y) + A Fy (X) Hoy (4) LOX, + Xq, ¥y + YQ) 
+ Ao (X41) ¥3y (¥y) L(x, + Xo, Yq) + DY (xy) Y 49 (¥y) 
+ Pg (X1) Yoo (yy) L(x, ¥y + Yq) + 2g (Xy) ¥ 39 (¥4) L(x, yy)] ax, dxy dy, dyy 
11 
. J J [- b py (xy) +B Gy (xq) 4 14 (Vy) + 2 Gy (Xy) Hyy (Vy) (1 - 294) + BGG (Ky) Hy (94) 
+ 


a Ps (x) y 22 (y,) (2 xy -l)+a %3 (x,) 439 (y,) L(x,, y,)] dx, dy, . 


* 
In this case, in order for I to maximize K(v,¥°), he should choose 9; (x,) = 1 for all 


* * 
xX, Such that Cy (y »X4) > €, (y »X,) , j = i. Player II will achieve his minimization of 
i j 
* 
K(¢ ,¥) in a similar way except that he will be interested in the smallest cy (o’, y,) for 
ij 


fixed j, and he must do this separately for each j. 


The optimal strategies will, of course, depend on a and b and even the form of the 
optimal strategies depends on the ratio a/b. With the exception of small values of a/b (when 
the optimal strategies are the same as in Model J), the set of optimal strategies for at least one 
of the two players contains more than one member and the entire set of optimal strategies will 
be given here. 


: * sd * * = * 
Since ey? » 94) =b J Po (x1) ax, ’ “0, 9 »Y4) = 0, and oa »Yy) 


= a(1- 2y,) J 93 (x) dx, , 


1 yy <3 0 Vy <3 
+ * * 
¥ 44 (9) 2 0; Yo (¥y) = ; % 34 (¥y) = 
0 Vy >¢ 1 Vy >$ 


* * * . * 
Now co X)2-b, Co, ,X,) = - a/4, and Co, ¥ »Xy) = b J Yyo(y,) dy, 


+ a(2 x1 - 1) J Yoo (y,) dy, +a J V30 (y,) L(x,,Y}) dy,- This indicates that the first range 


to consider for a/b is 1 < (a/b) < 4. Then, since Cy (", x,) is nondecreasing 
3 
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1 Xy <a 
* « * 
91 (x,) = 0; Po (x4) = ; $3 (x4) = 
0 X; >a 


Substituting this into K(y,~) yields 


a(1 - a) ¥, <@ 

* * a * 

Cy (@ »¥y) = b(l-a); cy (yo ,yy) FS aa(l-a); cy (Y , yy) = 
12 22 32 og 

a(1- 2y, + @) y,>a 


If a is such that aQ@ < b, then 


2 
1 y,< ++ 


Vio (¥y) = 0; V9 (¥,) = ; ¥30 (V4) =< 





2 


a 
There remains the requirement that Cy (¥ ,a) =-a/4, or 
3 


2 
site ae -a(t- 


which is satisfied by a@ = 0.567 as indicated above. This then is the solution for 1 < (a/b) 
< (1/0.567) = 1.764. 


If a is such that aa = b, then 


Vi9(¥y) + ¥50 (V4) ” 
1 y,> 


* * 
and ¥ 12 (Y,) . ¥o9(¥,) are arbitrary for y, < (1 + a2)/2 except that they add to one. There 
* 
is yet the requirement Cy (Y , b/a) = -a/4, or 
3 
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er 140 
i . 2 , 2 . 
b J Vig(y)) dy, +a (22- 1) a —- S ¥i2 (vp) ay, /- a\1-—* /=-F 
ae 3a b? 
al (a-b) J ¥y9(y;) dy, = —- b- — 
4 a 
_ 3a° - 4a - 4b? 





, * 
> J ¥42(¥) dy, = 3 
4a“ (a - b) 


x 
and this restriction must then be put on y 12 (¥,) (and therefore on Vo0(¥,) ). It is possible to 
do this if and only if a/b is such that 





2 
1 + be 
p< 3a°-4a%-4p> a? 
4a7(a - b) 2 
which is true for 1.764 =< (a/b) =< 2.920. 
If a is such that aa > b, then 
1 y, < 8B 0 yy, <8 
* * * 
¥ 429) = ; Yo9(¥4) = 0; Y 39 (¥4) = . 
0 yy, >B 1 yy >B8 


* 
This implies that c (y",x ) is constant for x, < 8. But comparison of c (~ ,y,) and 
25 1 1 Vio 1 
* * * 4 
“Yao? »Y4) shows that a < £8 and thus Co, »X4) = Co, »X4) for 0 = x; = B. This 


will be true if and only if 


or B = 3a 


bB - a(1- 8) = - ~ 4(a+b)- 


a 

4 

The equality of c y*,x ) and c @*,x ) for 0 =< x, < B gives more freedom in the 
%3 1 Po 1 1 


* 
choice of ~ than was originally supposed, i.e., 


arbitrary Xy < B arbitrary xy < B 


* * x 
P 4 (x) =0; P 9 (X}) = ; 0 3 (x) - 


In order to insure that these strategies are optimal, however, it is necessary to require 
* 
¢ Gi = 
vio ( Y >| ) 


c (yp) and c (o", 0) <c (y",0). These conditions are that 
¥ 32 ¥12 ¥ 92 
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2 2 
x » 
f Po (x) dx, » £3 ee and J Po (xy) (2x, - 1) dx, <- b(a_+ 4b) ; 


2(a +b)” 2(a +b)” 




















and this is always possible if a/b > 2.920, at least by taking 09 (x4) = 1 for 0 < Xy 


< (a2 + 2b2)/2(a+b)2 and zero elsewhere. 
In case a/b = 4, the analysis is very similar to the last case above except that now 

* * 
21 (x4) is also arbitrary for 0 < xX, = B = 3/5, and now the conditions on ~ are 


J (opm) + ogy) )dxy= 2, and f (v(x) + 90%) ) (2x, - 1) dxy = - 


=. 
25 25 





which are easily satisfied. 
* 
In case a/b » 4, the analysis is once again similar except that now Po (x4) = 0 and 


* * 
only 4 (x4) and Pg (x4) will be arbitrary on 0 < xX, = 8. Also £ has a different form since 


now the condition is c | (",x,) = ¢,, (¥,x,) for 0 = x, = 8 which yields bf - a(1 - 8) 
3 Y1 ‘ 








= -b or B = (a-b)/(a+b). The conditions c (o",B) =c¢ (*, B) and c (0*, 0) i 
. ¥12 32 ¥12 
= ty (y ,0) are now | 
22 
f ©, (x4) ax = fecey’ and f{ o. (x4) (2x - 1) dx, <= -b/y (a=)? ] 

oe ft \aeb a“ 8 i” a a+b) 


which are easily satisfied. 


* * * * 
In this case, since 9, =0,c Y , zc Q , Bec Y , = 0. Hence 
2 val y,) vor! y,) var! y,) 


* 7 * 
Vi (y,) . %91 (¥,) , and % 94 (¥) can be arbitrary in this case except for the condition 
* * 
“09 (y » X4) = Co, & »X4) or 


b f Vip) dy, +a J ¥31 (V4) (1 - 2y,) dy, -«<*: 


A summary of the optimal strategies in all cases follows. In cases I - IV, 


. 

1 W4< ; 0 yi, < 5 Ms 

* 4 * * iS 
441 (9) = 0; %o3 (¥,) = ; % 34 (yy) = oF 
1 1 

0 >= 1 = Fe 

yi > » Yi > 2 i 


£. 
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Hepa 


Tass 
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1<2< 1.764. 
b 





¥i0(¥)) = 0; 


* 
91 (x,) = 0; 


¥12(94) + V9 (V4) = 


* * 
4 (x,) =0; Po (xy) = 





Ss 
& 
: 
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2.920 <2<4, 
b 








4 (x,) =0; 09 (x4) = ’ 





















yy < 0.661 Y 
* a ; * 
Vo (¥,) = ; ¥39 (y,) "7 
0 y, > 0.661 ¥, > 0.661 
1 xy < 0.567 Xy < 0.567 
* * 
Po (x4) = ; Pg (x4) = 
0 Xy > 0.567 1 Xy > 0.567 
< 2.920. 
rc r 
2 2 2 2 
: a a +b 0 y,< a” +b 
1 2 1 2 
2a 2a 
* 
; ¥ 39 (Y,) a“ 
a” +b a? + b@ 
0 Yy > 9 1 Y4 > 9 
L 2a L 2a 
3 3 
x al = 
J ¥40(¥,) dy, = 3a 9 i. * 
4a“ (a - b) 
( b 
; 3 (x,) " ° 
0 b b 
ay >= 1 —y >= 
a L a 
3a vy < 3a 
4(a + b) 4(a +b) 
* - * 
; Yo9 (¥4) =0; ¥ 39 (Y¥4) =< 
% 3a Y,> 3a 
4(a +b) 4(a +b) 
cr 
3a 3a 
arb x, < —"— arb x, < —*"— 
1 4(a + b) : 1 4(a +b) 
; Pa (X,) = 
3a 3a 
0 x eo. 
el 4(a + b) 
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and satisfying 
a + 2b" 


f G5 (x,) ax, = = 
2(a +b) 


J op (x) (xy - 1) dx, = - P@+ 4D) 
2(a +b) 


Case IV. 


; ¥59 (V4) 20; Va9 (V4) = 
3 


1 ~~ = 
[ sea 


5 5 


* * ‘ 
Yy (x,) = < . Yo (x) = < : 0 (x,) = 
3 
0 ae 
1-5 


and satisfying 
9 


; a 
J P4 (x4) dx, + J Po (xy) dx, = on 


J v% (x) (2x, - 1) dx, + J 09 (x4) (2x, -® dx; <- . 


25 


CaseV. 25 4, 
b 


* a * 
% 11 (94) » 424 (94), and ¥3, (y,) are arbitrary except for satisfying 
. * 4 * 
bf 44 (yy) dy, + a f %51 (yy) (1 - 2y,) dy, = -b. 


% <7 
»? Yaatyy) 205 ¥39(%) = 


9 (x,) aa 


PAS Tao Sat Roce ae 
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and satisfying 


and 
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J of (mp ax, = (22) 


. b | -b 
J 9) (2x, - 1) dx, < - ‘ - - 


STUD POKER 


ante of 
Player 


In this game both players are dealt two cards, one face up and one face down, after an 
b units. Then I has the option of betting a-b units or folding. Following a bet by I, 
II has the option of seeing the bet or folding. If II sees, both players receive a third 


card and the player with the higher maximum card wins the pot. 


The hands are < X41, X9, Xg > for Player land < Yy> Yor ¥g > for Player II (x), y, are 


face down; Xo, Yo face up; X3, Y3 the third cards, if any). The strategies are completely 


described by 


9 (x, Xo, Yo) = probability I will bet holding < X41, Xq > if II's "up card" is Yo 
1-o(x,, Xo; Yo) = probability I will fold holding < X1,X_q > if II's "up card" is Yo 


7 (y,, Yo: Xo) 
1-y (y4, Yo» Xo) 


probability II will bet holding < Yy> Yo> if I's "up card" is Xo 


probability II will fold holding < Yy> Yg> if I's "up card" is Xo - 


Let 
1 r>s 
L(r, s) =< 0 r= s. 
—- r<8 


Then the payoff is given by 


K(9, 


Let é 


8 43 2 £3 

)= of - bil - —(X1,X9,Yo) f + DY(X1,X5,¥o) {1 - ¥( Xo) 
eeergre.. 0(X45X9,¥q) | 9 (X1,X9,¥o) {  (¥4sVo,Xq) } 

+a 9 (X4,X9,Yo) V(Y4,Y9sXq) L (max (x1,X5,Xg), max (Y4,Yo,¥9)) | dx, dx dx, dy, dy dys . 


* 
= max (xj, Xo) a max (Yj, Yo) and denote cy (gy 1 V4» Yq» Xo) by cy and 


* 
Co »X4,Xo5Vo) by C+ Then 


% — DN * 1 * 
-b J P (x,, Xo; Yo) dx, = an J Y (x1, Xo; Yo) dx, +a J x? ? (x, Xo, Yo) dx, Xo <1 
o U] 


x 








- * 2 x 1 * 
~& J ~ (x1, X95 Yo) ax, +a x J Y (x1, X9> Yo) ax, +a J x? Y (x1, Xo, Yo) dx, Xo >1 
L ° 9 
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* 2 gs 1 2.* 
2b-b fy (Y4,YooXq) dy, +aé f Y (Y¥4,YosXq) dy, - <! yy y (¥1,YosXo) dy, Yo< g 
oO 


V2 1 
* . * * 
2b-b fy (¥15¥q»Xq) dy, - aye J % (¥1,¥osXq) dy,-a f yey (¥y,Yo:%_) dy, Yo >é - 
oO Y9 





If a-b is sufficiently small, it seems reasonable that I should bet whenever § > Yo or 
whenever Yo < k where k is a parameter between zero and one. When  < Yo and Yo > k, 
I should bluff, i.e., bet with some probability between zero and one. With this in mind, let 














1 if yy <k_ or Yo < & 
* 

Y (x,, Xo, Yo) - 
? if yy >k and yo>& 


Then 





a 4 3 r 
obez-ar’. Yo < kK, Xy <0 , 
-bed-5 an, Yo >K, Yo <%y< Vy | 
72 2 2 s a3 
~b JP Gy %q9¥q) dx - b+ byg- ayy J (xy, %qr¥q) dm +5 - 5 92» 
oO 


Yo >k, Yo > X9; Yo > Vy 





. * : * 3 
-b J i (x1, Xo, Yo) ax, - b+ by» s ay? J Y (X15 Xo, Yo) ax, ” ay? (Y4-Yo) = * 331 ’ 
oO oO 
| 
| ¥g > Bs ¥g > %q> Yq < Vy ie 
| 4 
— ea 
-beSsF ax ‘ Xo > 1). Be 
L 3.3 2 2°" é 


The requirement on the size of a/b for this case now clearly stands out as a/b < 3. Then 
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(1 Yo <k, Xo <7 
Yo >k, 


< 
i) 
f 


j 1 Cm<F 

* 2 1 
Y (y,, Yo: Xo) = 4 1 

} 

‘ Xo > n, Xo > a 


Xo > 1; Xo < a 


where - b + (a/3) + (2/3) aa® =0 . The form of ¥” (¥4s Vos Xo) in the remaining regions is 
not yet clear. 
Since bluffing by I is expected when Yo > gs Yo > Kk, c. must be zero in this region, 
1, 
y 
rt 2,2) + 2 
2b - b J ¥ (Y15 Yo, Xo) dy, - aYo J ¥ (y;; Yo: Xo) dy, -a ‘ yy v7 (¥15Yo; Xo) dy, = 0 
2 
when Yo > k, Yo > Xq> Examination of cy in this region shows that it is a function of Yo 
x 
and X» only provided that y; < yy. If for some yo, x» this were positive, then y (yy, Yo Xo) 


= 0 for yo > K, Yo > X9, Vy < Yo and 


1 , 1 

; Pe 3 

C= B= b Jv (Vys¥qo%Xq) dyy- a J Y_H (Vy,¥q9%Xq) dy, = b+byy-F+F yp > 0, 
Yo Yo 

a contradiction. Thus cy < 0 for Yo > k, Yo > Xq, Vy < Yo- But cy for Yo > kK. Yo > Xa, 

¥y > Yo is strictly less than this (holding Yo and Xo fixed) and hence strictly less than zero. 

Thus 


a 
y (Y1»Yo.%q) = 1 Yo >k, Yo > Xo Yq > Yo: 
Substituting this back into ec. yields 


b- 2+ by, +2 yp 
? Yo >k, y > Xo - 
2 2 2 2 
b + ayo 





V2 
. 
V (¥4,¥qs Xo) dy, = 
oO 


* 
Any ¥ satisfying this would give an optimal strategy, but the simplest is 


a a 3 
» b- 3 + byy + 3.¥o 
v( X,) = 3 3 k x 
¥4°Yo> 2 3 ’ Yo > K,; Y9 > 2? Yo > Fg 
byy + ao 





To insure that this is less than one it is sufficient to require 


(1) b- 2-2 ak? <0, 
3. 3 
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* * 
The next step is to determine the y that maximizes against this y_ , hoping that it 
_ 
will be the original Y . Now 


fy-8+4ag5, Yo < X%g» Xo < a 


b+ bxy- S42 8°, Yo < Xg» Xg > M, Xp > Xy 


b+ bxp + axy (xy = Xp) - 2 +S xt, Yo < Xgs Xg > @, Xy > XQ 


bo S +i age, Xo < Vo < X1> Yo <k 


y y 
P22 « 22 2 

b+byy-bf W (¥4,¥o Xp) dy, + ax? | W (¥y,Vqs Xo) dy, + axy (x, - Yq) - 5+ x 
oO oO . 





Xo < Vg < X4> Yo >k 


Yo>&s Vg <k 


Yo>é, Yo >k. 





The first five expressions are positive and to insure the positivity of the sixth it is sufficient 
to require 


nee ; er, 
pe ae. 09 ONE) 


Sy 


(2) e-2.2@" =8. 
3 3 


Cee: 


Then the original o maximizes against y* . (1) and (2) imply k is the solution of 
b - (a/3) - (2/3) ak® = 0. 
The condition cy = 0 for Yo > # Yo > Xg, Yo > Vy yields 


Se 3 a 
= 4 b 
3 2° 3 

by» + ays 


by» > 





* 
Q (X1,%X9,Vo) = Yo > &> Yo >k, 


. + 
where the simplest version has been chosen as was the case with Y . This adds the condition 


3 .a ; 
B*>°e2° & kK=<Yyo=1. 
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‘“ a.3.a 
{(Yo) = bY5 "37273 - b. 


Then {(1) = 0, £(0) < 0, f'(1) < O imply f has exactly one root between zero and one. 
Thus {(Y¥5) = 0 fork ~ yp 1 if and only if 


3 3 2 
3 b 2 b . 6 b 9/b 3b 1 
f (k) = 0 <> bk- ak" = O<> > =k © (2) = « (8) 7a +e--7=0 
Let 
a 
g(x) = x 4% +o* 4° 


Then g(1) = 0, g'(1) = 0, g"(1) » 0, g(0) < O imply g(x) > 0 for 1/3 < x < 1 if and 
only if g (1/3) = 0 and g (1/3) = 1/27 > 0. 


Hence the condition is satisfied for 





te bei, 
3 a 
the range under consideration. 
In summary, if : < 3 
1 if Yo <k or Yo ae 
*( ) 
Op x } 9 = 
7 1 *2 V2 by, - 2 yp +2-b 
3 if Yo >k and yo > & 
bYo + ayo 
C 0 if X» >, Xy» >k 


a a 3 
o-3 eS 


byy + ayo 


v (y,, Y9> Xo) se 4 





if Y9 > k, Yo > Xg> Yo > Yq 





L 1 otherwise 


where k = 
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When . = 3, it seems reasonable that I should bet less often, namely when & > Yo 







and é > k. Accordingly let 


1 Yo<é,; é >k 














. 
Q (x,, Xo» Yo) = 
|? Yo > é or € <k. 





( y 
v2 
-b + bys - (b+an’) { ~ (X41, Xo; Yo) dx, - an? (n - yo) 
fe) 


Xo <Yo> Y9 >k 


’ Vo <Xq <j) %_>k 





Ks 29,7 x 
-b{ o (X15 Xp, Yo) dx,;-b+bk-an J{ o (x1, Xo, Yo) dx, 
oO oO 


alae dx, +2-85° k 
JO LP (Xp Vo) xy +S - SK, Xp << 
n 


k 


Kx 
. J Q (x,, X9; Yo) dx, 
Oo 


k . 
-b/ ¢@ (x1, X55 Yo) dx, - b + bk - an 
oO 


2 a_a.3 
-an w-Mer+ee Xp <K < yy, ¥9 <k 


k Xx 


; * “2 * 
-b{/ @ (x1, X95 Yo) dx, - b + bk + ax? J @ (x1, Xo, Yo) dx, 
Oo oO 


an 3 . 
+a f xp¢@ (X1,Xp, Yq) dx, +2-2k ‘ Xp > 1; Xo <k Pe 
Xo 3 3 
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Suppose that % = 0 when Yo > Xg» Vo > k, Yo > Y, aS was the case when a/b < 3. Then 


a a 3 
* oe ae 
Q (x1,Xo,Yo) = 





3 ’ Xo < Yo» Yo >K; Vo > X1> 
byy + aYo 


* 
where ¢ again is taken as a function of Yo alone for simplicity. The numerator may be 


written as § - b) (1 - Yo) + : (v9 - ¥3) = 0 and the fraction will be less than one provided 





(3) 2-b-2 ak? <0, 


* 
In the region Xo < cK, Yo < k , it is reasonable to suppose that WY will be one or 


* 
zero according as 7 > k or 7 < k, this corresponding to the behavior of ~ and also being 
a natural break point in ty as given above. This will be the case if 





k 
* - 7 * ’ Be. 
f(y ,k) =- ateeali Et Y (X1,Xq, Yq) dx, - b+bk “| ~2s = 0 (xX) <k, yo <k). 


* + 
This can be solved if k satisfies (3), for ¢ = 0 on this region makes f(y ,k) > 0 and 


* * * 
y = 1 onthis region makes f(y ,k) = 0. Supposing ~ to be chosen on this region so that 
* 


{(v ,k) = 0 yields 


0 if Xo > or 1 < k 
* . 
Y (Y4> Yo» Xo) =< ? if Yo > X9> Vo > k , Yo > Vy 
1 otherwise. 
Substituting this into c,, gives 
- 
b+bk - 24243, <&<k 
3 3 Yo g 
b+bk + ag 2(é -k) - 242 28, Yo <k <x »Xq <k 
3. 3 2 1 
» Ye Yo 


: * *x 2 3 
sis Mak (ys Yop %)ayy +a” J W (¥ysYq.%Xq) dy, +a i ~a9- 744 


ay 


Xq < V9 <Xy, Vg >k 

7) 

b+bx sat 2(é -x,) - 24223 S— > Xy >k 
2 ?-3*3* e> For"? 


b+bk - 24243. < Yo <k 
3 3 g Y9 


Re 


V2 
: * 

b+byy - (b+ay3) | W (¥4;¥o,%Xq) dy, - : += Yo, Yo>&5¥qg>k. 
oO 





ry 


See, eee 


Y 
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* 
Since ” is between zero and one when Yo > a4 Yo > a 
which implies 


¢ must be zero in this region 


. b+byo - ; +. y 
y (¥4 Yq» Xq) = 3 ’ Yo > X95 Yo > Vy> Vo >k. 
by» + ayy 





To make this positive, it is sufficient to require 


(4) b+ bk -2424% =0, 
3° 3 


* 
This determines ~” except in the region Yo < k, € < k and it agrees with the original 
* 
choice of 


The range of ; must now be divided once more. Let a satisfy 


2_p-4aa? 
3 


3 


= 0 


and £ satisfy 


ps obs be - 24% => 0eba - aa®=06 


‘ 4 


+3 > 09 =—4 < 4,822. 
b 


Consider first the range 3 = 4.822. To satisfy (3), it is sufficient to require 
. If k > a, then 


b+ bk - 24242 5 b+ ba +2 a3 => b+ bp - 24223 = 0 
3° 3 3 7 3 


* 
which implies (x1, X95 Yo) = 1 when Yo<k, € < k which in turn implies 
* 
i(@",k) =- b+ 2-2 ak <o, 


a contradiction. Thus k = a and 


3 


b+ bk - 24 k =b+bp-242p%=0 


a 
3 
with strict inequality holding except for a/b = 4.822 and a/b = 3. Thus (4) is satisfied and 
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i = 
~ (X1,X9,Y9) = 1, Yo <k, & <k. 


This in turn implies that {(o ,k) = 0. (When a/b = 3, k is zero and this region doesn't 
exist anyway. When a/b = 4.822, oe must be one in this region in order to satisfy 
{(v',k) = 0.) 

If a/b > 4.822, (4) <>k = B. If k > B, then ¢ (x,,X9, Yo) = 1 when yp <k, E<k 
which implies : 


* 
i(k) = - b+ 2-2 ak <0, 


w 


x 
a contradiction. Hence k = £, (3) is satisfied, and in order to satisfy f{(Y , k) = 0 


2-b+bk-2k° 





* 
Y (x1, Xo Yo) 4 ’ Y9 <k, € <k. 


bk + ak® 


This is legitimate since Cy is zero for Yo < k, € < k bythe choice of k. 


In summary, 


+= = 4.822 


s 
b 





| 1 if Yo <kK or yo < é 





if Yo>k, Yo > & 


(- 0 if X» > or n <k 





*x 
y Yq» Yq XQ) = if Yo > X%g» Yo > Yi > Yo >k 





otherwise 
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if yy <k and : <s 
* 
Q (x1, X95 Yo) = 


E and yy > k 














3 
* b + bys - $,3 y . 
e 2 2 : 
YW (VysVq.Xo) =< 3 — if Yo > X9, Yo > ¥y> Yo >k 
by» + ays 
L 1 otherwise 
where k satisfies b + bk - 3 + 3 k? = 0 
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PRINCIPLES OF DYNAMIC WEAPON SYSTEMS PROGRAMMING *t 


M. E. Salveson 
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| Over the last several years, there has been considerable interest in 

the problems of effective programming, decision-making, and progress 
evaluation, in the development and production of weapon systems. The 
need for solving these problems was great when budget ceilings were 
less restrictive. At the present time when budgets are more stringent, 
the need is indeed urgent if defense is to be most adequate. The prin- 
ciples and methods described in this article were developed as guides 
to the development and installation of a management system for pro- 
gramming, decision-making, and progress evaluation in weapon systems 
development, procurement, and production in the U. S. Air Force. The 
principles, however, have applicability to programming, decision- 
making, and control in other governmental and in industrial and busi- 
ness organizations. Some of these principles and methods already have 
been applied in other situations; some are new. However, the amount 
that is new is relatively small; we are confident the system is feasible 
and economic within the existing state of the art, and that it would con- 
tribute to efficiency in the development and production of weapon systems. 











INTRODUCTION 

One of the major tasks confronting the United States during the extended cold-war 
period is deciding from day to day and year to year on the continuing best allocation of funds 
for development and production of alternative, competing weapon systems in the Military Serv- 
ices arsenal. This is a difficult decision or cycle of decisions which top military and govern- 
ment officials must make. The decisions involve major issues, including the virtual life or 
death defense of the United States. The decisions are extraordinarily complex — they relate to 
intricate, large-scale weapons with long lead times; they affect economic activity in many sec- 
tors of the nation so that national economic policies also become involved in the decisions. 

The decisions are costly — they involve approximately 18 billion dollars per year, an 
amount approaching the prewar annual federal budget and the Gross National Product in earlier 
years. The decisions affect virtually all spheres of our civilian, commercial, and industrial 
activity, including mining, manufacturing, research, development, and transportation at all 
levels: subsurface, land, sea, air, and interplanetary. These decisions even affect educational 

*Written under University of Michigan Signal Corp Contract DA-36-039 SC-64627 (1959-60). 
Based on researchunder that contract for the U. S. Army under Contract Nonr 233(02) for U. S. 
Navy (1954), and under IBM Corp. Contract for U.S. Air Force (1958-59). The views expressed 
herein are the writer's and do not necessarily represent the views of any of the military serv- 
ices or any of its officers or employees. 
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Aeronautical Systems Center, USAF, for his understanding support and direction; to Mr. Paul 
A. Wilbor, Deputy for Production, ASC, USAF, for guidance and support; to Mr. Edwin A. Boyan, 
now Management Data Systems, Ballistic Systems Division, USAF, for initiating the original 
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79 









80 M, E, SALVESON 


policy at all levels of learning. For example, a certain program decision which led to lags in 
our space probes precipitated a crisis in education and a searching re-examination of our 
national educational objectives, methods, curricula, and resources. 

The cost of the decisions is not only in the weapon systems directly concerned, but also 
in the delay of the systems they indirectly affect and in the cost of recovering any time lost in 
making the decisions. Thus, a delay in a decision often necessitates a large expenditure merely 
to recover the time lost. These decisions also require extensive coordination, ranging from 
the several phases in weapon systems build-up: development, design, tooling, etc., to diverse 
but related subsystems — ranging from ground support, to fuel and fuel manufacture, to per- 
sonnel, and to using the operational system itself. It is to be expected that often it may take 
longer to make the decision than to make the weapon system — a finding of one of the recent 
Senate investigating committees. There are so many necessary commands, functions, and 
levels of coordination that it is sometimes surprising that decisions ever are made, especially 
in face of the speed of changing conditions. 

To add to the difficulties of the military decision makers, there is an increasing rate 
of technological innovation and weapon systems obsolescence. Indeed, deterrency now is 
measured by the adequacy and rapidity of the continuing "flow"' of always new and superior 
weapon systems. The programming methods need to be adapted to speeding and coordinating 
the decision, activities and events in the flow. Many systems often are obsolete before com- 
pleted or tested; existing ones must be phased out and new ones introduced with alarming fre- 
quency. Thus there is no opportunity to ''catch-up"; the military decision-maker becomes 
enmeshed in a complex web of plans and decisions which demand that he act, but which may be 
too difficult or complex to permit him to do so with ease or clarity, considering the informa- 
tion required and available. Large-scale, high-speed electronic computers have been adver- 
tised as a help to him, but thus far they have tended to produce more paper and detail than his 
overtaxed digestive processes can accommodate. 

One measure of the need for convenient, rapid methods of reprogramming can be 
gleaned from the changes or cancellations in weapon system programs. These changes usually 
arise from unexpected technological advances, but each change requires the development of a 
new program to accommodate that change and the reallocation of funds which it permits. Dur- 
ing fiscal years 1957 and 1958, the total dollar value of weapon systems abandoned prior to 
completion for these reasons was in excess of two billion dollars.! The weapon systems 
included, for example, the Navaho, F-108, the Neptune, B-70, and high-energy fuel. Need for 
change and reprogramming also can be due to overruns, retrofits, etc. It seems reasonable 
to attempt to develop a modern, scientifically sound programming system to aid the responsi- 
ble decision-makers, including military executives, Government executives, industry, and 
Congress. It is the purpose of this paper to discuss such a system. 

This report seeks to clarify some aspects of this complex decision process and to pro- 
vide a basis for making the decisions more effectively and rapidly. The particular decisions 
with which it is concerned are those which determine the compromise between total cost of all 
weapons, cost of each weapon, time of delivery of each weapon, rate of expenditure limitations, 
and the continuity of the deterrency capability of the Armed Forces of their "umbrella function" 
vis-a-vis the potential threats to the United States. It also is addressed to the related problems 
of progress measurement in weapon systems program completion. 


lFor details, see Aviation Daily, Jan. 25, 1960, p. 140. 
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NATURE OF PROGRAMMING DECISIONS 

} Time-Cost Relationships 

One major step in aiding the weapon systems programming decisions is to devise suit- 
) able mathematical models of the phenomena about which the decisions are to be made and rep- 
resenting the variables of interest to the decision-makers. These decisions are complex and 

ly involve many conflicting objectives; it would be misleading to imply that they can be repre- 
sented by elementary models. But, to the extent that partial models can be devised, they facili- 
tate these complex decisions, and accelerate analysis of those factors not included in the deci- 
sion model. Thus, a two-variable decision model is developed, which does not purport to 
represent the full complexity of realistic programming decisions, but which is nevertheless 
useful in aiding or evaluating these decisions. 

During cold war conditions, weapon system cost is, of course, a major factor. The 
time of availability is one of the other major factors. Ordinarily, these two variables are not 
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interrelated in procurement or production decisions or negotiation. However, researches over 
the past few years have indicated that, in fact, these two variables can be related. Consider 
the following sequence of examples. 

1. Suppose the only "production" is to transport an object from one location to another 
by a specified vehicle. The time and cost of transportation are related; if the fuel cost alone 


is considered, the cost ordinarily increases as the cube of the speed above minimum cruising 

velocity. Thus, the shorter the time, the higher the cost in cubic relation. 

‘ 2. If other costs are included in the cost of transport, such as crew wages and depre- 

ciation, then there will be both an economy to speed and a diseconomy. Under this formulation, 
a least-cost speed could be found, and the cost would be greater on either side. If there is 

also a value on earliness of transporting the object, then the minimum point will shift. 

3. The writer described earlier” how, in manufacturing, the cost of manufacture varies 
as a function of the time of manufacture. This relationship has been studied further and addi- 
tional examples have been developed. Consider the construction of a building: Our recent 
studies have shown that there is a minimum-cost time for construction of a given building and, 
equivalently, a minimum cost crew which will vary with the kind of building, its structure, 
features, etc. For example, one building contractor who constructs houses in the 1200- to 
1800-square-foot size bracket has found the least cost crew to be four men plus a working 
foreman. This cost increases if the crew is either increased or decreased. In the 2000- to 
3000-square-foot size bracket, another contractor had optimized on six men. 

Thus, if a given house (or other job) takes 1000 man-days with the minimum cost crew, 
it does not follow that 1000 men could complete the job in one day or one man in a thousand 
days. An experienced contractor can estimate quite closely, however, the additional cost 
involved in building the house with fewer or more workmen as well as the increase or decrease 
in the time required to erect the building. If, then, a buyer wishes to complete his building 
early, the additional cost for each length of building period can be estimated. 

Similarly, in manufacturing, a least-cost manufacturing time can be estimated closely 
for a given job and manufacturing facility. If expediting is required in order to save time, 
additional costs will be incurred. These costs include special handling between machine tools, 
S ; using more costly tools and tooling, using higher cost labor to save time, using more 





2M. E. Salveson, "On a Quantitative Method in Production Scheduling,'' Econometrica, 
Oct. 1952. 
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expensive material and processes, etc. In addition, there are indirect costs such as the cost 
of displacing other work. There are also contingent costs which arise. For example, suppose 
minimum time manufacture is sought. Then, one machine tool will be made ready to begin 
work as soon as the job has been completed on the preceding tool; but — since completion time 
is always a random variable with more or less known distributions — there sometimes will be 
idle time while the job is being finished on the preceding tool. This idle time incurs a cost. 

Similarly, in assembly operation, a larger crew can be used in order to accelerate the 
assembly rate. However, experience is the same as in the building construction example. If 
assembly crews are increased above minimum cost size, percentage utilization of their time 
decreases and the costs increase. 

The same relationship holds also in other phases of manufacturing, including in tool 
design and fabrication, in test and equipment evaluation, etc. Again, these have been observed 
in several industrial establishments. , 





Other phases of the weapon system programming also lend themselves to a similar 
analysis and reveal similar patterns in the relationship between time and cost. For example, 
development, design, and prototype test and experimentation follow similar patterns. 

Development work can be foreshortened by many of the same opportunities described 
for manufacturing. That is, it is possible to increase the number of persons in the task force; 
but there usually will be an increase in the percentage required for coordination and in the 
percentage of idle time due to waiting and unavoidable delays. Similarly, higher skills can be 
used, more expensive equipment and etc. Recent studies reported by the DuPont Company 
have confirmed these findings and their general applicability.° 

It would be convenient if in every case it were possible to program a weapon system or 
set of weapon systems for production in that amount of time which has the minimum cost. 
However, it seems a reasonable assumption that any potential enemy will not program his 
weapon systems to be so convenient to the United States. Indeed, the opposite hypothesis 
seems to be more appropriately the case; that he may program in order to maximize our 
inconvenience or disutility. This policy appears to be operative in every phase of the cold war 
in which he has any latitude for maneuver. This being the case, one of our major tasks is to 
meet threats with minimum disutility and in such a manner as to maintain continuity in the 
umbrella function. 

In order to maintain our military posture, each new threat must be countered with a 
counterthreat or countermeasure — and at the time the threat becomes effective. Seldom will 
an unfriendly nation conveniently time his threats so as to coincide with the time of our mini- 
mum cost. Indeed, his strategies may introduce weapon systems which demand countermeas- 
ures in quick time. In this see-saw of threat and countermeasure, programs are changed reg- 
ularly, some are accelerated, others dropped, retrofits introduced, etc., all requiring balancing > 
of timeliness and cost. 


Limitation on Resources and the Coordination Cycle 

Programming studies have shown that each original program or significant change 
must pass through a "'coordination cycle.'' This cycle briefly consists of the following usual 
steps. First, there is the introduction of a need for the program or change, for example, due 


3" Better Plans Come from Study of Anatomy of an Engineering Job,'' Business Week, 
Mar. 21, 1959. 
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to a technological breakthrough either by a potential enemy or by ourselves. This break- 
through obsoletes a weapon system or subsystem and introduces new program requirements. 

In order to establish the new program, the second step is undertaken. This is done by a top 
military executive or committee, such as Weapons Board, USAF, making a first approximation 
to a new program on time and cost of the new weapon required. The new program affects exist- 
ing programs by consuming funds or resources available to them. Then each of the cognizant 
functional and weapon system managers must determine the adjustments to his program neces- 
sitated by the change, including in time and cost. This adjustment is required because there 
usually is a restriction on total resources available, and the sum allocated to the whole pro- 
gram must not exceed that total. Thus, if a new weapon system is introduced and is required 
by a given date, the resources for it can be obtained only from cutback or elimination of other 
weapon systems, either in delivery date or quantity, or both. 

The third step in the coordination cycle is for each of the functional and weapon system 
managers to report to their common superior, whether he be the President of the United States, 
the Chief of Staff of one of the services, or the commander of a large laboratory or agency. 

On this report, the chief common executive and the subordinate managers confer to determine 
whether the first approximate allocation was close enough to that desired. At this time, each 
manager usually pleads his case and usually seeks to advance his weapon system. In popular 
parlance, the "squeaking wheel" principle may operate, with the allocation being made some- 
what proportionately to the relative intensity of the pleadings as well as according to the 
urgency of the threat. This is expected, since the decisions always involve judgment, and 
humans are affected in one way or the other by social pressures. If, on determining the con- 
sequences of the first approximate allocation, they are reasonably acceptable to all managers, 
the chief executive approves the allocation and the decision is affirmed. 

However, as the Senate subcommittee observed, this process of coordination consumes 
so much time that several unexpected consequences might occur. There might be intervening 
changes which obviate the original adjustment and require now a different adjustment. Thus, 
the cycle begins again before it was completed on the preceding change. More or less money 
may be available accordingly, as there is a new administration, a new congress, an off-election 
year, a new recession, a new communist atrocity or threat, or a new Russian satellite in orbit 
about the world or the sun, or on or about the moon, while one of our may be on its launching 
pad. Any of these will change the previous balance or tendency toward balance. Thus, there 
is not only a cycle but a continuing recycle. 


Model of Programming Decisions 

One of the routine continuing activities of these military executives has become that of 
programming weapon systems under these incessant and multitudinous changes. Our purpose 
here is to propose a system and method to aid these executives by increasing the speed and 
power of the analyses for their decision processes. We believe that their decisions could be 
aided by suitable techniques, including special facilities, which quickly would portray the con- 
sequences of any change upon an overall weapon system program. These techniques would 
need to account for the increases or decreases in cost as programs are accelerated, decel- 
erated, or abandoned; for the difference in relative priority of different programs; for the 
limited total amount of resources available; and for the relative contribution to the continuity 
of weapon systems deterency or umbrella function over time. For the purpose of these deci- 
sions, accuracy to the second decimal point (cents) is unnecessary and unrealistic. Hence, the 
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technique can work with approximate data, if only because decisions are made and used now on 
the basis of such approximate data. 

In order to sharpen ideas, the problem can be stated within a mathematical framework 
as shown in the appendix. The reader who does not wish to concern himself with mathematical 
detail may omit the appendix without loss of continuity. 

This appendix develops a mathematical model of the decisions under the assumptions 
that time-cost curves can be generated for the weapon systems and that funds and time are 
available in limited supply. These T-C curves show the relationship of time to develop and 
produce a weapon system (or component) the cost thereof. In general, the cost increases as 
the elapsed time is shortened. . However, if time is increased above the minimum, the cost 
tends to increase in some relative proportion. 

Three programming situations are analyzed. The first assumes that all weapon sys- 
tems have equal priority. The second assumes that there is unequal priority and the third 
assumes that there are arbitrary fixed dates by which some weapon systems are required, or 
there are fixed funds available for them. 

For each case, a computational technique is described which finds the best program for 
the assumptions. In addition, a method is described which selects for any overall program the 
best time and, or, cost to allow for each phase within each weapon systems program, such as 
for development, design, tooling, etc. Other techniques remain to be developed, such as for 
programming under constraints on rate of expenditure. 

One large question underlying this formulation of the weapon systems programming 
problem is whether, indeed, the time-cost curves are available or can be generated. Experi- 
ence cited earlier indicates they can. A reliable method of convincing the skeptic is first to 
get him to estimate the cost for a proposed project. After the time and cost are agreed, ask 
the contractor to accelerate. He will then provide all the reasons for the need to increase cost 
to reduce the delivery time. At least two major weapon systems producers were found which 
are using the T-C concept for their own internal control. DuPont's use is in engineering and 
construction. Indeed, the Department of Defense already uses the basic ideas of the T-C curve, 
but apparently has not thus far incorporated them formally into a programming model. For 
example, on February 1, 1960, President Eisenhower approved an additional $113,000,000 to 
complete the Saturn Rocket program one year earlier. 


Decision Facilities 

In order to provide a usable system for military executives, it is necessary to ''match" 
their characteristics of human perception to the computer which is used. In the initial stages 
of this study, Boyan and Salveson (1958) achieved this ''manager-machine" match with an 
analogue computer feeding solution to print out on an x-y plotter. The analogue is especially 
suitable because its inputs readily can be taken from human sensible devices, such as time or 
dollars-calibrated potentiometers and its outputs are in a form that quickly can actuate human 
sensible media, such as x-y plotter to draw a Gantt Chart. Thus, the executive can make trial 
inputs of his available data and decisions by setting dials calibrated to represent times 
required, money available, relative priorities, etc. An experimental computer system was 
programmed for demonstration purposes at Wright Patterson Air Force Base. It would com- 
pute and print out a Gantt Chart indicating the optimal time of completion and availability of 
each weapon system under consideration and for evaluation of military judgment. This system 
used preliminary data; operational application would have to be based on a system and model 
with actual data. The feasibility was proven, however. 
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No effort was made in this study to devise a method for determining, measuring or 
comparing the military value of different weapon systems. A prime reason for making no 
effort is that it is vastly difficult at the present state of development of value theory. Instead, 
this evaluation is allowed to remain the intuitive or judgmental contribution of the military 
executive, but the exercise of this judgment is facilitated by the model and system. He can 
readily evaluate the worth of any given trial allocation and program and decide to accept, 
reject, or modify it. He can do the latter by a new setting of the human sensible dials. Depend- 
ing upon the computer, the plotter, and the programming, the total time per decision iteration 
is no greater than 30 seconds. This would help reduce total reprogramming time which is now 
12 to 18 months in many programming decisions. One possible further application of the method 
would be in "war gaming" the weapon system programming to evaluate rationally the strategic 
moves which two nations, say blue and green, might make to each other's program moves. 
Thus, the implications of failing to program an ICBM, lunar probe, manned space vehicle, or 
other can be tested, together with program cost and timing. 

Although no effort was made to measure or compare the value or priority of different 
weapon systems, it should be noted that an implied result of a decision, of course, is a com- 
parison and ranking of values or priority of alternatives. In addition, the comparisons are in 
easily communicable form. That is, the marginal cost of the last month saved can serve as a 
decision guide for the innumerable financial and scheduling decisions which must be made by 
subordinate and delegant decision-makers. In evaluating any decision on whether or not to 
authorize a change or expenditure, the delegant decision-maker can determine whether the 
cost per month of saving (or avoiding losing) the expected amount of time is equal to, less 
than, or greater than the marginal cost as determined in the original decision. If the cost is 
equal to or less, then it is a decision consistent with the original; if greater, it is inconsistent 
and would require verification. Thus, decentralized execution becomes consistent with cen- 
tralized planning. 

The system described here was designed to be available to the military executive con- 
tinuously in order to be of best use to him. Indeed, it is designed to permit continuous or ad 
hoc updating of the status of work-in-progress, of relative priority of weapons systems, of 
funds available, committed, flexibility for shifting funds between accounts, new estimates, etc., 
for reprogramming. 

In order to apply such system, it would be necessary to conduct further research on the 
structure and character of the time-cost curve in different weapon systems and phases thereof. 
Such study might have double benefits: It could lead to the data required for the above system, 
and it also might give better understanding of the decision processes and thereby aid in increas- 
ing their efficiency. For example, one general officer reports the Russians are able to develop 
a weapon system in one half the time and cost that we require. Hopefully, our lead time could 
be reduced by accelerating the programming processes. The following describes the concepts 
for such system. 


General Concepts of Programming System 

In the day to day management of a large-scale weapon systems program, there are 
additional problems. It is well known that in maintaining continuity in the umbrella function, 
the future weapons are based successively upon new scientific knowledge and unfolding tech- 
nology. Such advances into the unknown inherently involve uncertainty and risk. 

The uncertainties are many. For example, there is the inability to predict precisely 
that a given proposed system, subsystem, or component can be developed and designed to 
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function in accordance with operating requirements. Thus, a program may involve expendi- 
tures of time and money for a weapon system, all of whose subsystems may operate satisfac- 
torily, except for one or a few subsystems or components. There is an inability to predict the 
time or money (scale and composition of effort) required to complete the development and 
design of a given weapon system and its subsystems and components. Thus, a given project 
may require more or less time and money than initially programmed. 

Progress in science and technology frequently is discrete and often characterized by 
"breakthroughs."' This is the case both for progress in America and its allies and in the com- 
munist nations. A consequence is that in a relatively short period of time, often less than the 
programming cycle of a major weapon system, new scientific knowledge or technology may 
make obsolete a programmed weapon system. The breakthroughs may be either by friendly or 
by opposing scientists or engineers. The consequence is the same: existing programs must 
be altered or abandoned and new ones developed to meet the new opportunity or threat. 

There are other, less frequent, courses of program perturbation. Although not prac- 
ticed by reputable weapon systems producers, some contractors deliberately or inadvertently 
through lack of experience may understate time or cost estimates. The consequence is the 
same however; whenever actual values are significantly different from estimates, they perturb 
the program and necessitate adjustments to the program and to time and cost estimates. 

Each program adjustment is subject to the earlier described overall financial and time 
limitations and must be processed in the reprogramming cycle. In addition, whenever a pro- 
gram for a weapon system has commenced according to a specified time, cost, and performance 
plan, it involves inflexible outlays for facilities, procurement, etc. Any acceleration or decel- 
eration will result in increases in cost which affect the overall program. The increases are 
due (1) to the partial loss of the program outlays already made and (2) the need to replace those 
program outlays in developing the new program, and (3) relatively constant or invariant over- 
head. These new outlays usually require "undoing" or discarding some of the previous arrange- 
ments and providing new. Thus, the effect of changing programs, whether to accelerate or 
decelerate, typically increases weapon systems costs. 

Because the over-all programs are temporally and financially constrained, a change in 
time or cost of one element or weapon system affects the money available to the others. The 
resulting dynamic weapon systems program control problem is how to reoptimize a whole pro- 
gram continuously as perturbations occur in time, cost, science, technology, or etc., in any 
program phase of any individual weapon system therein. 

System Requirements - To handle this reprogramming problem requires an effective 
management and progress reporting system designed to meet the unique meeds of weapon sys- 
tem programming. This system should satisfy certain requirements. 

1. It should be based on the structure of the specific activities and milestones in the 
program. 

2. It should obtain and correlate information on the intensity of effort and cost of the 
activities so timed to accomplish the full sequence of milestones in the program. 

3. It should provide for the flow of information on progress and status toward achieving 
the milestones, together with revised time-cost estimates for any adjusted activities. 

4. It should include a rapid procedure for reprogramming under each program change 


and for communicating information on each reprogram to the managers responsible for 
affected phases. 
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5. It should be adapted to the over-all system within which the programming is done 
and to the decision needs of program managers and engineers. This over-all system embraces 
Congress, Bureau of the Budget, Defense Department, Military Services, and the contractors 
and subcontractors. 

The nature of each of these requirements and suggested methods for satisfying them 
are described hereinafter. A greatly simplified block diagram of the proposed system is 
presented in Figure 1. In the demonstration model, a Reac 302 and 201 with dollar and time 
calibrated dials were used, together with a Mosely x-y plotter for printing the Gantt Chart. 
The T-C curves used were taken from Lockheed data. 

1. The structure of the activities and the milestones are basic to program formulation 
and control. An inherent characteristic of weapon systems programs is "structural inter- 
dependence" among the many tasks to be performed. For a simplified example of structure, 
when a Service develops and produces a desired new weapon, the general operating require- 
ments must be specified before design competition can be undertaken. Thus, the sentence 
"G. O. R. specifications must precede design competition" is a statement of the structural 
relationship between the two activities. The two milestones resulting from these activities 
are (1) the statement of specification for G. O. R. and (2) the preliminary design proposals 
from competing contractors. In smaller detail, the design of the main wing spar must precede 
the design and manufacture of tooling for producing the spar. Thus, (1) wing spar design must 
precede (2) wing spar tooling.* 

Another class of relationships also is important in weapon systems programming. The 
design of a complex system or subsystem often requires simultaneous interactive design of the 
subsystems and components. For example, the design of a complex bomber requires simul- 
taneous consideration of the performance, weight, and space characteristics of its subsystems 
and how these affect each other and, in turn, how they integrate with each other to yield the 
over-all performance of the aircraft. For example, of the circular interaction, the desired 
bomb load affects fuel capacity and range; guidance and navigation equipment affects bomb and 
fuel load; range and speed affect the amount, kind, and weight of navigation equipment. Thus, 
these must be considered simultaneously in the design process. In large degree it is the com- 
plexity and intricacy of the simultaneous design of modern weapon systems that have created 
the need for the weapon system concept and weapon systems contractor method of procurement 
and production. 

The weapon system contractor thus has responsibility for integrating design character- 
istics of all subsystems and components and for procuring and assembling these in the com- 
pleted weapon. Yet, the Military Services can not and do not release themselves from the 
important responsibilities of evaluation and control of the evolving weapon system design. 
Prudent weapon system management requires that the Services' representatives con‘: :bute to 
and participate in deliberations and decisions on design compromise among the weapon system 
components and subsystems. Final performance characteristics and satisfaction of G. O. R. 
are responsibilities of the Services which cannot be delegated. Thus, the weapon system con- 
tractor needs to determine and elucidate the effect of compromises among subsystem design 
and performance and as these affect over-all performance and satisfaction of G. O. R. Then, 


4For the graphic and algebraic method of handling these relationships, see M. E. Salveson, 
"The Assembly Line Balancing Problem,'' Journal of Industrial Engineering, 1955, or see the 
Appendix to this article, 'Principles of Dynamic Weapon System Programming." 
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Figure 1 - The proposed system 
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the weapon system project officers have the information for rendering decisions which inte- 
grate subsystem design, system design, and system performance. 

For purposes of programming, ''simultaneous design" of subsystems in the design of 
the weapon system means that the activities and the milestones "must be concurrent" (rather 
than "must precede," as in the earlier above examples). Simultaneous activities usually begin 
with a common milestone and lead simultaneously to one or several concurrent milestones. 

For example, in the above activities, they might begin with preliminary estimates of range, 
speed, and allowable equipment load, and terminate in specification of each equipment's weight, 
operating characteristics, and space requirements in accordance with current state of the art. 
The simultaneous activities can be represented graphically and algebraically by the same tech- 
nique for representing sequential activities. With structure stated in explicit form such as this, 
the contributors to the program and design decisions have a better basis both for program for- 
mulation and control and for coordination of design decisions. These data are one of the infor- 
mation feedbacks from contractor to military in the proposed system. 

2. All activities required to complete the weapon system can be conveniently and 
accurately portrayed as above for program planning and control purposes. As will be seen, 
these data are essential to the several succeeding steps in the program. In particular, they 
help determine both (a) the time-cost curves for achieving the milestones set forth in the pro- 
gram under different time spans and (b) the cost of allowing departures from the structural 
relationships in order to reduce the elapsed time of a program. These can be used to develop 
the time-cost feedback information in Figure 1. 

There are many ways in which the time can he reduced or extended in a given program. 
For one example, in manufacturing one contractor used multiple tooling in order to reduce 
the time required to fabricate or assemble a specified number of parts or assemblies. The 
additional tooling tends to increase the cost of the program, ceteris paribus. In design, the 
elapsed time may be reduced by overtime work, by using key or higher paid personnel to speed 
work otherwise delegated to lower paid personnel, etc. 

A second method of compressing the elapsed program time is to take a calculated risk 
to reduce the time required to complete a sequence of activities. For example, it may be 
decided to begin tool design and fabrication some specified number of day ; prior to completion 
of spar design. Thus, tool design would be based only upon preliminary engineering designs. 
The cost of using this method of acceleration is in the expected cost of having to discard tool 
designs because of probable changes in engineering design. This expected cost is the product 
of the probability that there will be such an engineering change by the cost of any tool designs 
which must be discarded. In general, the probability of engineering change increases as the 
overlap increases. Thus, if final engineering design release is awaited, the probability of sub- 
sequent change is low; if tool design is begun early in the engineering design, the probability 
of change is high. Limited evidence indicates that the cost of change also tends to increase as 
the overlap increases. This is due, apparently (a) to the increased investment which results 
from a longer period of time in which to carry out the follow-on activity before late changes 
can occur. This increases the investment in the follow-on work. Also, it is due to (b) the 
increased magnitude of change which can be made in the output of the preceding activity dur- 
ing its earlier periods as contrasted with those during its later periods. This is only a tend- 
ency, however, and is reversed if the period is long enough so that there are significant changes 
in technology or state of the art. 
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Engineering design superficially appears to have a time-cost curve in which cost 
monotonically increases as a function of increasing time allowed for design. This relation 
often has been observed, but it actually involves the following phenomena. If engineers are 
allowed more time, they design a different weapon system from that if they are allowed less. 
The difference will manifest itself in different operating characteristics, cost of manufacture, 
reliability, etc. This relationship can be handled easily in the model described here. For 
example, the program manager would seek that length of time for engineering and manufactur- 
ing, such that marginal cost of increased engineering is equal to the marginal cost of manu- 
facturing, etc., as a function of amount engineering. 

As a practical aid to programming, the utility of computing expected cost as a function 
of amount of overlap for individual pairs of activities is doubtful. For any individual pair of 
activities, a change in one which affects the other is usually a one occurrence event (or small 
number of occurrences). However, the law of large numbers can apply quite usefully to these 
more or less randomly originating changes, if applied to a collection or large sequence of 
activities. It is believed that computing the expected cost of overlap for single pairs does not 
give information of sufficient value to use in programming such pairs. An exception might be 
especially high cost items. Otherwise aggregative behavior appears sufficient for most pro- 
gramming purposes. In view of the universal problem of engineering change control in com- 
plex systems, much data are latently available that would be useful to this end. 

Another datum for measuring the status and progress is the cumulative man-hours or 
dollars expenditure ''S'' curve. It has been found in many laboratories and companies that in 
well programmed and managed weapon systems, the characteristic cumulative curves fall 
within a narrow range of values. Thus, this datum for any program can be used to evaluate 
whether or not a program is in control and progressing as planned. But, that a program falls 
within the range does not necessarily assure its proper progress. Alone, the cumulative curve 
cannot detect a failure to achieve programmed progress; it must be related to progress as 
measured also by achievement of programmed milestones, and time and cost estimates for 
future milestones. 

3. Uncertainties are always involved in carrying out the activities to achieve any mile- 
stone. Hence, program progress reporting is essential to control of program time and cost. 
As a program or program phase evolves, it is necessary to determine the status of the pro- 
gram, compare it with planned status, evaluate departures, and reprogram according to cur- 
rent requirements, including time, cost, and funds availability. For this purpose, it is neces- 
sary to report also the time-cost curve for the program as it then exists, considering the work 
completed, work to be completed, and the changes, if any, in the remainder of the program. 

Our experience in industrial studies indicates that estimates of probabilities of com- 
pletion of specific activities provide little useful information. Instead, it indicates that it is 
better to make experience-based estimates of time and cost (including expected cost from 
overlaps) and to use these for reprogramming. There is a real self-discipline in having the 
person responsible for the work participate in making such estimates and then be responsible 
for achieving them. This self-discipline is vitiated if probability distributions of range of 
expected completions are used. Psychologically, the person responsible has a keener sense of 
responsibility to his performance, if he is considered to have control over his unique work 
assignment rather than if his work output is considered to be a random variable defined over 
some probability distribution. On the other hand, for economical and realistic programming, 
the programmer of a complex of activities should include the stochastic element in program 
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formulation, but then he should treat each person responsible for an activity as if the outcome 
were under that person's control and determination. Many of these concepts were developed 
by E. A. Boyan in 1944 while at the Radiation Laboratory, M.I.T.° One facet of the over-all 
system described here has been nicely developed independently by Malcom et al. © It is PERT 
(Performance Evaluation and Review Technique). PERT corresponds in part to the feedback 
portion of our system. Indeed, PERT is a computerized version of E. A. Boyan's Target Com- 
mitment Scheduling. PERT utilizes an algebra of ordered relations similar to Salveson's, ’ 
which is presented in the appendix here. Only in regard to PERT's probability statements for 
expected task completions and scheduling and load commitments do we feel PERT could be 
improved for progress reporting. It is understood that PERT increased performance efficiency, 
as also did Target Commitment Scheduling. It is conjectured that the searching analysis of 
program performance per se required by both techniques induced much of the improvement. 
This analysis is the same in PERT and TCS. It is believed that the ''commitment" feature of 
TCS will have a stronger tendency to induce achievement to schedule than the "probability" 
feature of PERT. In any case, the favorable experience of both TCS and of PERT lends 
credence to the more comprehensive system proposed here. 

Status reporting for program control requires a compact system of notation for report- 
ing milestones and activities. Indeed, the magnitude of the military development, design, and 
production programs, makes this an essential requirement if the information volume is to be 
kept within manageable proportions. (The matrix notation for structuring activities within a 
program provides this notation, but it also conveniently relates program status to program 
structure for reprogramming. Our experience indicates that reports of percent achievement 
of milestones are almost meaningless. Instead, they should state the milestones achieved and 
the time, cost, and program to achieve the remaining milestones.) 

In view of the old aphorism, more true with modern weapon systems, "that for want of 
Se ee. , the kingdom was lost," the system should be capable of identifying and report- 
ing departures of even small items which can adversely affect the broader program. A com- 
pact notation requires, of course, some regularity in methods of program description and in 
the encoding and decoding procedures for transmitting information on the program. This is 
accomplished if standard program structures are used throughout the over-all program. At 
the same time, the frequency and detail of reporting should be consistent with the self-discipline 
principle. It is desirable that departures which do not affect a program be unreported, both to 
reinforce the self-discipline of seeking to return to schedule and to avoid such mass of infor- 
mation at higher echelons that it becomes unintelligible. 

4. It is important that contractors maintain their autonomy and management preroga- 
tives. The flow of information from the programming centers back to contractors thus must 
be designed to assure an optimum balance as to detail of corrective instruction, volume of flow 
of information, and sufficient delegation of authority to allow decentralized decision-makers to 
be effective. One possible method of achieving this goal would be to communicate the marginal 
costs of saving an additional day (or month) as computed by the reprogramming decision. That 
is, under a changed circumstance, a weapon system may be required earlier and hence it would 


, See E. A. Bryan, Lecture Notes, Course 15:71, M.I.T., 1946, Boyan terms these concepts 
"Target Commitment Scheduling." 

D. G. Malcolm, et al, ''Performance Evaluation Review Technique,'' presented at meeting 
of Operations Research Society of America, Washington, D. C., 1959. 
Op. cit., 1955. | 
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be allowed a higher cost in the programming computation. By indicating such numbers, as 
well as dates, to military representatives and contractors, they would be enabled autonomously 
to develop the new work plan and program for that greater or lesser priority. Consequently, 
by providing them decision criteria rather than decisions, the contractors can maintain their 
autonomy and yet make decisions consistent with and in support of the over-all program. 

The frequency of reporting and response also must be optimized so that there is neither 
loss of control nor too frequent program correction. There is presently opportunity for unde- 
tected over runs. Yet, the remedy should allow for the self-correcting tendencies of contractors. 

The compact matrix notation permits the convenient communication of the outgoing 
instructions which again are in accordance with the structure of the program and the assign- 
ment of responsibilities therein. 

5. It has been recommended by some authorities that government, military and indus- 
trial executives should sit down around the conference table as a method of improving the 
decision making processes in weapon systems planning and programming. While direct person- 
to-person communications is always desirable, this method does not appear to be practical as 
a correction of the basic programming problems. In the government there are dozens of offices 
or committees which would need to be represented at the table. In the military, there would be 
many more. There would be hundreds of contractors or potential contractors who would have 
the right to be present. In addition, each person at the table would have to be backed up by his 
staff specialists; legal, financial, estimating, contracts, engineering, etc. Of course, each staff 
specialist probably would require his files and his administrative assistants. Thus, the table 
would not be a table, but an auditorium-office building. 

The practical requirement is that the programming system must operate within the 
offices and locations of established government agencies, military commands and contractor 
corporations. Each must continue to make those inputs which are indeed required to any pro- 
gramming decision. Of course, the decision procedures which are used, the methods and con- 
tent of communication, and the facilities for these may be changed. Also, improved bases of 
coordination may be achieved. For example, in a major procuring command, such as the AMC 
Aeronautical System Center, a comprehensive integrated programming system might facilitate 
providing the information for more immediate coordinated reprogramming of all weapon sys- 
tems and subsystems under the impact of a major fund or other change. This could coordinate 
not only the schedule and financial requirements on the weapon system directly affected, but 
also on the weapon systems indirectly affected. In the large complex programs which such 
centers administer, the coordination and paper work load is high: the first ADPM installation 
in one such center is to be used for "document control," i.e., to control the "paper work," the 
flow and processing of documents involved in this process. 

Additional Applications - If it can be developed economically and reliably, there are 
several additional uses which could be made of a model of a complex program. These would 
include (a) estimating optimal "slack" for stability and control of the propagation of disruption 
introduced by a perturbation to a program step; (b) generating time cost curves for programs; 
(c) determining the effect of selected priority schemes upon the flow, production, and delivery 
of weapon systems; (d) determining effect of selected amount of various resources upon the 


flow, production, and delivery of weapon systems; and (e) evaluating adequacy of the program- 
ming system etc. 


The internal operation of a program can be represented as a multistage, multichannel, 
parallel-processing queuing problem with arbitrary interdependent queue disciplines, with 
















er 


ors. 


ff 





2 

4 
& 

a 
3 
# 
re 





ee 





PROGRAMMING DYNAMIC WEAPON SYSTEMS 


patterned but variable processing, and with transient loading. It also has sequenced birth 
processes. The characteristics which lead to this formulation are as follows. 

Usually the weapon system program, whether in industry or military, is handled by a 
functional organization which services several programs. Each functional department may be 
considered a servicing facility in the queuing theory sense. Usually there is a sequenced flow 
of items to be processed through these departments. Each department thus may be considered 
one service stage in the sequence. To service completely any arriving item usually requires 
several successive stages. The service times in each stage may be assumed constant (although 
this can be relaxed), while the holding or delay time at any stage may vary with the state of the 
system, i.e., the load. 

In many service centers there may be several persons, equipments, or other facility 
for servicing. Thus, many parts of the system may have multichannel capacities. The incom- 
ing items can be distributed over these chaunels in any arbitrary manner, e.g., to balance load 
on the several channels, to assign items in particular categories to selected channels except 
when they are overloaded, etc. 

Often an item may require parallel processing of interdependent subitems. For exam- 
ple, to ship an item may require preparation of packing slip, invoice, bill of lading, and factory 
material release. These may be processed in parallel by different channels in the organization, 
either simultaneously or independently. 

The queue disciplines may be arbitrarily determined and represent priority schemes 
which may be adopted to expedite one weapon system or another. These disciplines affect the 
holding and over-all flow times of the variously prioritied items in the total program. There 
is interdependence among the servicing. That is, for example, there may be a limited amount 
of resource, Say money or material which is assigned to items as they are serviced. Depending 
upon the assignment policy (for instance, first come first serve), the servicing time of later 
items may be influenced by earlier items. Thus, some may be held until the resource is 
repienished, or not serviced at all. 

There usually is a pattern in the processing of each class of items. One class may 
first be processed by the Deputy Chief of Staff for operations, by DCS/Development and, finally 
by, DCS/Controller before released to an external organization. Other classes may have dif- 
ferent patterns over these same offices. In some cases the patterns may be variable, depend- 
ing on decisions at a service facility or upon the loaa at the facility. For example, if a negative 
decision or service is rendered by DCS/O, the item may never reach DCS/Controller; or, if the 
screw machines in a factory are overloaded, a job might be assigned to a subcontractor or to 
turret lathes for alternate servicing. 

The servicing loads in a program are transient in many phases. During one phase, 
there may be an increasing load in engineering followed by a decline to a steady state for the 
balance of the program. This also applies to tooling, subcontracting, and manufacturing as 
well as to operations in the military portions of a program, such as in source selection, price 
negotiation, and contract administration. The loading and holding times during these transient 
periods are of critical interest and can be studied via the model. 

The sequenced birth processes arise from the fact that an item being serviced at one 
stage leads to origination of several items for processing at other stages. For example, a 
development directive may give rise to several follow-on actions in offices of AMC and ARDC 
where each is concerned with one phase of implementation. 














M. E, SALVESON 


Insofar as any one servicing facility is concerned, its incoming individual items of input 
are usually randomly timed, but with more or less predictable mean rates during short periods 
of time. 

Naturally, the queuing problem which the above characteristics generate, is complex. 
However, similar organizational problems have been handled by an equivalent formulation, and 
there is no reason to expect that this could not also be handled. The "solution" of the problem 
surely must be by simulation to determine optimum values for such characteristics as "slack" 
(excess capacity), priority schemes, and resources, as well as for generation of operational 
characteristics such as families of time-cost curves, delivery schedules, adequacy of the 
programming system, and efc. 


Weapon System Organization 

Characteristically, both military and industrial organizations are organized by function 
(finance, production, contract, etc.) and by weapon system, or product. Within any of these 
larger organizations, program decisions usually must be coordinated over all these functions 
and weapon systems. It would be an important contribution to the decision process if a decision 
aiding facility, such as described earlier, were available to quickly compute and portray the 
organization-wide implications of a programming decision in any part of the organization. In 
the absence of such system and facility, one method of trying to solve difficult decision prob- 
lems is to alter the "form" of organization. 

One form may attempt to solve the complexities of multiple weapon system programs 
(or product programs in industry), by organizing weapon system project officers (or product 
managers). Each Weapons System Organization (WSPO) then acts as the expediting and coor- 
dinating corps for its assigned weapon system. Functional responsibility may or may not be 
assigned to the WSPOs, depending upon organizational philosophy of the military service (or 
industrial enterprise). If it is, then there is an advisory functional staff to control quality of 
functional decision by WSPO or Product Managers. If not, the functional departments make the 
decisions subject to coordination and expediting by the WSPOs. 

In either case, the "problems" clearly are not eliminated by the form of organization. 

If the primary mode of organization is by function, the functional decisions are usually well 
made; but weapon system coordination decisions are usually poorly made and progress is 
slowed. If organized by weapon system, each weapon system program is usually well coor- 
dinated internally, but interweapon system program coordination and functional decisions tend 
to suffer. There always remains the problem of testing the impact of a programming decision 
simultaneously on the several weapon systems and functions. 

We have found that proper design of decision processes and communication facilities 
are as essential to full program coordination as is the organization form. Indeed, by selection 
of these processes and facilities, the programming efficiency of the organization becomes much 
less dependent on the form of organization. To quote one anonymous military executive, 'When 
we are unable to handle the problem, they reorganize our inability." Now, of course, clear 
lines of authority and channels of communication are essential to effective organization. But, 
organizing one way or the other, does not change the difficulty of the decisions. A complex 
decision is not made simple by subdividing it into many small partial decisions. This only 
increases the amount of coordination required, and thus makes it more difficult. The resolution 
is in increasing the power of the decision processes, facilities, and techniques so as to cope 
with the complex decision. These comments on the organizational problems are most 
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sympathetic of the difficulty of the tasks facing military organization planners and executives; 
they are not intended to be critical, but only illustrative. The U. S. Military Services comprise 
probably the largest operating organization in the world, different from and more complex than 
commercial businesses. They are highly technical and highly dispersed; they are carefully 
scrutinized by many other groups: congressional committees, civic bodies, governmental 
agencies, and etc. Our intent is to illustrate the difficulty and complexity of their decision 
problems so that adequate understanding and facilities can be accorded them. A particular 
difference is that a large industrial corporation can be decentralized into departments each with 
its own budget, program, product line, and markets. Decisions made in one department have 
little effect upon the others. This is not the case in the Defense Department. Programs and 
budgets in one service, fund series, or project affect and are affected by those in other serv- 
ices, series, or projects. Much wider and more complex coordination is required in the mili- 
tary than in industrial organizations. This increases the difficulty of programming decisions, 
and minimizes the value of product or weapon-center organization structures. Facilities for 
coordinating programs and for making decisions remains the essential requirement of effective 
programming. 


CONCLUSIONS 

Most of this analysis has been in terms of time of delivery and cost. There are, of 
course, other important program criteria: weapon system performance, quantity, reliability. 
The two-dimensional, time-cost curves can be generalized to include these other variables, 
but at some cost in computational efficiency. The result, however, is to establish the basis for 
more rational trade off rates between these weapon system variables. 

One problem of including "performance" in the decision variables is to quantise this 
attribute. Arbitrary measures of performance can be devised combining such performance 
variables as speed, range, ceiling, armament, etc., but they would require considerable test 
and validation by experienced military officers before they could be used with the same degree 
of confidence that could be attached to the other variables. Another improvement would be to 
include the expected annual and total cost of each weapon system for its full life, from develop- 
ment to discard from the operational inventory. 

The principles of weapon system programming described in this article imply a pro- 
gramming system based on system concepts and methods derived from analysis of the present 
programming practices. While our experience in developing other management systems 
assures us that these principles will lead toward more effective and economical programming, 
it also suggests that such system should be developed on an evolutionary basis and should be 
implemented through parallel operation, test, and refinement. There remains much analysis 
of operations, model building, and testing; and of facilities design which must precede develop- 
ment and application of these principles. 

The least that can be said at this time about developing any such programming system 
is that our government and military executives urgently need it. Surely it should be based on 
scientific methods of systems design and on facilities and processes which are as modern and 
effective as those used in the design and construction of the weapon systems which it programs. 
In addition, it should have large carryover value, as is often the case in the research and devel- 
opment studies for weapon systems. Other branches of the Federal Government have similar 
programming problems, as well as the larger state governments and industrial organizations. 
Thus, the development and implementation of such system could pioneer other systems and 
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lead to greater efficiency in Government and industry. We are happy to report that a system 

of these characteristics is being developed in the U. S. Air Force for more effective weapons 
systems programming and control. It is known as WSPACS—Weapon System Programming and 
Control System. (Copies of the Project Guide to the Air Force ''Weapon Systems Programming 
Control System (WSPACS)" are available to persons in the Defense Agencies and Industries. 
This project Guide describes program development and the installation of a Programming and 
Control System. Write to: Dr. D. E. Danielson, Headquarters, Air Material Command, Wright- 
Patterson Air Force Base, Ohio. Attention: MCPA.) 


APPENDIX 


Three weapon system programming cases are considered herein — (a) equipriority, 
(b) unequal priority, and (c) fixed dates, or costs, among the weapon systems. The following 
are the concepts of the decision processes and the notation used. 

Let T be time and T,, i=1,..., n be the time at which the i” weapon system is to be 
available. tj is the time required in the j 


phase of the jth system. = tj = T; assuming no 


J 
overlap. Let M be the total amount of money for a program, M; the amount for the jth weapon 
system and mij the amount for the jo phase of the jth system. 


M =).M, =). m,. 
aM, Ta 


Inasmuch as the weapon system race seldom, if ever, allows approaching minimum cost, the 
problem is simplified, without loss of realism, by assuming that all time and cost values lie in 
the left-hand side of the concave (upwards) time-cost curve. In an earlier study, an exponential 
function was used. The realistic problem involves curve fitting to empirical data. In this 
analysis, this left side of the T-C curve is approximated by a hyperbola of the form: 


dij ot rr 
ey" ta, * Sy 


ij ij -,N 


where dij» ej » and hij are constants om to define + peace closely approximating the 
left side of the time cost curve for the j~ phase of the i weapons system. By inspection, 
ej is the minimum cost regardless of the time allowed. dij and hij are parameters of mini- 
mum time and of slope. 

In the case of one weapon system, with several phases, the problem can be stated in 
one of the following ways. (A) If M; is fixed, what are the best T; and ti and resulting m,; 
for that i and for j=1,..., m. That is, if M; units of money are allocated to weapon system 
i, what is the time to allow for its completion and the amount of time for completion of each 
phase. Of course, given either the time-cost curve or function, the time, . can be read 
directly. The best time for each phase is not available by direct inspection but must be com- 
puted. One way is to solve for that set of points at which 
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This solution begets that allocation which will allow to each phase an amount of money which 
will give the same marginal cost of last unit of time saved for all phases. This is, of course, 


id the minimum cost allocation of time and money to phases. (B) If T; is fixed, the amount of 
ig time for each phase for least cost program can be computed by solving the equations: 
| dm, ; 
. =r, j=l,...,n 
it ; dt; ; 
subject to 
: >» SS = 
q 7 i 
y r is the slope at which the optimum solution is found and is not known in advance. 
In other more general situations, there are several weapon systems, each with several 
| phases. The simplest case is equal priority among them. Even in this case, however, the 
BS meaning of "equal priority" is not clear, unless defined in some quantitative way. One such 
e : definition is available in the solution to the several weapon system case. If, intuitively, each 


10 a weapon system has equal priority, then this can be interpreted to mean the marginal cost of 
the last day saved for each is equal for all. Thus, in the equal priority case, if a given amount 


on Z 
a of money M, is available, then the "best" allocation of that money to weapon systems is 
: obtained by solving 
dM; 
—-=r; i=l, »m 
i aT. 
i 
. and 
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This case would represent the practical situation wherein it was known or suspected that the 
opponent is developing advanced weapons systems which ours are designed to counteract or 
neutralize, and where each has equal destructive potential but where we have no knowledge of 
the completion dates relative priority assigned by the enemy to his weapon systems. 

Of course, once the equipriority solution is obtained, it is necessary to solve within 
each weapon system for the minimum-cost allocation of time and money to each phase. 

At this point, it is emphasized that the above definition of equipriority is arbitrary but 
meaningful. That is, if two projects have the same priority, this does not mean intuitively that 
the same amount of money should be applied to each. Rather, it seems more sensible to say 





that the same amount of money would be spent in order to save one more day in each. This is 
the meaning of the above solution. 

If it is known that the enemy is not giving equal priority to his alternative weapon sys- 
tems, it is possible to allow for a strategic weighting of his inequality. It may be desirable to 
emphasize development and build-up of those systems to which he is giving the least priority 
or vice versa. But, whatever strategic decision is made on relative priorities, it can be 
reflected in the relative marginal cost of the last day saved in each weapon system. In this 
case, a weighting factor is included in the system of equations 
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where a; is the index of relative priority of the jh weapon system. 


Another alternative is the case in which the enemy's intentions are known for at least 
one or more weapon systems. If the threat is to be countered with no discontinuities in the 
umbrella function, then the allocation of funds need provide that sum to the countering weapons 
which will produce them by the time desired. In this case, the weapons’ costs for those whose 
dates are fixed can be determined by direct inspection of the time-cost curve. Then, if total 
funds are fixed, the amount remaining for other weapon systems can be determined by deduct- 
ing this cost from the total. The allocation can be made with equal or variable priority among 
these remaining weapon systems. 

The following describes a method for analyzing and controlling the program progress 
for the build-up and delivery of a weapon system. A first illustration is given based upon an 
example with which all men (military and civilian) will be familiar; namely, dressing. It then 
is extended to show applicability to another simplified structure which is more typical of 
weapon system program activities and milestones. 


ANALYZING THE PROGRAM 
The first step is to make a systematic analysis of the program and all its required 
activities and milestones. The programmer enumerates all the program steps, decisions, 
requirements, and estimated time and cost for each. It is not necessary at this step to arrange 
them in any particular sequence, although it is helpful to try to follow a beginning-to-end path. 
In the case of dressing, milestones and activities are presented in Table 1. 











TABLE 1 
Milestones and Activities 
Milestone Activity Activity Time 
0 Awake -- -- 
1 Showered Taking shower 2 
2 Dried Drying 5 
3 Socks on Donning socks 4 
4 Shoes on Donning shoes 8 
5 Undershirt on Donning undershirt 2 
6 Shirt on Donning shirt 1 
7 Tie on Donning tie 6 
8 Coat on Donning coat 2 
9 Undershorts on Donning undershorts 3 
10 Trousers on Donning trousers 4 
11 Dressed -- -- 
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DETERMINING PRECEDENCE RELATIONS 
It is necessary to determine the "precedence relations" for the milestones and activi- 
ties in the program. ''Precedence relation" is a name we use for the relation defined by the 
statement "must precede."' For example, when you dress (assemble clothes onto yourself) you 
must (or at least it is conventional), put on your socks before you put on your shoes. Hence, 
we say "socks must precede shoes," in the sequence in which we assemble clothes onto our- 
selves. For the purpose of explaining this concept and the subsequent illustrative problem, let 














s be 
a us consider the dressing problem more fully here. It is very similar to any typical program 
3 i problem and illustrates the concepts in terms familiar to everyone. 
3 Instead of writing "socks must precede shoes" let us now develop a method for writing 
ES this notation more compactly. Use the notation: 
a 
4 To express this relation read from left to right, ''socks must precede shoes."' By analysis of 
a the precedence relation between our other elements of clothing we could write for all other 
dressing operations: 
(Undershirts Shirt Tie 
C Undershorts Trousers ) 
e 
; ) If we add the operations "shower, dry, and dressed,'' we could write these altogether as 


presented in Figure 2. We call this a "precedence graph."' It summarizes the precedence 
relations. It states, for example, that we must put on our undershorts before our trousers, 

our socks before our shoes, etc. It does not say anything, however, about socks in relation to 
undershirt, socks in relation to trousers, shoes in relation to undershorts, or shoes in relation 
to trousers. We assume, therefore, that we may put on these "unrelated" clothes in any 
sequence, only so long as in the over-all sequence we always select, (a) socks to precede 
shoes, (b) undershorts to precede trousers, etc. Hence, we could use any assembly sequence, 
involving socks, shoes, undershorts, and trousers, from among the following possible alterna- 
tive sequences (reading from left to right). 


1. C Socks )}-{ Shoes Undershorts Trousers ) 












































2. ( Socks Undershorts Shoes Trousers ) 
3. (_Undershorts Socks }+{ Trousers Shoes ) 
etc. 


From the preceding sequences, we derive a further principle about the precedence relations. 
That is, for example, "socks must precede shoes, but (certain) other tasks can intervene." In 
the above possible sequences, we interspersed "undershorts and trousers" with "socks and 
shoes," with full freedom as to sequence so long as socks preceded shoes and undershorts 
preceded trousers. In situations such as this where one item need neither follow nor precede 
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another, we say that the first is "unordered" with respect to the other, and vice versa. Items 
which are unordered with respect to each other are written (as in Figure 2), by showing them 
on different parallel lines. 








i 4 
[socks} [SHOES -— 
fe) | 2 
5 6 7 8 W 
[ awake }—[ srower }—fory}+{unoersniat}—sniat}—{rie}—{coar}+-{oresseo| 
[ UNDeRSHoRTS }—{TROUseRs}#—_—<1 


Figure 2 - Precedence graph 





























Another method of recording the precedence relations which some engineers and pro- 
grammers have found more convenient than the precedence graph, is the "precedence matrix." 
It is constructed as follows: Make a square table with as many rows and as many columns as 
there are milestones or activities in the program. Number the rows from top to bottom and 
the columns from left to right. Let each number correspond to the number assigned to one 
milestone, as for example, the numbers on Figure 2. To record the precedence information 
onto the matrix, read into the matrix from the left side along a selected row and enter where 
that row intersects with any column: 

(a) A minus one whenever the item represented by the number of the row "must fol- 

low'' the item represented by the number of the column. 

(b) A plus one whenever the row item "must precede" the column item. 

(c) A zero whenever the two items are unordered (including when the row and column 

are for the same item). 
To illustrate, let us make a precedence matrix for dressing. Number the milestones, as shown 
in Figure 2, and as follows: 


1. Shower 5. Undershirt 9. Undershorts 
2. Dry 6. Shirt 10. Trousers 
3. Socks 7. Tie 11. Dressed 
4, 


Shoes 8. Coat 


The precedence matrix then is made by listing these numbers along the left and along the top, 
and entering a zero, a one, or a minus one as appropriate, in each position in the matrix. We 
always assume an item is unordered with respect to itself. The precedence matrix for the 
dressing problem appears in Table 2. The matrix is read as follows: take, for example, Row 
4, which represents shoes. The entry -1 under Column 1 says that "shoes must follow shower," 
the entry -1 under Column 2 says that ''shoes must follow drying,"’.... andthe 1 under Column 
11 says "shoes must precede 'dressed','' before the person can depart for the office. 

The matrix has certain other properties which we should explore now. First, note that 
all of the +1's are above the main diagonal (a line drawn from the upper left to the lower right), 
and all the -1's are below that line. This is as it should be when the precedence graph is both 
drawn and numbered correctly — that is, so that the numbers on the precedence graph are 
assigned to the items in such a way that they always increase from left to right along any 
precedence line. If the precedence graph is not drawn before the precedence matrix, it will 


Bsus 





MAEM RAR EA POSEN I" 


PROGRAMMING DYNAMIC WEAPON SYSTEMS 


TABLE 2 
Precedence for Matrix Dressing Problem 
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not be possible before hand to number the items in the convenient left to right sequence. Then 
the plus ones and the minus ones will be spread throughout the matrix with little pattern, 
except that they are symmetrically located about the main diagonal; that is, wherever a plus 
one appears above, a minus one will appear below. (Mathematically, this is called anti- 
symmetrical.) 

In order then to develop the precedence graph from a precedence matrix which is not 
already conveniently numbered, the following steps can be used. For illustrative purposes, we 
use the dressing problem but with some of the rows and columns interchanged, and we use 
letters to indicate rows and columns instead of numbers. In order to arrange the parts and 


operations in an abitrary, random order, renumber and list them alphabetically in the matrix 
as follows: 


(a) Coat (e) Shoes (i) Undershorts 
(b) Depart ({) Shower (j) Tie 

(c) Dry (g) Socks (k) Trousers 

(ad) Shirt (h) Undershirt 


Then construct the precedence matrix by taking each part or operation one at a time and 
noting for every other part or operation whether that operation is unordered, must precede, 
or must follow the other part or operation being considered. Enter as appropriate 0, +1, or 
-1, respectively, in the matrix. Table 3 is a matrix developed in accordance with this pattern. 
In order to assure that our evaluation of the precedence relations has been individually con- 
sistent, we can check by determining that the following conditions hold: 

1. If at the intersection of, say, the jth row with the jo column there is a plus one, 
then at the intersection of the ij row with the jth column there must be a minus one. 

2. If at the intersection of the jth row with the jth column there is a minus one, then 
at the intersection of the je row with the jth column, there must be a plus one. 
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TABLE 3 
Matrix Based on Table 2 
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If for every row and column, these conditions are satisfied, then the matrix is at least 
not inconsistent insofar as individual relations are concerned. We then can drop the part of 
the matrix below the main diagonal, and are ready to construct the precedence graph from the 
precedence matrix. 

The precedence graph is prepared from the precedence matrix as follows. Take the 
triangular matrix from above. Translate each precedence relation as stated in the matrix in 
graphic form. Start with the first row and work across noting for each column the precedence 
relation until the row is complete. Then proceed to the next row, etc., until all rows have been 
completed. This procedure is illustrated below. 

In developing the precedence graph from the precedence matrix, we will use a property 
of the "must precede" relation which is called "'transitivity.'' That is, for example, if "'r" 
must precede "'s"’ and "'s"' must precede "'t,"’ then, also, "'r'’ must precede "'t."" Because of the 
transitive property of the relation, we will find in the following steps for developing the preced- 
ence graph that some of the relations which have been graphed are already implied by an 
earlier statement. For example in the analysis for Row b, the first item is c - b, i.e., c must 
precede b. But, in the analysis for Row a we already had c - a and a - b. Hence, we have 
from the transitivity principle, c - a - b and thus c precedes both a and b. Hence, it is not 
new information when in the analysis of Row b, we find c - b and we can disregard it. 
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From the Matrix we determine the 
pairwise relations and state them 
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These pairwise relations are adjoined to 
form as much of the entire matrix as we 
can determine atany point inthe analysis. 
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Implies no new information 





7 Implies no new information 
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If we take the final graph, developed from the above process, and renumber as follows: 











We now have the same precedence graph (with addition of milestone 0) that we developed 
earlier for this problem. With this numbering we could develop the precedence matrix with all 
+1's above the main diagonal. 







Methods for preparing the precedence graph both by direct analysis of the program or 
by first preparing a precedence matrix have been presented here. Both methods if done cor- 
rectly give the same results. Which should be used depends upon the programmer making the 
study and which he finds better. In general the precedence matrix may be unnecessarily time 
consuming because it requires each and every relation be enumerated. Its advantage is that it 
is a systematic method of handling the often intricate relationships among milestones and has 
the virtue of assuring a complete analysis. On the other hand, if the precedence graph is 
developed directly from the program, it will save time because the transitive principle will 
automatically exclude the need for writing many of the redundant relations. Some program- 
mers have had difficulty in making the graph without aid of the matrix and hence we have 
developed both. 

Once the graph has been developed it should be checked for certain further inconsist- 
encies. For example, there should be no situations in which 
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XY 


either directly or indirectly. This would result in a problem of "which came first, the chicken 
or the egg?" Also, there should be no lines within the graph which do not terminate at a mile- 
stone, i.e., every job should have a purpose. When the precedence relations have been firmly 
established, it is possible to proceed to the next step. 


DETERMINING TIME AND COST OF ACTIVITIES 

Time values for activities are determined by conventional estimating procedures. How- 
ever, there are certain particular problems which arise. For example, if one first puts on his 
trousers and then his shirt, he will have to unbutton the trousers for tucking in the shirt tail. 
The method which we use for handling this interdependency is to neglect the effect of lack of 
independence up to approximately 5 percent variation in the activity time value. If the assembly 
operation for any part has a variation in its time value of more than, say 5 percent, because of 
differences in the amount of interference from other parts or because work must be undone 
then redone (as buttoning trousers twice), we arbitrarily exclude carrying out that activity at 
that greater time value by prescribing a "must precede" precedence relation which does not 
allow the activity to be carried out in this less efficient sequence. 

In the dressing problem for example, if it were determined that it takes more than 5 
percent more time to put on the shirt if done after putting on trousers than if before putting on 
trousers, we would include the precedence relation ''shirt must precede trousers."’ From a 
mathematical or computational point of view, this is an important assumption. It makes the 
computational problem much easier. However, it can be relaxed in shortening the over-all 
program cycle, and introduces additional costs. 

From the precedence graph it is noted that any milestone can be uniquely identified by a 
single number, say i=1,..., I and that each activity can be uniquely identified by the pair 
(i,j)i=2,...,1; j=1,...,1 where i is the milestone accomplished by the activity and j is 
the milestone before the activity. Each activity so defined requires a (more or less known, or 
estimated) amount of time to complete. This number can be designated oy aij and is the 
amount of time required to achieve milestone i after completing milestone j. It is assumed 
also that there is a cost ci associated with each value of aii. In general as as decreases in 
magnitude, Cij increases and that relationship is expressed in a table or with some suitable 


function 
P 1 
eyo t (at) 


e.g., { may be hyperbolic, negative exponential as earlier suggested, or other suitable function. 
The objective in programming is first to determine that path of activities, P*, for which 
. = max, (i, j)e P*. A path is a connected link running from the first milestone to the last, 


J 
e.g., in Figure 1 one such path is 


o> as 
i 


(1,0) (2,1) 6,2) 6,5) (7,6) (8, 7) (11, 8) 
the length of the path would be 


™ o*"o,.5*%,9 7,8 °"9.8* "9° 8 
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milestone i immediately follows milestone j. 
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The programming problem differs from the dressing problem, of course, and one such differ- 
ence is that in dressing (or any other activity carried out by one person only), the program's 
elapsed time is the sum of times of all activities. In w/s programs, however, each activity 
may be performed by a separate person or facility so that the elapsed time is the sum of times 
in the longest path. (This is modified, of course, when different programs are served by the 
same facility. In such cases, a queuing element is introduced, a feature of Boyan's target 
commitment scheduling but apparently not of Malcolm's Performance Evaluation Review 







In order to illustrate the problem of finding the longest path, assume a program pre- 
cedence graph as in Figure 3. The matrix in Figure 4 contains (a) 0's to indicate those activi- 
ties which do not exist or (b) a positive number in any element of a; ;) which indicates the 
amount of time required to achieve milestone i after milestone j and, also to indicate that 
















ip, i=l 2 3 4 5 6 
ar 
2 apy 0 
3 a3 0 0 
4 %, 0 0 0 
5 8 952 %3 54 9 
6 0 0 %3 9 +O 0 
7 0 0 0 0 as 0 
8 0 0 0 0 0 age 
9 0 0 0 ag4 0 agg 
10 0 0 0 a 
11 OO 0 0 0 o oO 
12 «0 0 0 0 oO 0 


Figure 4 - Matrix based on Figure 3 
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Using this matrix, the longest path can be found in many different ways, such as, tracing 
out all paths graphically and summing their activity times. The difficulty of this depends upon 
the number and complexity of the routes. A more formal method, suitable to larger programs 

) x and to machine computation, follows. 

Use two matrices fa; ;) and (b; j) as follows. aj) is the precedence matrix with 
activity times. (b; i) is computed in the process. Perform these computations on a; ;) by 
rows from top to bottom and from left to right. (a) Read the entry Bint in (a; ))- If Ajj = 0, 
enter in (b; i) at (i'j'); go to the next entry in @;,). If ay, > 0, go to row i" =j' in (b; :) and 
read largest Div: in that row. Compute Divs = inj + max Din . Enter sum Dis jn in ,) at 
(i'j'). (b) Go to next entry in a; ;) and repeat until all entries have been so operated upon. The 
length of longest path will be shown by the largest number in the last row or (b; i). The path 
can be printed out by printing out the pair (i*j*) in (b; ;) for the largest entry, bays Then, 
take the row i** = j* and print out the pair (i**, j**) with largest Di xj . Continue until first 
entry (1,1) is printed out. In case of ties, chose arbitrarily, or successively print out all 
equally long paths. 

The matrix (b; ;) for the assumed problem is shown in Figure 5. 

To make further explicit, this problem is numericised and solved, as shown in Figures 
6, 7, and 8. 

The longest paths are 


(12, 10), (10, 8), (8,6), 6,3), (3,1) = 8+9+7+6+43 = 33 
: (12, 10), (10, 7), (7,5), 6,2), (2,1) = 8+6+84+7+4 = 33 
(12, 10), (10, 7), (7,5), 6, 4), (4,1) = 8+6+8+5+6 = 33. 


The longest paths are used in programming for determining the elapsed time and for 
computing with the function 


c. =f (a) 


the increased cost of the program if some element(s) aij in the longest path are shortened so 
as to reduce elapsed program time. Of course, as the longest path is shortened, others 
emerge as the longest and these also may then be shortened as the program is compressed. 

In most practical situations, there is more than one program operative at any time and 
these are processed through the servicing facilities which perform the activities. Thus, there 
is also the problem of computing the expected holding times under different load conditions 
and of simulating the several programs to determine their over-all behavior as program 
variables are adjusted. Thus, the programmer is interested in determining expected and 
range of holding times for activities in the longest paths in each programming and determining 
the expected and distribution of likely outcomes for each program and for all programs as one 
over-all program is simulated under different time compressions, different priority schemes, 
servicing capacities, ''slack" allowances, etc. The program structures provide one basis for 
this evaluative procedure. 
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Figure 6 - The problem numericized 


@;,) = 6 0 0 6 0 0 0 





| 9 0 0 0 4 0 2 0 0 0 
10 0 0 0 0 0 0 6 9 7 0 
11 0 0 0 0 0 0 0 4 0 0 0 
12 0 0 0 0 0 0 0 0 0 8 3 0 


Figure 7 - The problem solved (ajj) 
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Figure 8 - The problem solved (bij) 
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LOGISTICS: CONDITIO SINE QUA NON FOR NATO DEFENSE* 


H. E. Ecclest 


INTRODUCTION 

The dominant position of logistics in the NATO Defense System is clearly indicated by 
the fact that of the 20 principal committees of the NATO council, nine of them deal specifically 
with logistics matters and five more deal with broad economic matters closely related to the 
industrial mobilization aspects of logistics. 

The nine specifically logistic committees are: Armaments; Infrastructure; NATO 
Pipeline; Military Budget; Planning Board for Ocean Shipping; Planning Board for European 
Inland Surface Transport; and Medical, Petroleum, and Manpower Planning. In addition, there 
are committees for: Economic Advisors; Annual Review; Civilian Budget; Industrial Planning; 
and Food and Agriculture. Thus, when one speaks in general terms of NATO logistics, one is 
speaking of the tangible heart of the North Atlantic Treaty Organization Defense System. 

Obviously, it is impracticable to place all the recognizable logistic activities under a 
single operating head, entitled "logistics.'' However, sound concepts of logistics and particu- 
larly, an understanding of how logistics is related to strategy, to tactics, and to economics, 
are essential for understanding NATO Defense. The present military situation contains many 
contradictions and paradoxes which create true dilemmas for the political and military leaders 
in NATO. To surmount, or solve, some of these problems, one must establish basic concepts, 
understand cause and effect relationships, have a philosophy, and describe the terms used. 

First, the problem should be examined both from the perspective of the military com- 
manders who are responsible for combat operations and from the perspective of the political 
leaders who are responsible for the Affairs of State. To do this, one must subordinate techni- 
cal detail in order to identify the central issues and to understand their relationship. 

Second, no commentator—no matter what his experience or other qualifications— should 
attempt to prescribe specific solutions to the problems presented. Solutions are the preroga- 
tive of those who are charged with responsibility for the conduct of affairs. The commentator, 
however, can analyze situations and point out principles and perhaps shed light on perplexing 
matters. 

For the purpose of this discussion, the following descriptive terms and related concepts 
are offered. 

Strategy is the comprehensive direction of power to control situations and areas in 
order to attain broad aims. This includes the selection and time-phasing of the minimum spe- 
cific objectives whose collective attainment will accomplish the broad aims. Strategy, being 
comprehensive, concerns both military and political leaders. 





*Manuscript received 3 March 1960. 
TRADM, USN, Retired. 
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Tactics is the immediate employment of forces and weapons to attain the specific 
objectives selected by strategy. Normally, it primarily concerns the military leaders; how- 
ever, with the advent of modern weapons, political considerations more and more influence and 
limit tactics. 

Logistics is the provision of the physical means by which power is exercised by organ- 
ized forces; it is the creation and sustained support of combat forces and weapons. The major 
components of logistics include infrastructure, transportation, ships, planes, weapons and 
equipment, supply, engineering, maintenance and repair, personnel and medical. The logistic 
process is both the economic element in a military system and the military element in an eco- 
nomic system; thus, logistics forms the bridge between the industrial-economic system of a 
nation and the actual operations of combat forces. Therefore, logistics inevitably becomes a 
major concern of both the political and the military leaders. In its economic or preparatory 
sense, it is largely civilian; in its operational sense, it is primarily military. 

In modern times, it is impracticable for the term ''Command" to have the same mean- 
ing as it once had in the person of an ''emperor general,'' such as Napoleon. ''Command," of 
necessity, has now become diffused through the whole national and international system. It is 
a blend of civilian and military, councils and committees, premiers and supreme commanders, 
and Commanders and Staffs. Nevertheless, it is important to say that Command has three 
major responsibilities: to create combat forces, to support combat forces, and to employ 
combat forces. 

From these concepts, several important corollaries and implications develop: 

1. Strategy and destruction are not synonymous. Strategy uses destruction when and as 
it is appropriate to the establishment or maintenance of control. 

2. Economic-industrial factors limit the creation of combat forces. 

3. Logistic factors limit the employment of combat forces. 


4. Command transforms war potential into combat power by its use and control of 
logistics. 


CONTRADICTION AND PARADOXES 

NATO is an attempt to build an effective military alliance without sacrifice of national 
sovereignty; however, the great scope and speed of modern conflict require a sacrifice of sov- 
ereignty to attain effectiveness in allied combat operations. While the maximum standardiza- 
tion of weapons and equipment has long been a goal of NATO, the improved economic conditions 
in these countries, by stimulating individual industries and national production, are creating a 
greater diversity of NATO weapons and equipment. This same economic improvement and the 
emotional pressure for greater expression of national sovereignty are drawing more nations 
into atomic weapon development at the same time as there is a greater recognition of the 
extreme dangers of widespread atomic or nuclear weapon capability. 

Logistic efficiency demands greater concentration and, in some instances, centraliza- 
tion. The threat of nuclear weapons demands greater dispersal and decentralization. As 
advanced technology designs and builds more sophisticated weapons, actual conflict is employ- 
ing more primitive weapons. For example, the nations of the world h-~e acquiesced in the use 
of assassination, terrorism, and mob violence in the various recent limited conflicts. Part of 
today's paradox lies in the conflict between (1) the need for rapid reaction and for great mili- 


tary flexibility and (2) the need for a new and broader concept of "command," as related to 
defense systems. 





CONDITIO SINE QUA NON FOR NATO DEFENSE 113 


While many persons have recognized the need for flexibility in modern military sys- 
tems, few have analyzed its sources. Flexibility requires a variety of forces and weapons that 
may be used in various combinations—appropriate to the situation, actually confronting the 
military commander. The forces and weapons he employs must be appropriate to the nature 
and degree of strategic control, required by his particular situation and by the further situa- 
tions that may result from his actions. Thus, it is obvious that the combat forces he employs 
must be versatile, both in equipment and personnel. It is essential that logistic support forces 
be equally versatile in their equipment and personnel. This rules out the "attractive fallacy" 
that short-term conscripts can provide effective military forces today. 

But even the best of equipment and men cannot be variously and flexibly employed, 
unless the tactical commander's logistic support is as responsive to his immediate needs as 
are his tactical units. This implies that the over-all military commander should be able to 
(1) allocate and reallocate major logistic resources, (2) exercise movement control over 
forces and supplies, and (3) control enough transportation of all types to insure 
immediate, positive action in critical areas. Does this mean that the military commander 
must have absolute control over all logistics? This of course would be absurd and impossible. 
It does mean that the problem of ''command control of logistics"' should be very carefully ana- 
lyzed and that various solutions should be tested by realistic war games, until the critical ele- 
ments and areas are clearly identified. Then the division of authority can be made with 
greater understanding and wisdom. 


THE NATURE OF THE LOGISTIC PROBLEM 

The objective of all logistic effort is to achieve maximum sustained combat effective- 
ness. Thus, the primary question is: What kind of combat action must be prepared for? If 
one considers that the NATO forces should act primarily as a temporary shield or "trip wire,"' 
then it is suitable to provide a relatively static, logistic support system. However, if one con- 
siders that the NATO forces may engage in protracted operation over wide areas, then a 
mobile, flexible logistic system is required. 

This distinction is particularly apparent, if one considers the requirements of limited- 
war actions and the difference between the type of operations SACEUR may conduct, as opposed 
to the type SACLANT may conduct. Since both SACEUR and SACLANT must operate under the 
same general NATO policies, the necessary compromises—between the military principle that 
a commander should control his own logistic support and the NATO policy that logistics is a 
national responsibility—pose many dilemmas. 

Some persons who recognize the need for more command control in logistics in cases 
of actual combat, are reluctant to grant this control to the NATO commanders, until the need 
is actually demonstrated by events. This attitude fails to take into account the realities of 
logistic planning and coordination. Vast amounts of information must be readily at hand, in 
usable form, before good planning and coordination can take place. During peacetime, it may 
be acceptable—even though inconvenient—for the NATO planners to go to their opposite num- 
bers on national staffs for logistic information, but in time of modern crisis, this is too slow 
Senior commanders, making urgent strategic and tactical decisions, must have up-to-date 
logistic information and planning factors immediately at hand. Their headquarters, facilities, 
and their staffs must be prepared in advance to process and manipulate these data with great 
speed; otherwise, valuable time will be wasted or wrong command decisions will be reached. 
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Good logistics alone will never accomplish a strategic concept. Bad logistics, however, 
may frustrate the most brilliant concept. Although there are many studies of supply systems 
and other technical aspects of logistics, there are few studies of the behavior of a logistic 
system. 

In his recent book, The Logistical Support of the Armies (Vol. II), Ruppenthal has made 
some very perceptive observations of the logistic system that supported the Allied Forces in 
Europe during 1944-45. We find on page 509: ''The movement of supplies from ship to Army 
depot, ... entailed a series of highly synchronized functions, the failure of any one of which 
could have a resonant effect reverberating along the entire line of communications."' Thus, 
combat forces which were created at great effort remained idle because faults reverberated 
through the logistic system. This quotation illustrates the two basic facts previously men- 
tioned: (1) "Economic factors limit the combat forces which can be created; logistic factors 
limit the combat forces which can be employed."' (2) ''Command transforms war potential into 
combat power by the manner in which it controls and uses the logistic process." 


x x x 


If one views the problem of NATO defense as that of creating a static, nuclear-armed, 
push-button combat force, which would (on signal) loose nuclear devastation on preselected 
targets from heavily protected missile sites, then logistics (in the full sense of the ideas 
herein expressed) is not too important. The problem simply becomes that of apportioning 
available funds and trained manpower among a relatively few weapons, in a relatively few con- 
crete emplacements, an elaborate if very vulnerable, electronic detection-and-command signal 
system, and a small screen of national superpolice. 

If, on the other hand, one considers that NATO defenses should include substantial 
naval, air, and ground forces (with the capability of a varied, flexible combat power) to be 
employed in an unpredictable and variable, unfolding combat situation, then one needs both 
staying power and mobility. The logistic control which command must have to employ such 
power effectively is entirely different from the logistic authority necessary to create and use 
the static, inflexible power of the first type. 

Admittedly, present political and economic attitudes and situations make the practical 
application of sound military concepts and principles exceedingly difficult. But, one should not 
ignore these principles and concepts, simply because they pose difficult problems; this is par- 
ticularly true if one expects a defense system to represent strategic realism. If the objective 
of retaining traditional concepts of sovereignty overrides the objective of collective military 
security, well and good; but in such event, one's security expectations should be correspond- 
ingly reduced. The two cannot be maximized simultaneously—one, or the other, must yield. 
This unhappy fact makes it mandatory that the objectives be thoroughly weighed and analyzed 
in the light of basic principles, rather than merely stated in general. 


NATO ACCOMPLISHMENT 
It is only natural that an observer, considering the complex problems confronting 
NATO, should be unhappy at its imperfections and eager to see its problems specifically 
solved. Also, it is easy to emphasize its shortcomings and to neglect its accomplishments. 
Two vital general facts are paramount: NATO has protected the security of its members. 


NATO has stimulated the economies of its members and has developed excellent methods of 
international cooperation. 

























































CONDITIO SINE QUA NON FOR NATO DEFENSE 


In the last 10 years, NATO’s specific accomplishments are significant, both in its 
mobilization base and in its operational capabilities. Among these improvements have been: 

Ec The NATO military commanders have been permitted to state their military require- 
ments to the member nations. They have been granted the right to prescribe the supply levels 
to be achieved and to receive periodic reports as to their actual status. 

The establishment of the NATO Planning Board for Ocean Shipping (PBOS) and the 
a Planning Board for European Inland Surface Transport (PBEIST), whose members can be 
| expected to exercise considerable authority in their own countries in times of emergency, will 
permit a high degree of coordination and efficient use of transportation. These two committees 
and the NATO pipeline committee, all reporting directly to the council, by their very existence 
belie the concepts either of NATO being a merely political group or of NATO forces merely 
7 being a “trip wire” for push-button thermonuclear retaliation. 

The NATO Maintenance Supply Services Agency, established in 1958 at Chateauroux, 
should not only reduce costs and increase logistic effectiveness, but it should stimulate the 
unification of the economies of Western Europe. 

The chief logistic accomplishment in NATO has been in infrastructure. This program 
is significant for (1) its obvious accomplishment in the construction of vital facilities and (2) its 
less obvious accomplishment in the psychological field and in developing planning, budgeting, 
and financing procedures. The major obvious accomplishments in the physical area are 160 
- airfields, 8800 km of POL pipeline, 2 million bbl of POL storage, and 4200 km of signal network 
il (either constructed or under construction). 

However, the program may be even more important for its other effects. Infrastruc- 
ture provides tangible proof of the successful operation of the alliance, thereby giving a valu- 
able psychological lift to all of NATO; all of this tends to stimulate greater integration of other 
effort and promote the vital sense of interdependence. Scarcely less important has been the 
development of budgeting and funding procedures which are gradually being extended to other 
areas of NATO work. Although the need for unanimity among 15 nations naturally slows plan- 
ning and construction, the budgetary controls (developed in the last 7 or 8 years) ensure that 
all expenditures are well screened and that the money is well spent. This, and the fact that all 
{ construction must be on the basis of international competitive bidding, builds confidence and 
stimulates the over-all economy. 

These great advantages, however, do not mean that all is clear and easy—quite the con- 
trary. There is always a conflict between the national interest and the NATO interest. 
Furthermore, there is always a problem of the influence of the commercial interests influenc- 
ing the writing of specifications. 

These interests conflict with NATO military commanders’ desire for (1) the best equip- 
ment available and (2) complete compatibility among equipments and systems. The fact that 
these various interests clash at every level—together with the principle of unanimity—make it 
easy to stay on dead center, to do nothing. This situation also severely tests the character and 
skill of the negotiators representing the various nations. 





CONC LUSION 
The most important and puzzling problem lies in the principle that logistics is a 

national responsibility. When one considers this problem in the light of the last 10 years' 

experience, one conclusion seems clear. 


H,. E, ECCLES 


The basic principle that logistics is a national responsibility is sound. However, the 
rigid and universal application of this principle is unsound. Variation can and must be made 
when military necessity so requires. Each exception to the principle should be judged on its 
own merits. The military commanders' increased authority in infrastructure, the development 
of NATO Maintenance Supply Services Agency, the development of the lightweight strike fighter, 
the maritime patrol aircraft, all represent sound departures from the rigid application of the 
principle. Perhaps the opposition to the idea of an integrated common stockpile would dimin- 
ish, provided the NATO nations could be assured that the national forces who draw from such a 
stockpile had already contributed their full share to the formation of the stockpile. 

Finally, the problems posed in this paper are neither simple, nor easy to solve. To 
deal with them effectively, we must evolve a new and comprehensive theory of modern human 
conflict—for it is only through a knowledge of fundamental principles and a knowledge of 
cause-and-effect relationships that wise action can be taken in the midst of the extraordinary 
dilemmas, paradoxes and contradictions, posed by the interplay of politics, strategy, econom- 
ics, logistics, and weapons of today. In all of this interplay, logistics provides the foundation 
for effective action. This action must deal with complex and unpredictable situations. The 
factor of logistic flexibility is vital; and, so in essence, logistic effectiveness becomes the 
conditio sine quo non for a realistic defense system. 








PROBLEMS 

It has been suggested that this journal might serve as a medium for publishing prob- 
lems in the areas of Logistics—with the idea that other persons might be interested in those 
problems and might submit comments. Readers are invited to submit brief statements on 
applied and theoretical problems in Logistics. Address letters to Managing Editor, Naval 
Research Logistics Quarterly, Office of Naval Research, Washington 25, D. C. 


AN ALTERNATIVE SOLUTION TO 
THE "LOST AT SEA" PROBLEM* 


Brian Gluss 


Armour Research Foundation of Illinois Institute of Technology 
Chicago, Illinois 





A problem posed by Bellman and considered by Isbell is as follows: 
Suppose one is a mile from a straight shore with no means whatsoever 
of ascertaining its direction. What is the optimum path to follow so as 
to (a) minimize the maximum distance travelled in reaching shore, 

(b) minimize the statistical expectation of the distance travelled, or 
(c) maximize the probability of reaching shore within a given distance 
travelled? Isbell found the solution to (a); in the present paper a 
sequence of approximations to the optimal policy for (b) is considered. 











INTRODUCTION 

A problem of considerable mathematical and geological interest is that posed and 
defined by Bellman as the "Lost at Sea'"' problem. Suppose it is required to trace a path from 
a point O to a straight line (or shore) in the same plane known to be at unit distance from O, 
no other information being available as to where the line is. The problem is to "optimize" this 
path so that the line is reached in the shortest possible distance travelled along the path. Three 
optimization criteria spring to mind: (a) minimization of the maximum distance travelled, 
(b) minimization of the statistical expectation of the distance travelled, and (c) maximization of 
the probability of reaching the line within a given distance travelled along the path. The solu- 
tion to (a) has been found by Isbell [1]; the object of the present paper is to present a sequence 
of approximations to (b) by finding suboptimal paths of two types. 





*Paper presented at the Sixteenth National Meeting of the Operations Research Society of 
America, Pasadena, California, November 1959. 


117 















The paths determined may be of use 
as approximations to optimal geological 
search paths. If a sample of earth contains 
elements of a material that is being sought, 
geologists are quite often able to give a 
rough estimate of how far away a layer of the 
required material is. These paths may then 
be used. Admittedly, the situation is differ- 
ent, in that information is acquired as one 

w proceeds along the path and takes new sam- 
ples; however, fragments of the material out- 
side the layer are oftenhighly discontinuous. 











Figure 1 - Type of path considered 


CONSTRAINTS FOR SUBOPTIMAL PATHS CONSIDERED 

The paths to be considered will be confined to those of the type shown in Figure 1. 
Assuming an arbitrary tangent AFW to the circle with center O and unit radius, the path 
comprises (i) a straight line from O to the tangent which it meets at A, (ii) the tangent AB 
that touches the circle at B, (iii) an arc BVX of the circle that leaves the circle tangentially 
at X, and finally (iv) a path XW, on or outside the circle at all points, that meets AFW at W. 

Such a path is suggested by the solution to the problem of minimizing the maximum 
distance [1], and given a path of type BVW it can be shown rigorously that OAB is optimal 
relative to it. 

The minimization problem may now be split into two independent parts: (1) minimiza- 
tion of OABV, and (2) minimization of VXW. This may be seen from Eq. (1). If d is the 
expected total distance travelled along the path, and d’ the expected distance travelled along 
VW given that the shore is not reached before V, then 
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where @ is the angle between OF and the radius from the shore. Hence d is minimized upon 
6 when @ = 33° 56’ , from which it follows that for this value of 6 


(2) d= 2. 8422+ a'/2. 


FIRST APPROXIMATION 
The first approximation considers paths XZW of the type shown in Figure 2. This 
comprises two straight lines XZ, ZW, where XZ is the tangent at X that meets Ox at Z 
(where XO Z = 5), and ZW is at an angle 8 to Ox. It is now required to determine those 
values of 8 and 6 that minimize d'. Note that the paths of this type cover a very wide range 
of paths, because of the two degrees of freedom allowed. Hence minimizing over f and 6 
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Figure 2 - First approximation 


should give a very good approximation to the best path; this is to some extent borne out by the 
fact that the second approximation reduces a’ only very slightly. Then 


25 
eee © | ™ 64 tan 2) da 
T 0 2 2 


(3) 
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Simplifying, the minimum is found with the aid of a digital computer to occur when 
6 = 50.3° and B = 77.0 , giving d' = 1.2746, and consequently, using Eq. (2), d = 3.4795. 
For Isbell's solution, the expected distance is 3.6276, which is approximately 4.2-percent 
greater. 

This suboptimal path, which appears in Figure 2, gives an approximate idea of what the 
optimal path really is: this may be seen by smoothing out this path. This suggests the second 
approximation. 


SECOND APPROXIMATION 
Assume now that the path XW is, as shown in Figure 3, of the form 


2 
(4) X = a) + a,y + acy 
and that the curve leaves the circle at a tangent. Then 


(5) a, - 2a, = 0, 


Since it is obvious that an optimal curve should be perpendicular to FW at W. 





E (d) = 3.4691 





Figure 3 - Second approximation 


ee ee 
(6) ay + a; COS Y + ay cos’ Y = siny , 
since (sin Y , cos Y) is on the circle. Also, because of tangency at Y , 
(7) a, + 2a, cosY = -cotY. 


It is easily seen that 
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(8) ce oe J s(y,v) da, 
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where s(y', y) is the distance travelled along the path to the point at which y = y . Further- 
more, 


cos Y 
(9) sly, v)=7 + f. V1 + [h(1 + y) |? dy , 
y 





where h = ~ 
1 + cos Y 


From Eqs. (8) and (9), it is now required to find the value of Y that minimizes a’. 
Using a computer, it is found that this value of Y = 33° 58" , and hence that d' = 1.2538 in 
this case, when it follows that d = 3.4691. This is an improvement of approximately 0.3 per- 
cent over the best path of the first type. Since Y = 33° 58' , it follows from Eq. (3) that the 
best path of the second type is 


(10) x = 1.5108 - .8115y - .4057y”. 
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Note that coincidentally 6 and Y are almost identical for this suboptimal path, and 
hence that the path involves, after reaching the circle at B, going almost exactly half-way 
round the circle before leaving it at X. 
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[1] J. R. Isbell, ''An Optimal Search Pattern," Naval Research Logistics Quarterly, December 
1957. 





Professor Klein has reported an error in his paper published in the Quarterly, Vol. 4, 
No. 4, pages 269-286. Algorithm 4, page 284 and the remarks pertaining to it on page 285 are 
incorrect. An error in the proof (which was not included in the paper) was recently uncovered 
and a counter example was provided by Mr. A. Veinott, Jr. of Columbia University. 


“The Naval Supply Depot, Philadelphia, is endeavoring to fill two Operations Research 
Analyst positions in Grades GS-12 and GS-13. These are interesting and challenging positions 
involving the direction and coordination of logistics research and development projects in the 
Military Industrial Supply Agency. Interested persons may apply by letter or in person to: 


Consolidated Industrial Relations Office 
Naval Supply Depot 

700 Robbins Avenue 

Philadelphia 11, Pennsylvania” 


“The board of U. S. Civil Service Examiners Aberdeen Proving Ground, Maryland 
announces a Regional Examination for appointment as a Mathematician, GS-9 through GS-15.” 


MATHEMATICAL MODEL FOR DECISION MAKING IN A CRISIS SITUATION 


Where 


x =v? Bi + sr + VT,s 
Pp 


Rank of commander whose equipment is deadlined 


Temperature at Clearfield, Utah when "gig" is received (for east coast stock 
points substitute temperature-humidity index). 





NOTES 


Degree of sarcasm in personal letter sent to depot commander by old classmate 
whose equipment is deadlined. 


Impact on the supply system if missing item interferes with a Varsity athletic 
program at requisitioning activity. 


Time (in hours) since supply officer was last ''chewed out" by his commander for 
a "stock-out." 


Pride of stock clerk in never having a "'stock-out." 


Probability of doing the wrong thing 
or 


Holmes Law: The probability of not having the right part in the right place at the right time is 
equal to the sum of the squares working on the problem. 


LCOL Fenwick W. Holmes, USMC. 


EDITOR’S NOTE: This article may not meet all the requirements usually necessary for publi- 
cation in the Naval Research Logistics Quarterly. However, it does make one become conscious 
of the fact that there are many causes and effects to be considered in the process of providing 
Military Logistics support, although Holmes Law may not accurately reflect them. 
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NEWS AND MEMORANDA 


Readers are invited to submit to the Managing editor items of general 
interest in the field of logistics 


SUMMARY 
PROCEEDINGS OF THE 1959 COMPUTER APPLICATIONS SYMPOSIUM 


Armour Research Foundation of Illinois Institute of Technology 


The Proceedings of the 1959 Computer Applications Symposium contain the 14 invited 
papers presented at a 2-day meeting in Chicago, on October 28 and 29, 1959, and the edited 
text of the associated discussions. The first day's papers were devoted to business and 
management applications, and those on the following day stressed engineering and scien- 
tific applications. 

Several themes were each represented by a group of papers. One of those was the 
development of methodology of general interest to exploit the capabilities of specific machines, 
ranging from large to small. Another was the extension of a general-purpose computer's 
resources by means of specialized input-output and display equipment to handle graphical 
information. A third was the impact of electronic data processing on the handling of massive 
clerical procedures such as inventory transactions. The largest group of papers, and a con- 
siderable portion of the panel discussions, were addressed to the question: Now that we have 
the machines, how do we communicate with them? The answers centered around the pros and 
cons of automatic programming and language design. Finally, there was a lively report on the 
current status of computer technology in the USSR. 

Titles and authors of papers presented are as follows: 


Business and Management Applications 


Shareholder Record-Handling with the Aid of Character-Recognition Equipment, J. M. Wells 
Around the World in Eighty Columns, W. E. Hanna, Jr. 

Cost Reduction Through Integrated Data-Processing, R. F. Hamaker 

Some Aspects of Computer Technology in the USSR, S. N. Alexander 

Experience and Plans for Marketing-Research Operations, R. B. Wilson 

A Modern Approach to Inventory Control Utilizing a Large-Scale EDPM, W. F. Harvey, Jr. 


Current Developments in Common-Language Programming for Business Data Systems, 
E. J. Albertson (155 + x pp., $3.00). 
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Engineering and Scientific Applications 


Linear Programming on the Bendix G-15 Computer, J. R. Wall 


The Design and Use of the APT Language for Automatic Programming of Numerically Con- 
trolled Machine Tools, D. T. Ross 


A Quasi-Simplex Method for Designing Suboptimum Packages of Electronic Building Blocks, 
R. H. Glaser 


The International Algebraic Language and the Future of Programming, C. Katz 


Training for Engineering and Scientific Applications via Compilers, Interpreters, and Assem- 
blers, W. F. Atchison 


Scientific Design Procedures Utilizing a Small Computer, T. I. Harris 


FORTRAN Experience and Remote Operation by Non-computer Specialists, F. Engel, Jr. 


TIMS MEETS IN BRUSSELS IN 1961 


The 8th Annual International Meeting of The Institute of Management Sciences will be 
held in the Palais des Congrés at Brussels, Belgium, 23-26 August 1961. The following topics 
will be included in the program: 

1. Programming under Uncertainty 

. Simulation and Gaming 

. Rational Investment Decisions 

. Organization Theory and Analysis 

. Subjective and Objective Probability 

. Reality and Theory in Management Science 
. Adaptive Systems 

8. Behavioural Sciences. 

Arrangements will be made for meetings of the colleges, and special sessions are also being 
planned, including one which will be devoted to reports on research from the newly founded 
International Center for Management Science at Rotterdam, The Netherlands. 

The registration fee is $12 for TIMS members and $15 for others. The registration 
will be closed on April 1, 1961. Prospective participants who fail to apply before that date are 
advised to write to Dr. Jacques Dréze, Chairman of the Local Arrangements Committee. 
Address: Department of Economics, University of Louvain, Louvain, Belgium. 


INTERNATIONAL STATISTICAL INSTITUTE 


The 33rd Session of the International Statistical Institute will be held in Paris from 29 
August to 7 September 1961. 





RECENT PUBLICATIONS 


EDITOR’S NOTE: Listing of a publication in this section is for record and reference only 
and does not constitute an endorsement of point of view or advocacy of use. 


TIME SERIES ANALYSIS, by E. J. Hannan. Methuen, London: Wiley, New York. 
147 pp. $3.50. 

One of Methuen's Monographs on Applied Probability and Statistics, this book provides 
a broad and current analysis of discrete stationary time series. Developed from the author's 
notes for a graduate statistics course, it well deserves to become a standard text at this level. 
Of course, it cannot attempt to be complete or detailed in 147 pages (it doesn't mention multi- 
ple time series), but it does give a very readable survey. References to recent work make it 
valuable not only for introduction but also for getting caught up. 

Chapter 1 gives the probabilistic background and spectral theory. For those unfamiliar 
with infinite dimensional Hilbert spaces, the author introduces them via a finite dimensional, 
circular process, which greatly clarifies both spectral decomposition and linear prediction. 
Chapter 2 discusses the estimation of the correlogram and of the parameters of autoregressive 
schemes. Chapter 3 gives the estimation of the spectral density (spectrum). It explains why 
the periodogram is inadequate and how, by a suitably weighted average of it, the useful esti- 
mates of Bartlett, Tukey, etc., are obtained. Chapter 4 gives tests for exact periodicities 
(jumps in the spectral distribution), tests of the length of autoregressive schemes, and more 
general "'goodness of fit'' tests of the spectrum. Finally, confidence limits of the spectrum are 
outlined. Whereas the first four chapters deal with the random part of a time series, Chapter 5 
describes how to extract the "deterministic" part (trend) by regression when the residuals 
(error, random part) are correlated (this correlation is the essence of time series analysis, of 
course). This chapter is nicely handled with matrix notation. 

There are a few points which this reviewer disagrees with, or would qualify. For 
example, in Chapter 3, page 74, it is suggested that in estimating a spectrum having only one 
relative maximum, "one will almost certainly do better by fitting some such finite parameter 
scheme as an autoregression. For this reason it may be simpler to plot the correlogram first. 
If this appears to have a very simple structure one may fit an autoregression to begin with and 
subsequently examine the spectral density of the residuals (after regression), thus avoiding a 
considerable amount of computation.'' Now the correlogram is often useless in trying to deter- 
mine accurately the type or order of the scheme, even if it is simple. Furthermore, once the 
correlogram is calculated, it requires only a little more extra computation to obtain the spec- 
trum, and this is usually automatically done as part of the computer program. 

Another example of qualification is the test for periodicities in Chapter 4. The author 
doesn't sufficiently emphasize that the only way to differentiate between an exact periodicity 
and a sharp peak in the spectrum is by an infinite amount of data or by a priori information, 


for example in communications engineering or astronomy, where a known signal may or may 
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not be present. Scientists generally lack such a priori information and should be discouraged 
from an undue effort to determine exact cycles. 

The above comments are not intended to obscure the fact that this is overall an excel- 
lent little book, at a price that should make it popular. 


(Reviewed by Thomas Wonnacott) 


LES CHOIX ECONOMIQUES, by P. Rosenstiehl and A. Ghouila-Houri. 

The processes of the real world are dynamic in nature—a fact that has gradually pene- 
trated the ranks of mathematical model makers of economic, industrial, and military problems. 
Static or steady-state processes, significant and useful as they are, must always be regarded 
as limiting cases of the more accurate descriptions involving time variation. Furthermore, 
they must be interpreted in these terms. In return, it is often the case that they can more 
profitably be treated, both analytically and computationally, from this point of view. 

In considering multistage decision processes, a number of mathematical approaches 
are currently available. Among these are the general theory of stochastic processes, the 
theory of dynamic programming, and the many varieties of simulation or ''gaming"' techniques. 

The book by Rosenstiehl and Ghouila-Houri, with the aid of a number of collaborators, 
is an excellent introduction to the modern field of sequential decision processes. The text is 
interesting and well-written, with copious references to contemporary work. All those inter- 
ested in the areas of mathematical economics, operations research, systems analysis, and 
applied mathematics in general, will enjoy reading this book, and derive some benefits from it 
in several different ways. 

The first chapter appears in both French and English, and the succeeding chapters are 
in one language er the other, with brief summaries in the alternate tongue. The second chapter 
discusses dynamic programming; the third presents the rudiments of the theory of stochastic 
processes with applications. 

The fourth through eighth chapters are devoted to the theory and practice of simulation. 
There are discussions of traffic control (both vehicular and air), logistics problems, and ques- 
tions of river control. Chapter Nine deals with war gaming and Chapter Ten with business 
gaming, of the type introduced by the model constructed for the American Management Associ- 
ation by Bellman, Clark, Craft, Malcolm, and Ricciardi. 

This brief description of the contents does not at all indicate the careful, thoughtful, 
and erudite presentation of the material in the various chapters. It is a pleasure indeed to see 
books of this calibre appearing in this new field. 

(Reviewed by Richard Bellman) 
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The NAVAL RESEARCH LOGISTICS QUARTERLY is devoted to 
the dissemination of scientific information in logistics and will publish 
research and expository papers, including those in certain areas of 
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to improve the efficiency and effectiveness of logistics operations, 
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